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Abstract 



Results on the correlations of low density classical and quantum Coulomb systems 
at equilibrium in three dimensions are reviewed. The exponential decay of par- 
ticle correlations in the classical Coulomb system - Debye-Hiickel screening - is 
compared and contrasted with the quantum case where strong arguments are pre- 
sented for the absence of exponential screening. Results and techniques for detailed 
calculations that determine the asymptotic decay of correlations for quantum sys- 
tems are discussed. Theorems on the existence of molecules in the Saha regime are 
reviewed. Finally, new combinatoric formulas for the coefficients of Mayer expan- 
sions are presented and their role in proofs of results on Debye-Hiickel screening is 
discussed. 
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I. COULOMB SYSTEMS 



A. Introduction 

The emergence of thermodynamics and of various physical laws from the more funda- 
mental levels of atomic theory and statistical mechanics were high points in our education 
in physics. But by the time we reached textbook treatments of the complex world of the 
Coulomb interaction, for example, formation of atoms and molecules and equilibria between 
them, we learned to demand less from theory and to be more content with reasoning by anal- 
ogy consistent with thermodynamics. We know experimentally that atoms and molecules 
form, that they can behave like mixtures of ideal gases, so phenomenologically they have 
chemical potentials and we no longer insist that this emerge in some limit from A^-body 
quantum Coulomb systems. We also know that Coulomb systems can have screening and 
dielectric phases. But when screening is expected to occur, it is quite a common practice to 
rather bluntly replace the Coulomb potential by a screened potential inherited from mean- 
field theories without inquiring too deeply into the legitimacy of this change. Likewise, there 
does not exist a first principle theory of the dielectric constant that does not presuppose the 
existence of atoms and molecules. At best, assuming that atoms form, one uses several laws, 
such as the Clausius-Mosotti formula, whose microscopic foundations have to be elucidated. 
In all these cases, this is a sensible attitude since most of non-relativistic physics is reputed 
to be hidden within the A^-particle Coulomb Hamiltonian. There are however basic facts 
that can be cleanly formulated as limiting theories and shown to persist near the limit as 
well. 

This review is devoted to these types of results on Coulomb systems at low density. 
At low density the most famous properties of Coulomb systems are related to screening. 
The classical Debye-Hiickel theory and its quantum analogue the random phase approxima- 
tion have been supplemented in recent years by additional information established both by 
rigorous proofs and by closer inspection of resummed perturbation theory consistent with 
aspects of the rigorous proofs. There is also a precise understanding how a dilute assembly 
of quantum nuclei and electrons can form gases of atoms and molecules in the Saha regime. 
These results are physically well understood but the challenge to derive them has shown 
that standard perturbation theories miss some effects; indeed it turns out that quantum 
Coulomb systems do not screen, in the strong sense that the classical system does. This is 
a principal focus of this review. 

The first rigorous results were at the level of thermodynamics, the celebrated stability 
of matter theorem of Dyson and Lenard and the existence of thermodynamic behavior 
established by Lieb and Lebowitz. They hold for all states of nonrelativistic matter. The 
importance of the stability result is that it shows that quantum mechanics and electrostatics 
is a consistent simplified statistical mechanics of nature that needs no intervention from more 
fundamental levels of description. 

The next level, which is the focus of this review, is to ask for information on correlations: 
here the variety of physical phenomena in question is so vast that we cannot hope for too 
much generality: thorough consideration of the low density behavior is the first and the 
most natural step in this direction. 



I-l 



There are three parts. The first is devoted to classical systems. Here a significant achieve- 
ment is the proof that the Debye-Hiickel theory is a limiting theory and Debye screening 
holds close to this limit, that is, at sufficiently low density (or high temperature). This is 
the subject of chapter II. 

Frequently proofs of such results are very long with tedious parts devoid of physical 
insight, but there are also components where a physical insight becomes a precise argument 
or where a useful new identity such as a resummation of perturbation theory appears. The 
philosophy here and throughout the review is to state the theorems and to give only those 
ideas in the proofs that conform to this specification. Many of these arguments are founded 
on the Sine-Gordon transformation, described in chapter II. We hope that our applications 
will raise appreciation for this transformation. We hope to display enough arguments without 
overwhelming technicalities to satisfy the physicist curious to see the way between physical 
concepts and their full mathematical realization. 

In the second part we consider quantum mechanics. In between there is a natural transi- 
tion chapter (chapter III) which revolves around the absence of screening in classical dipole 
systems. We present the simple proof of this lack of screening. This chapter prepares for 
the study of the quantum case because it turns out that the Feynman-Kac formula reduces 
quantum equilibrium statistical mechanics to a classical formalism in which both monopoles 
and higher order multipoles appear. 

The proper quantum mechanical effects on screening are exhibited in chapter IV in the 
simplified framework of the semi-classical approximation. The well known Feynman-Kac 
representation of quantum statistical mechanics is recalled and quantum statistics are ne- 
glected. Perturbing around the classical gas makes immediately clear why intrinsic quantum 
fiuctuations always destroy Debye screening. This appears to have been first suggested by 
Federbush, see (Brydges & Federbush 1981, page 428, section 5). One also discovers that 
the remaining screening mechanisms are more subtle than in the classical revealed 
by the different long distance behaviors of various types of correlations. This chapter is 
intended to be a gentle baptism before the confrontation with the full quantum mechanical 
setting in chapter V. 

Here the mathematics is much harder and less has been proved: in chapters V and 
VI the style is that of formal perturbation theory. In these chapters, following mainly the 
work of Alastuey, Cornu and collaborators, one exploits the full power of the Feynman-Kac 
representation including the quantum statistics. In this language, the quantum Coulomb gas 
is similar to a certain classical gas of fluctuating multipoles and all the efficient techniques of 
classical statistical mechanics are at hand. In particular, one can perform the usual partial 
resummations needed to cure the Coulomb divergences and end up with well defined Mayer- 
like diagrammatic rules. This offers a new perturbative scheme for the many-body problem 
which is particularly suited for low density expansions, in contrast to the standard Feynman 
perturbation theory with respect to the coupling constant. 

The formalism is applied in chapter VI to the determination of the long distance behavior 
of the correlations and to the equation of state at low density. The (formally) exact results 
presented here provide the most detailed information available now on these questions. Some 
of them are new, in particular the explicit formulae for the tails of the correlations at low 
density. Also one recovers and completes the equation of state already obtained by Ebeling 
and coworkers by different methods. 
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In chapter VII we report on a nice development, initiated by Fefferman, on atomic 
and molecular phases at low density (the Saha regime). One considers a joint limit, the 
atomic or molecular limit, of vanishing density and vanishing temperature by fixing the 
chemical potentials (negatives) and letting the temperature tend to zero. Lowering the 
temperature enhances the probabilities for quantum mechanical binding whereas lowering 
the density favors ionization. It turns out that in this limit each value of the chemical 
potentials determines in principle certain chemical species and the system behaves as a 
mixture of these ideal substances. Ionization equilibrium phases are obtained in this way. 
In this asymptotic sense, it gives a precise meaning to the notion of atoms and molecules in 
the many-body problem. 

The third part consists of a single chapter VIII that is more technical: rigorous results 
at low density require control of the convergence of various types of cluster expansions. The 
proof of Debye screening has stimulated the development of resummation techniques which 
may be useful in other situations. For example there are new combinatoric formulae for 
Mayer coefficients. Without aiming at completeness, this last chapter provides the dedicated 
reader with some arguments left out in the main course of the text and gives an updated 
summary of some of the available tools in this field. 

The considerations in this review are essentially concerned with our three dimensional 
world. There has of course been dramatic progress in the one and two dimensional solvable 
models of Coulomb systems, which we do not include in the review, apart from isolated 
comments included to demonstrate some contrast with three dimensions. Let us also mention 
that basic new results on the stability of Coulombic matter in presence of magnetic fields 
and quantized radiation field have been obtained recently (references in I.C). 

The history of proofs of fundamental results, for example that statistical mechanics 
predicts phase transitions, suggests that difficult results can simplify over time and serve 
as consistency checks on what comes next. Some but not all of the derivations in this 
review arc far from this stage: nevertheless, the existence of the subtle tunneling corrections 
in the Dcbyc asymptotic regime revealed by the Sine-Gordon representation is not seen 
by conventional Mayer graph summations. The suggestive view of the quantum gas as an 
assembly of random multipoles displayed by the Feynman-Kac representation and the related 
lack of screening are much easier to miss in the conventional many-body perturbation theory. 
The physically very natural molecular limit considered in the main theorem of chapter VII 
should be introduced and discussed in any text book on quantum statistical mechanics. 

We conclude this introduction with a list of some questions that seem to us of interest, 
to be pursued in the spirit of this review and perhaps by similar techniques. 

CLASSICAL: Debye regime 

A more detailed understanding of boundary conditions and charge expulsion. 
Results on correlations along walls, surface charges. 

Understanding the improved Debye-Hiickel theories of Fisher et al. as limiting theories. 
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Dipole systems 

More definitive results on the dielectric constant and clarification of the literature in 
the light of such results. 

QUANTUM: Low density 

Prove absence of exponential screening or at least finiteness of each individual Mayer- 
like diagram of chapter V. 

Correlations in non-uniform systems, e.g. along walls. 

Coherent treatment of para and diamagnetism, boundary currents. 

Slightly imperfect molecular gases 

Lowest order corrections to the molecular limit, atomic correlations and van der Waals 
forces. 

Possible existence of dielectric phases with true atomic dipoles. 
Compatibility of macroscopic electrostatics and statistical mechanics of quantum charges. 



B. The classical Coulomb gas 

We consider several species, a = 1, . . . , 5, of charged particles with charges e„ of both 
signs, confined to a box A in M:'. If there are N particles in total sitting in an external 
electrostatic potential 4> then the potential energy is given by 

N 

Cqi) • • • ) ^Ni Cajv) + ■ ■ • + ^ ] CQ.0(rj) (1-1) 

1=1 

where 

N 

[/(ri,ea,,...,rjv,e„jv) = ^ ea,e«.1/(ri - r^) (1.2) 

i<j=l 

V is the two body Coulomb potential energy 

= (1-3) 

and the eUipsis indicates potentials due to other forces, that are discussed below. 

Suppose that denotes the number of particles of species a. Then the equilibrium 
statistical mechanics is determined by the partition function 

oo cS / njVq r f ( ^ \ 

Sa = Sa(/30) = J] n / rfri... / rfr^exp -/3C/-/3^e„,0(r,) (L4) 

{jV„}a=l \ i=\ ) 

and Za is called the activity of species a. In particular the pressure P is given by 
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/3P= lim ^InSA. (1.5) 

lAHoo |A| ^ ' 

The activities in the partition function are related to the chemical potential /Xq of species a 
by the standard formula 

^'^ ^ (27rA2)3/2 ^^-^^ 

where is the de Broghe length given below in (1.7). This formula results from the classical 
limit of quantum mechanics. 

In three dimensions a purely Coulomb Boltzmann factor exp(— /3C/) is not integrable 
because of the Coulomb singularities at = r^: additional forces are necessary for stability. 
We suppose there is a length scale A where the Coulomb potential gives way to these other 
forces. We used A because it suggests a de Broglie length 

A, = h{i5lm^fl'' (1.7) 

where is the mass of species a. For example, for individual point particles, classical 
mechanics gives way to quantum mechanics when their wave functions overlap, which occurs 
at separation 0[\a)- For ions, A would instead be the radius of an outer orbital. For 
Fermions, one attempts to model the Pauli principle together with electrostatic repulsion 
by classical effective short range forces. We will either represent these short range forces by 
a hard core of radius A or we will smooth the Coulombic singularity by replacing |r|~^ by 

Ko,A(r) = |rri(l-exp(-|r|/A)) (1.8) 

Neither choice is likely to be a realistic description of the physics at this length scale, so 
these are only claimed to be useful when one can show that there are predictions which are 
not very dependent on how the short range force is chosen. 
The notation Voo,\ conveys the equation 

Voo,x ^Voo-Vx (1.9) 

where we set 

Vyiv) = |r|-iexp(-|r|/A); ^^(r) = l/|r| (1.10) 

The potential V^^x is the Green's function for (— A+A^ A^) /(47r) with the boundary condition 
appropriate for a system confined to a box A with insulating walls. We shall also consider 
systems confined in a box with conducting grounded walls. In these cases VA(r— r') is replaced 
by Vx{r,r') which then denotes the Greens' function with zero boundary conditions. 

Consider the case where the eUipsis in (1.1) represent a hard core of radius A. By 
Newton's theorem each point charge can then be spread out into an equivalent charge on a 
sphere of radius A without changing the Coulomb interaction energy. If self-energies of the 
spheres are added into the Coulomb energy then it is positive because it can be written as 
the integral of the square of the electric field. Since the self-energy of N spheres is | e^/ A 
this argument shows that 
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Stability: U > -BN; B ^ 0{-^) (1.11) 

Thus the introduction of the hard cores has made Sa finite (Onsager 1939). 

The same stabihty bound holds in the other case where the additional terms in the ellipsis 
in (1.1) smooth out the Coulombic singularity. The total potential energy with self-energies 
added is 

2^^"^^". |,^_,.| (1-12) 
which is positive because it equals, in Fourier space, 

The self energies are again ^^^e^/A so we obtain (1.11). Analogous arguments using the 
appropriate eigenfunctions in place of exp(ik • r) show that stability bounds hold also for 
particles in an grounded container. 
We use the summation identity 

{Na} ot " Ar=0 ai,...,ajv 

to write the partition function (1.4) in the form 

N=0 ' k=l \ i=l / 

where £ = (r, a) unites the spatial and species coordinate so that (f){S) — ea(f){r) and the 
dS integration means 

s 

d£--- = Y^ / dr--- (1.16) 
a=i -^^ 

and z{S) = z^- 

We denote by a; = (^i, . . . ,Sn) a configuration of the particles, U{uj) = U{Si, . . . ,Sn) 
and z{u!) — z{Si) • • ■z{£n)- The grand canonical summation over these configurations is 



dw 



oo ^ „ N 



and the grand canonical average is, for zero external potential, 

(...) = ^ I dujz{uj)e-^^^^^ ■ ■ ■ (1.18) 
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A system in which the species occur in pairs of equal but opposite charge with equal 

activities is said to be charge-symmetric. 

The observable that measures density of particles of species a at r is given by 

TV N 

p{S,U;) = =J2^a,aA^ - Tj). (1.19) 

i=i i=i 

and n - particle distributions for points r^, a = 1, . . . ,n are defined by 

n 

p{S,,...,Sn) = {l[p{Sa)) (1.20) 

a=l 

These definitions coincide with the standard definitions of distribution functions in statis- 
tical mechanics when Si, . . . are distinct. If they are not all distinct (coincident points) 
then our definition inchides squares or higher powers of delta functions. At various points 
in the text we will use the phrase "does not include contributions from coincident points" to 
indicate that such powers of delta functions are omitted. To illustrate what "including con- 
tributions from coincident points" means, consider J dS f{£)p{£,£i). The integral smooths 
the delta function in (1.19) so that the diagonal terms 6{£, £j)6{£i, £j) in p{£, u)p{£i, u) 
give the "coincident contribution" /(^j)'^(^i' ^j)- 
We also introduce the charge density observable 

N 

c{r,u;) = '^eapa{^,a,u;) = ^ea/(r-rj). (1.21) 

a j=l 

and the charge distributions 

c(ri,...,r„) = (nc(ra)) (1.22) 

a 

More generally one obtains mixed charge - particle distributions by taking the average of 
products of density observables and charge observables (1.19) and (1-21). 

We shall also consider the corresponding truncated distributions (correlations) pr and 
ct defined in the usual way 

Pt{£i: £2) = P{£l. £2) - P{£l)p{^2) (1.23) 

and so on. 

Correlation functions of charge densities can also be obtained as the response to an 
infinitesimal external potential by variational differentiation of the logarithm of the partition 
function Sa(/90) with respect to — In particular, the charge density induced in the system 
by a potential perturbation at ri is 

and the truncated distribution of two charge densities is 
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S d 

By this definition, we find that 

cr(ri, ra) = (c(ri)c(r2)) - (c(ri))(c(r2)). (1.26) 

The physical properties of an homogeneous infinitely extended charged fluid are conve- 
niently described by its linear response to an external classical charge density c^^^{r). One 
defines the static susceptibility x{^2 — fi) as the linear response at r2 of the charge density 
of the fluid to the external density at ri 



X(r2 - ri) 



c r2 



(1.27) 



The function x(r) describes the shape of the screening cloud around an infinitesimal point 
test charge at the origin. In Fourier space, 

X(k) = J x(r)e-*^dr (1.28) 

is related to the Fourier transform S(k) of the charge-charge correlations CT(r, 0) (1.26) at 
zero potential (the structure factor of the fluid) by 

X(k) = -^S{k) (1.29) 

Finally, the total effective potential due to the test charge plus its screening cloud is 

V^^^(k) = (1 + x(k))V^(k) = e-'{k)V{k) (1.30) 

where we have defined the static dielectric function £(k) by £~^(k) = x(k) -|- 1. Equiva- 
lently we immerse two test charges (ja, ea)a=i,2 in the ensemble. Their (nonlinear) effective 
potential energy is the excess free energy due to the external charge distribution given by 

exp(-/3y-*f((r„ej,=i,2) 

_ (exp {-:3ciV * (jri) - Jc-A' * ^(rj) - :^( ic-2^'(ji - r-,))) 
{ex.p{-(3eiV * c(ri)))(exp(-/3e2V' * c(r2))) 

But then we take the test charges infinitesimal and find the leading term 

V^''{{ra,e,)a=i,2) ~ eie2y^^(ri -r2) 

V'^ir, - r2) = V{r^ - r^) - pV * ct * V{r, - r^) (1.32) 

which is equivalent to our previous formula (1.30). 

A well known criterion for the plasma phase of the Coulomb system is the vanishing 
of e-^{k) as k ^ 0. This implies with (1.29) that S{k) ~ (47r/5)-i|kp, k ^ 0, which 
is equivalent with the second moment rule (Stillinger-Lovett rule) for the charge-charge 
correlation 



;i.3i) 
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rfr|r|2cHr,0) = -^ (1.33) 
If the Coulomb interaction U is set to zero the partition function is 

.exp(/...(.,e— ). (1.34, 

We will call the resulting system an ideal gas, because the external potential is just a multi- 
plication of the activities by exp(— /?eQ,j0(ri)). In this case all distributions are computable, 
e.g., 

(c(ri))ideaw = S{-f3<j){v,)) "ideai,A(/90) = ^ e^.z^.e-^'^^'f'^'''^ (1.35) 
where we used (1.24). 



C. The quantum Coulomb gas 

We consider S species of quantum point charges (electrons and nuclei) in a box A in R^, 
with masses m^, charges Ca and spins Sq, a — 1, ... ,5. Each species obeys Fermi or Bose 
statistics and at least one species is Fermionic. The system is governed by the non-relativistic 
N-particle Hamiltonian 

Hk,n = - X] ^ A, + ^ eaie«.1/(rj - r^) + ^ ea,0(rj) (1.36) 

i=l °' i<j i=l 

acting on the A'"-particle Hilbert space appropriately symmetrized according to the statistics 
of each species. The Laplacians Aj in (1.36) have Dirichlet boundary conditions on the 
boundary of A. The external electrostatic potential is set to zero unless we warn the reader 
otherwise. 

We denote by the Hamiltonian for particles in infinite space. The stability of the 
Hamiltonian is the statement that if at least one of the particle species obeys Fermi statistics, 
then 

Hn > -BN (1.37) 

with B a positive constant independent of N (Dyson & Lenard 1967, Dyson & Lenard 1968). 
The proof has been greatly simplified by Lieb and Thirring with the use of the Thomas- Fermi 
theory. This is reviewed in (Lieb 1976). A new proof (Graf 1997) is based on a remarkable 
electrostatic inequality for the classical Coulomb interaction (1.1). For a review on recent 
developments on these questions, see (Lieb 1990). We mention here that the problem of the 
stability of matter in presence of magnetic fields (Lieb, Loss & Solovej 1995, Fefferman 1995) 
and the quantized radiation field^ (Bugliaro, Frohlich & Graf 1996, Fefferman, Frohlich & 



The matter is non-relativistic and an ultraviolet cutoff is imposed. 
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Graf 1998) has received much attention lately. 
Let 



^- = W1 Y.^'^-Y-Pa, Va = ±1 (1-38) 



TV I 

p. 



be the projection onto the symmetric {rja = +1) or antisymmetric {rja = —1) states of 
particles of species ct.The sum runs on all permutations Pa of the Nq, particles and {rjay" is 
the signature of Pq,; Pa acts both on the position r and the spin variable a of the particle. The 
projection onto the subspace of states of a many-body system with N — Na particles 
having the appropriate statistics is given by 

a p \ a / 

where the sum runs now on all permutations of the N particles that are compositions 
p = pa^ ■ ■ - Pas of permutations of particles of each species. 

Associating to each species a chemical potential Ha, the grand canonical partition function 
of the quantum gas is given by the sum 



{Nc} V V a J J 

-2:2:11 dr,...r. 



^ J: ii^pi^^ ^«pw}I (-/^^A,iv) |{ri, (T„J) (1.40) 

In (1.40), |{rj, cTj}) = J| |rj, cTj) is a product of states of individual particles that also diag- 
onalize the spin component along a fixed direction with quantum numbers ctq taking the 
values -Sa, + 1 , . . . , s^- 

The pressure is defined as in (1.5). Particle and charge distributions are defined as in 
the classical case, replacing the average (1.18) by 

(...) = ^TrP(e-^^^^ • • •) (1.41) 

Here the particle and charge distributions and correlations are defined as in (1.20) and 
(1.22) where now the particle coordinates in (1.19) are the quantum mechanical position 
operators^. 



^We use the same symbol r for the quantum mechanical position operator and for the argument 
of a correlation function. We do not consider spin correlations. All the quantum correlations in 
the review are purely positional with spin variables averaged out. 
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In contrast to the classical situation, the response to an external electrostatic potential 
yields a new type of correlation functions, the imaginary time Green's functions. They are 
obtained by variational differentiation with respect to the external potential. Of special 
interest is the charge distribution at r2 induced by an external potential perturbation at ri 
in the linear response regime 



_ind 



(ri|r2) 



6 



(c(r2)) 



0=0 



(HX'TrPexp(-/3i/A,jv) c(r2)) 



In (1.42) 



- / 0?TCT(ri,T,r2) 



CT(ri,T,r2) = (e^^^'^c(ri)e-^^^'^c(r2)> - (c(ri))(c(r2)) 



(1.42) 



(1.43) 



is the imaginary time displaced charge-charge correlation (or Duhamel function): it reduces 
to the static charge-charge truncated distribution at r = 0. 

For an infinitely extended quantum charged fluid, we introduce the spatial Fourier trans- 
form 5'(k, r) of cr(r, r, 0). The static structure factor of the fluid is then 



S{k) =S{k,T = 0) 



[l.U) 



and the generalization of the classical relations (1.29) and (1.30) to the quantum mechanical 
situation is 



X(k) 



,-1 



(k)-l 



An 



13 



dTS{k, t) 



:i.45) 



Notice that in the latter case, the susceptibility x(k) (defined as in (1.27)) is no longer 
proportional to the static structure factor. 

Equivalence of ensembles and existence of thermodynamic functions was established in 
(Griffiths 1968, Lebowitz & Lieb 1969, Lieb & Lebowitz 1972, Lieb & Narnhofer 1975). The 
analogous results for classical systems are corollaries. 

There are no results on existence of the infinite volume limit of distribution functions 
except for charge symmetric Bose systems with positive-definite interactions.^ In this case 
at any density and temperature for which the system is stable^ all particle distribution 
functions have unique infinite volume limits (Frohlich & Park 1978, Frohlich & Park 1980)^. 



^(1.8) is a positive-definite interaction but hard cores plus Coulomb are not. 

^For sufficiently large activity the grand canonical ensemble diverges. For the ideal gas this 

happens at the Bose-Einstein condensation. 
^Their proof only works for some choices of boundary conditions but this is probably not a physical 

limitation. 
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These results include charge-symmetric classical systems with positive-definite potentials. 
They rely on the Sine-Gordon transformation described later. 

The physical content of these results on the existence of the infinite volume limit is 
that the limit exists independently of the way the container is enlarged. For example, we 
could consider two sequences of increasingly large containers related by a fixed translation. 
Both sequences would have the same infinite volume limit, so the limit must be translation 
invariant. In particular the limiting distributions could not describe a crystal with a low 
density of defects. On the other hand one would expect to be able to select such a state by 
choosing a sequence of containers whose boundaries match the natural planes of the crystal, 
unless the system never forms crystals. Therefore these results are evidence that stable 
charge symmetric systems with positive-definite potentials are always in a fiuid phase. 

Further results on existence of the infinite volume limit of classical distribution functions 
at low density are included in the work on Debye screening (Brydges & Federbush 1978, 
Imbrie 1983a) reviewed in the next chapter, but these are low density results. 



D. Neutrality 

In systems with short range interactions the densities of different species can be varied 
independently, but Coulomb systems maintain charge neutrality 

^Paea = (1.46) 

a 

by expelling excess charge to the boundary (Lieb & Lebowitz 1972, Graf & Schenker 1995a). 
We introduce a vector notation 

e = (ei,...,e5), /x = (/^i, ■ ■ ■ , /^s) (1-47) 
for the set of charges and chemical potentials, and decompose /x 

A* • § 

fx = fx-\- ve, /X • e = 0, v = — (1-48) 

into its components perpendicular and parallel to the charge vector. Then, the charge 
neutrality is equivalent with the fact that the infinite volume limit of grand canonical pressure 
does not depend on the component z/ of /x along the charge vector 

P(/?, /x) = M + i/e) = P(/3, /x) (1.49) 

This means that different choices of chemical potentials /x^ or activities Za do not necessarily 
lead to different densities. It is common to break this redundancy by imposing 

pseudo-neutrality: '^^CaZa — ^- (1.50) 

a 

The system will be neutral regardless of whether pseudo-neutrality is imposed or not, but in 
section II. F we will argue that charges that are expelled to the boundary create a constant 
external electrostatic potential if] which renormalizes the activities Za 2;^ exp[— /Scq,/] in 
such a way as to restore the condition for the renormalized activities. Thus imposing the 
condition amounts to working with the renormalized activities. 
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II. DEBYE SCREENING 



This chapter is devoted to screening in classical Coulomb systems. We begin with a 
review of the original theory of Debye and Hiickel followed by an account of some rigorous 
theorems. Then there follow several sections that illustrate the main ideas used in the proofs 
with discussions concerning open problems. 



A. Debye - Hiickel Theory 

The mean field theory approximation for the dilute plasma phase (Debye & Hiickel 1923) 
is at the center of our discussion. We will briefly review their argument to prepare the way for 
further developments. For further background see (McQuarrie 1976, Hansen & McDonald 
1976). 

The activity Za has dimension length"^ and j3e\ and A have dimensions length. Thus if 
Id is defined by 

iD = (47r^e^^«/3)-i (2.1) 

a 

then //) is a length, called the Debye length. The activities in this definition of Id are 
required to satisfy the pseudo-neutrality condition (1.50). It is standard to define the Debye 
length in terms of densities by the formula 

«-i = (47r5]e^p„/3)-^ (2.2) 

a 

and we have used the notation Id instead of the standard because the two definitions 
are not equal, although their ratio tends to one in the Debye-Hiickel limit described below. 

The conclusions of Debye-Hiickel theory include that the pressure P and charge-charge 
correlation obey 

PiP - Pideal) ^ 



127r/|, 



CT(ri,r2) - - ((47r)^/5Zl,|ri - r^l)"' exp(-^-^) + {Ai^^llyHiv, - r^) (2.3) 

I'D 

The Debye-Hiickel argument suggests that these are really limiting laws, holding in a limit 
in which the number of particles in a volume /|) tends to infinity, while the number in a 
volume A^ tends to zero. Perturbation expansions (Mayer 1950, Meeron 1961, Lebowitz & 
Stell 1968, Hansen & McDonald 1976, Mayer & Mayer 1977) are also hkely to be asymptotic 
as opposed to convergent, but it is still natural to believe that the second result holds in a 
stronger sense: that without taking the limit there can be exponential decay of correlations, 
particularly if one is reasonably close to the limit. We shall refer to this stronger statement 
as Debye screening. 

The first law in (2.3) can be expressed in terms of densities by using po, — Zad(5P/dza 
and /SPideai — Yli^oi- One finds that 
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(2.4) 



The second result (2.3) was based on an argument that can be paraphrased as follows. 
Consider a grand canonical ensemble of particles interacting with each other and also with 
a charge cq held fixed at the origin. This is the same as adding to /?[/ the external potential 
term j3 J eoc(r, uj)V{r)dr. Then the potential at another point r will be the sum of V{r) and 
another potential (I) created by the charges in the gas. Dcbyc and Hiickel approximate this 
(p by assuming that the gas reacts to the combined potential / = CqV + as an ideal gas, 
which leads to the self-consistent scheme 

A/(r, Lj) = -c(r, u}) - eo5{r) Ri -(c(r))ideai,/3/ - eo5{r) (2.5) 

47r 

We wrote /(r.u;) to emphasize that the exact field is configuration dependent but having 
made this approximation we write /(r) from now on. Prom (1.35), 

A/(r) ZaCa exp i-PeJir)] - eoS{r) (2.6) 

a 

If Pf{r) is a weak potential then exp(x) ^ 1 + x and pseudo-neutrality (1.50) leads to 

[A - ln']f{r) = -iireoSiv) (2.7) 

which has the solution 

/(r) = eo|r|-iexp(-|r|//^) (2.8) 

The exponential decay of / is the origin of exponential decay of charge density correlations 
within this mean field approximation. Indeed noting that / is related to the effective po- 
tential introduced in section LB by / = CqV^^ we can obtain the second Debye-Hiickel law 
from (1.32) using the Fourier transform. 

In (2.5) the (c(r)) ideal,/?/ is the mean density of a spherically symmetric opposite charge 
"screening" cloud that forms around the positive charge at the origin. Note that within the 
linear approximation exp {x) — 1 + x our calculations imply 

J dr (c(r))ideai,/3/ ~ -eo (2.9) 

By Newton's theorem, such a cloud is equivalent to an additional negative charge at the 
origin that neutralizes the fixed charge up to a remainder exponentially small in the radius 
to the point r. Of course no single charge configuration of point particles can be spherically 
symmetric. The cloud is an ensemble average over configurations. This cancellation of 
charge is captured in "sum rules" which will be discussed in section 11. B. Related to this 
point is the fact that the integral over r2 of the right hand side of (2.3) vanishes. 
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B. Theorems on Debye screening 

In this section we will survey the rigorous results on Debye screening and the Debye- 
Hiickel approximation concentrating on three dimensional systems but with some remarks 
about other dimensions. 

The one dimensional Coulomb system without hard core has been solved exactly (Lenard 
1961, Edwards & Lenard 1962). Lenard's solution has Debye screening for charge density 
observables, but one dimensional systems have exceptional screening properties, because 
the linear potential —\x — y\ has no multipoles. For example, the electric field is piecewise- 
constant with jumps at the positions of particles so that outside an interval containing any 
neutral configuration, one can achieve a zero electric field. On the other hand if a single 
fractional charge is placed in a system of integer charges then the resulting potential can 
never be screened because the ensemble of integer charges can do no more than cancel 
integer parts of the corresponding electric field. One dimensional Coulomb systems have 
been reviewed in (Choquard, Kunz, Martin & Navct 1981). Since then, highly non-trivial 
results have been obtained on the correlations of classical charges confined to a circle and 
interacting with a logarithmic potential (Forrester 1992, Forrester 1993?), Forrester 1993a). 
Some comments on two dimensional systems are given at the end of this section. 

For three dimensions the results are complicated to state completely and the complica- 
tions are probably artifacts of the proofs. Therefore let us restrict ourselves to a simple case 
of a hard core plasma with two species, whose charges have rational ratio. Recall that 
the combination /5e^ has dimensions of length. By absorbing a fundamental unit of charge 
into 13 we assume that /5 carries the dimension of length while is dimensionless. To put 
it another way, we state results regarding z^, A as parameters, and keeping ta fixed. The 
results below are valid in a thermodynamic limit in which there are boundary conditions 
on the Coulomb potential. Details may be found on page 199 of (Brydges & Federbush 
1980); the boundary conditions are a little artificial but their essence is that the walls of the 
container are grounded, so that all electrostatic potentials vanish at the boundary dK. 

Theorem II. 1 Suppose that 



then all (truncated) n-point particle and charge density correlations decay exponentially on 
length scale I ~ Id- Furthermore, if E\{ri,r 2) is the grand canonical charge symmetric 
ensemble of integer charges with two fractional charges held fixed at ri,r2; then the two 
fractional charges are screened in the sense that 




(2.10) 



a 




(2.11) 




(2.12) 
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The Debye sphere assumptions say that the number of particles in a volume of size 
is small, so that the hard core is out of play, and the number of particles inside a sphere of 
size Id is large. For comparing these statements with the ones in the references it is useful 
to realize that 



so the first Dcbyc sphere condition is equivalent to (3 /Id being small. Theorem II. 1 is 
a combination of results from (Brydges 1978, Brydges & Federbush 1980, Imbrie 19836, 
Brydges & Kennedy 1987). 

It is shown in (Imbrie 1983a) that there is Debye screening in Jellium and in systems with 
irrational charges. Jellium is a limit in which all species with positive charges are smeared 
out into a background charge density by letting their charge tend to zero as their density 
increases. The results of Brydges and Federbush were not proved with enough uniformity 
to persist as this limit is taken. 

The Debye-Hiickel limiting laws should hold in the 



but published proofs (Kennedy 1983, Kennedy 1984) have the much more stringent condi- 
tions that /3/X is kept bounded and the system be charge symmetric. Similar results for 
quantum systems have been established by (Fontaine 1986). The original Debye-Hiickel re- 
sults included correction terms containing a hard corc^ length scale a. More systematically, 
hard core corrections have been obtained by resummations of Mayer expansions (Meeron 
1958, Meeron 1961)^ and by the Kac limit developed in (Lebowitz & Stell 1968) which 
develop corrections in kq. There is also a long history of intuitive extensions of the Debye- 
Hiickel theory to include corrections from hard cores and dipolar effects. These are reviewed 
and advanced in (Fisher & Levin 1993, Fisher 1994, Levin & Fisher 1995, Fisher 1996, Lee 
& Fisher 1996, Lee & Fisher 1997, Zuckerman, Fisher & Lee 1997)^. However consistency of 
these approximations is guesswork: is there an asymptotic expansion in Ka? What physical 
effects provide the next largest corrections to Debye-Hiickel theory? 

In the proofs of screening complications at the boundary have been avoided by artificial 
boundary conditions, but only laziness stands in the way of proofs for the case where the 
particles are in a conducting grounded container so that the Coulomb potential vanishes 
at the boundary. Particles with 1/r forces, i.e., particles in an insulating container, have 
not been completely analyzed, but Debye screening for a simplified model is established 



^ P{P - Y.a Pa) ~ -^cr(Ka) where a{Ka) = {^1 + kg - - 2 ln(l + ko) j . It may be 

better to consider a to be a length scale that characterizes "effective" short range forces created 

by cooperation between Coulomb forces and hard cores length as opposed to a genuine hard core. 

'This is compared with extended Debye-Hiickel theories in (Bekiranov & Fisher 1998). 

^There is a competing program (Stell 1995, Yeh, Zhou &: Stell 1996) based in part on the mean 
spherical approximation instead of Debye-Hiickel theory. 




a 



(2.13) 



Debye-Hiickel Limit: ZaX — > 0; Zalo^ 2a — > oo 



(2.14) 
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(Federbush & Kennedy 1985). This case is discussed some more in section II. G. It is 
difficult because screening fails at the boundary; there are power law forces between charges 
on the boundary. 

It should be possible to omit the pseudo-neutrahty condition (2.10), but to do it rig- 
orously may require insight into the boundary condition problem. In section II. F we will 
discuss these issues along with charge expulsion. 

We return to (2.9). The ~ involves approximations but the theorems known as sum 
rules say that if truncated distributions decay integrably or better then this is exact in the 
limit Co — > 0: the screening cloud around a fixed infinitesimal charge neutralizes it. More 
generally any fixed set of charges surrounds itself with a screening cloud in such a way that 
all multipoles are canceled. In particular, the function of r defined by 

CT(r|ri,e„i,...,r„,ec,„) = 
(c(r)p(ri, e^J • • • p(r„, ea„)) - (c(r))(p(ri, e^J • • • p(r„, e^^)) (2.15) 

which represents the excess charge density at r when particles in the system are fixed at 
ri, . . . , r„ has no multipole moments in the Debye regime 

J drgi{r)cT{r\ri,ea„...,rn,eaJ =0 (2.16) 

for all harmonic polynomials gi{r) of order / = 0, 1, . . .. In (2.16) the contribution of coin- 
cident points r = Fj is included. The case / = 0, go{r) = 1, called the charge sum rule, is 
expected to hold very generally in homogeneous phases of Coulomb systems (classical and 
quantum). These results are reviewed in (Martin 1988). Results on the effect of slow decay 
at boundaries are also covered in this review. In the next section, we present a derivation 
of the sum rules with the help of the Sine-Gordon formalism, extending the arguments of 
(Fontaine & Martin 1984) to charges with hard cores. 

There are also results (Alastuey & Martin 1985) that say that if truncated distributions 
are integrable and decay monotonically at infinity then they must decay faster than any 
inverse power. A different result of the same genre was obtained by (Federbush 1979). 

Two dimensions, no short range forces: in units where the Coulomb interaction between 
two unit charges is — Inr let j3i = 2,f3oo = 4 and define the sequence of intermediate 
thresholds /3„ = /3oo(2n — l)/(2n) with n = 1,2, . . .. The Coulomb potential is not stable 
in the sense of (1.11) but the instability is weak enough to permit thermodynamic behavior 
anyway (Hauge & Hemmer 1971, Deutsch & Lavaud 1974, Frohhch 1976), at least for P < Pi. 
At Pi the Gibbs factor exp —P1V2 for two oppositely charge particles is no longer integrable 
and the partition function diverges. The partition function for the Yukawa gas also diverges 
and, in the Yukawa gas Mayer expansion, the diagram with two vertices becomes infinite 
aX P = Pi = 2. However, all diagrams with more vertices remain finite. Similarly, at 
/?2, neutral diagrams with 2, 4 vertices diverge, all diagrams with more vertices remain 
finite, and so on. It has been proved (Brydges & Kennedy 1987) that for P < 4/3/3i the 
Mayer expansion without the infinite two vertex term is convergent for small activity. The 
natural extension of this result to the higher thresholds almost certainly holds for P < P^ ^. 



^An almost equivalent result appeared in (Dimock & Hurd 1993), but the authors have reported 
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These results imply that there are natural infinite renormalizations of the partition function. 
Presumably, this enables some version of the Yukawa and Coulomb gases to be defined in the 
range P G [Pi, Poo) as a renormalized limit of systems with regularized interaction, because 
dropping a divergent term in the Mayer expansion is the same as multiplying the partition 
function by an infinite factor. This has not been discussed clearly in the literature, but see 
(Gallavotti & Nicolo 1985 Speer 1986). 

It has been proven (Yang 1987) that the two species charge symmetric Coulomb gas of 
point particles in a grounded container has exponential screening if /3 <^ /3i and the activity 
are small. The grand canonical ensemble with insulating boundary conditions is harder to 
analyze because it has to be restricted to neutral configurations, whereas image charges 
provide neutrality for free in the grounded container case. 

As the last paragraph suggests, boundary conditions are more important in two dimen- 
sions than in three. For example it is known (Frohlich & Spencer 19816, Theorem 4.1) that 
a pair of oppositely charged fractional charges immersed in a two dimensional Coulomb gas 
are not screened from each other when the thermodynamic limit is taken with insulating 
container boundary conditions. On the other hand, although the details are not published 
anywhere, it ought to follow from the results of Yang that they are screened when the ther- 
modynamic limit is taken with grounded containers. This infiuence of boundary conditions 
does not happen in three dimensions. We will discuss this further in section II. G. 

Two dimensions, with stabilizing short range forces: The significance of Poo is that the 
Kosterlitz-Thouless transition takes place at a critical Pc ~ Poo^^- By stabilizing the inter- 
action with short range forces such as hard cores one can reach Pc without infinite renor- 
malizations and see the P > Pc dipole phase. The original argument (Kosterlitz & Thouless 
1973) was given a complete proof (Frohlich & Spencer 1981a) for the Coulomb system on a 
lattice. A more detailed analysis of the Kosterlitz-Thouless phase, based on term by term 
analysis of the Mayer expansion, is given in (Alastuey & Cornu 1992, Alastuey & Cornu 
1997c, Alastuey & Cornu 1997a). There is confusion over whether the thresholds Pn should 
be interpreted as the successive collapse into neutral clusters of 2, 4,. . . particles. (Gallavotti 
& Nicolo 19856, Speer 1986) proposed such an interpretation based on successive divergences 
in the Mayer expansion. This is analogous to the Yukawa gas alluded to above, except that 
now infra-red divergences cause the thresholds. To complete their argument one must show 
that there are bulk observables whose correlations have associated singularities. On the 
other hand (Fisher, Li & Levin 1995) argue that no singularities show up in thermodynamic 
functions and consequently the thresholds are not physical. 

The heuristic renormahzation group argument (Kosterlitz & Thouless 1973) predicts that 
exponential screening of particle distributions holds for all P < P^ii the density is sufficiently 
small. A complete proof of screening at low activity for P in this range would be a wonderful 
achievement. 

Exact formulas in two dimensions for the excess free energy of a single charge fixed at the 
origin are conjectured in (Lukyanov & Zamolodchikov 1997). See also (Smirnov 1992) for 
other interesting formulas. These formulas in principle could be used to express the excess 



a serious error invalidating their proof (Dimock 1998) 
^°/3c — Poo as the density tends to zero 
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free energy of a fractional charge in terms of the density ( or the activity) and temperature. 
The authors study the Sine-Gordon functional integral with boundary conditions that corre- 
spond to a grounded conducting container so that violating neutrality by fixing a fractional 
charge makes sense. Their conjecture translates, via the Sine-Gordon transformation to the 
two component charge symmetric gas with a purely Goulomb potential at least for (5 < 
With limitations on the size of the fractional charge their formulas apply for /3 G [/3i,/3oo)- 

In two dimensions, there are a number of exactly solvable models for the special value of 
the temperature that corresponds to starting with the two dimensional one component 
plasma (Jancovici 1981). These models (Alastuey 1987, Jancovici 1990, Jancovici 1992) 
exhibit "super-Debye screening", namely Gaussian decay of the correlations. They have 
played an important role in testing various refined screening properties in presence of walls 
and of different types of inhomogeneities. 

Although a characteristic of Debye phases is the rapid decay of the particle correlations, 
one should be aware that it is never the case for potential and electric field fiuctuations 
(Lebowitz & Martin 1984, Martin 1988). For instance, correlations of the components of 
the electric field E^j^ always have an asymptotic dipolar behavior 

{E,{v)EM) ^ -^^mV. (2.17) 

This observation has recently led to the development of the interesting view point that there 
are some universal properties of Debye phases that closely resemble those of a system with 
short range forces at its critical temperature (Jancovici, Manificat & Pisani 1994, Jancovici 
1996, Forrester, Jancovici & Tellez 1996). Actually, potential and field fiuctuations are 
critical in Debye phases, and this leads to universal finite size corrections to the free energy 
as in critical systems. 

C. The Sine-Gordon transformation 

1. The Fourier transform of a Gaussian is a Gaussian 

In chapter I we have defined a potential 

Ko,A(r) = |r|-i(l-exp(-|r|/A)) 

which is a Coulomb potential smoothed at the origin. We saw that Voo,x has a positive 
Fourier transform. For such potentials there is an identity that expresses the interaction 
in terms of external potentials — an exact version of the mean field idea of Debye- Hiickel. 
This is the Sine-Gordon transformation (Kac 1959, Siegert 1960). Kac and Siegert noticed 
that if self-energies are included, then the Gibbs factor 

g-l 1:5=1 e«,e„.y,o,x(ri,r,) ^ ^-f / dr / dr' c(r)y^.;,(r,r')c(r') 

is Gaussian in c(r). By a functional integral generalization (Simon 1979) of "the Fourier 
transform of a Gaussian is a Gaussian" there exists a Gaussian probability measure, intu- 
itively described by 
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(2.18) 



dpLx{(l>) = P(/)e-^^'^("-)^->("-)'^'' (2.19) 

such that 

g-f E5=l-"i-«i^oc,A(ri,r,) ^ y (^^^(0)e-^^'/'/*£W'^{r). (2.20) 

V^^;^ is the operator inverse of the operator whose kernel is Voo,\- In our case we find from 
(1.10) that it is the partial differential operator 

V-\ = — (-A + A^A^) (2.21) 

possibly with boundary conditions on A . Equation (2.20) says that the Gibbs factor is a 
superposition of Gibbs factors for external potentials It follows that if an ideal gas of 
particles in an external potential is integrated over the external potentials then the result 
is a partition function for a gas with two-body interaction Voo,x, including self-energies. 
Self-energies are equivalent to a shift in the activities. We define 

z^^\£) = ^(£:)et^«^-'^('"''") (2.22) 

Then the partition function with potential energy U (cu) built out of Voo,a and an external 
potential ip, 

Sa(/3^/V)= / da;z(a;)e-'^^(-)-^'''-^^^(-)<^'- (2.23) 

J A 

is given by 

Sa(/?V2^) = J dl^M 5ideal,A,.(>^) + ^])- (2-24) 

where 

Sideai,A,.(^)(/) = / «?c.zW(a;)e-/^^* (2.25) 

is the ideal gas partition function (1.34) with the renormalized activities (2.22). (2.24) is the 
Sine-Gordon representation. It says that the partition function for the interacting system 
exactly equals a (functional) integral over imaginary external fields of the ideal gas partition 
function. 

Gaussian measures have been studied in the mathematical literature and it is known 
that for this one the typical potential is a continuous function (Glimm & Jaffe 1987, 
Theorem A. 4. 4). This result holds provided A 7^ 0; when A = the instability of the 
Coulomb potential returns to haunt us in the shape of the typical (p becoming a distribution 
instead of a function. 

The right hand side of formula (2.19) is best understood as a mnemonic for the translation 
formula, also known as the Gameron-Martin formula. 
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(2.26) 



which is vahd for any function g for which J gV^\g exists^^. Translations 4> ~^ ~ig are 
legitimate as well whenever the integrand permits analytic continuation to a — i after the 
change of variable — > + ag. The translation formula is the reason why computations 
with the intuitive formula (2.19) are usually correct. 

It is not physically very reasonable to try to capture the very local fluctuations on length 
scales much smaller than /^i in a mean field picture. For example, the activity z^^^ is very 
different from z when A is small, but we shall see in section II. E that z is the correct effective 
activity, in the sense that an ideal gas with this activity is the best approximation in the 
Debye-Hiickcl limit. Also hard cores do not have positive Fourier transforms and cannot be 
represented by a Sine-Gordon potential. For these reasons we now consider a better class 
of representations in which the particle picture is retained on scales less than 0{Id) and 
the mean field or Sine-Gordon is used to represent interactions on larger scales. Define the 
Yukawa interaction with a decay length L > X, 

VL{r) = Irl-^e-l"-'/^; VlM = |r|-^ (e-l""'/^ - e-l""'/^) (2.27) 

We split the interaction T4o,a according to 

Voo,x = Ko,L + Vl,x (2.28) 

and use the Sine-Gordon transformation to represent just the part Voo,l- The result is 

where 

^LAzw{f)= I dujz^''\uj)e-l^^'^^^^-S^^^'^^'^'' (2.30) 

is no longer ideal; Ul is built out of two-body interactions a- If hard cores are used in 
place of the cutoff A in Vl,\ then [/l is built out of Vl and hard cores. 

The replacement of V by Voo,l is the same as smoothing the Coulomb interaction by 
introducing form factors for the charges so they are smeared out into spherical distributions 
of characteristic radius L. Note that 

T4o,L(k) = F{L'k')^^F{L^k% F{L'e) = (^q:^) (2-31) 

which shows that F is the Fourier transform of the form factor. 



In the sense /(| ^ -|- | XAg p) dr < oo. 
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2. Sum rules 



Sum rules provide a good illustration of aspects of the Sine-Gordon transformation that 
will reappear in the sequel. Here we sketch how the argument in (Fontaine & Martin 1984)^^ 
applies to a Coulomb system with hard cores. Wc will show that sum rules result from an 
invariance of the Sine-Gordon measure under infinitesimal translation by harmonic functions. 
Formulas arising from infinitesimal symmetries are called Ward identities in quantum field 
theory. While it is interesting that sum rules are Ward identities, it is also easy to derive 
them without the Sine-Gordon transformation (Martin 1988, Brydges & Martin 1999). 

We consider charged particles with hard core interactions of radius A in a box A whose 
boundary is grounded and impenetrable to the hard cores. We spread each charge in the 
ensemble uniformly onto the surface of its hard core sphere. By Newton's theorem the 
interaction energy U is unchanged by the spreading out. However the self-energy of each 
particle becomes finite and of the order of It can be added into U and the change 

in U is compensated by replacing Za by z^^ — Zq, exp(/3e^/(2A)). By the Sine-Gordon 
transformation 

di_i{(l)) = e-^/(^") /a(^^)'i?0 (2.32) 

where 4>\{r) is the average of over a sphere of radius A and center r and 5a(0) is the hard 
core gas in external field 0: 

-a(0) = II "at I / dri... j dvN cxp - ^ eaj(i)irj) - hard core (2.33) 

{N^}a=l \ j ) 

dfj,{(l)) is characterized by 

dlj{(l) + g)= ^^(0)el/(47r)/^0Ag-l/(8,r)/^(as)2. ^2.34) 

where g is any function which vanishes on the boundary with Jj^{dgy < oo. 
The complete expectation 

-^{^^^^=1 J A J A 

X exp {-pU - hard core) F (2.35) 
can be taken in two stages; the first stage is an expectation conditioned on 



{Nc}a 

X exp 



eaj(f){rj) — hard core^ F 



This paper uses the smooth short range regularization (1.8). 
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and the second is the remaining integration over 

E[F] = J e'-^^ACi^i/V) F I J e^--A(i/3i/^0) ^^^^^ ^^.sg) 

and these fit together by 

(F) = E[{F\ijj)] where -0 = iP^^'^4>x (2.37) 

To obtain a sum rule we make the change of variables cf) ^ (f) + tg in the integral 
on the right hand side and differentiate with respect to t at t = 0. The derivative must 
vanish since the left hand side is independent of t. The derivative of lnSA(0) with respect 
to 0(r) is — (c(r)|0), which accounts for the first line in the next equation. The second 
line arises because the derivative of {F\(p) with respect to 0(r) is the truncated expectation 
-(Fc(r)|0) + (F|0)(c(r)|0). 

+ E[{Fc{gx)\ip) - {F\xp) {c{gx)\^)] + contributions from (2.34) (2.38) 
The contributions from (2.34) are proportional to 

EiiFliP) J <f>Ag] - E[{F\'iP)]E[ J <f>Ag] (2.39) 

We would hke to choose to be a harmonic polynomial because then these contributions 
vanish. However this is not a legal choice of g because no harmonic polynomial can vanish 
on the boundary of A. Instead choose g to he a harmonic polynomial multiplied by another 
function h that is one in the interior of A and which goes to zero near the boundary. 
Now there arc non-vanishing contributions from J (pgAh and J (pdgdh localized near the 
boundary. Suppose that the container A is enlarged: there is competition between the growth 
of the harmonic polynomial g{r) for r near the boundary and the decay of the truncated 
expectation (2.39). If the system is in a screening phase where truncated expectations of 
local functions of (p decay exponentially then the polynomial growth of g is crushed as the 
container is enlarged and the contributions in (2.39) vanish in the thermodynamic limit. It 
is true that the conditional expectations (c(r)|^) are not exactly local in but if the Mayer 
expansion of (c(r)|'0) is convergent then (c(r)|0) is almost (exponentially) local in 0. This 
captures a physically natural condition that the system should not be so dense that there 
are long range correlations caused by the hard core. If the Mayer expansion converges at 
■ip = then it converges uniformly for all ip because ip 7^ just means that activities are 
multiplied by modulus one phase factors. 

Thus in the infinite volume hmit (2.38) simphfies to = E[{Fc{gx)\ip)] — 
E[{F\iIj)] E[{c{g\)\ip)] where 5^ is a harmonic polynomial. Now we simplify further using 
(2.37 ) and Newton's theorem (spherical average of harmonic = value at center) gx{r) — g{r) 
to reach 

= {Fc{g)) — {F){c{g)) g any harmonic polynomial (2.40) 

For example by choosing F = p(ri,ai) . . .p(r„,a„) (2.40) becomes the same as (2.16). In 
conclusion, the sum rules (2.16) for a system of classical charges with hard cores hold if the 
Mayer expansion for the hard core system without Coulomb interaction is convergent and if 
the system is in a screening phase. 
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D. Debye spheres 



In the previous section we divided the interaction into the Yukawa part Vl^x whose decay 
length is L and a remaining part, the slow part V^,l, which carries the slowly varying long 
range piece of the Coulomb interaction. In this section, and in detail in section VIII. B, the 
Yukawa part is shown to leave the gas essentially ideal, provided its range L is of order the 
Debye length Id and the Debye sphere assumptions (2.11) are imposed. The underlying idea 
is that the effect of the V^^x part of the interaction is to associate some of the monopoles 
into aggregates which can be regarded as new species so that the resulting multi-species gas 
is ideal. The Mayer expansion is used to express this. Instead of leaving out the long range 
part of the interaction, one can instead say that the gas inside a sphere of radius less than 
the Debye length is close to ideal. This illuminates the Debye-Hiickel argument in which 
the charge density was computed by an ideal gas formula: the mean field captures the part 
of the interaction with range larger than L and the interactions of shorter range are small 
enough to use an ideal gas calculation. 

Consider the partition function for a Yukawa gas in an external potential (p, 



The Mayer expansion is an expansion for the logarithm of the partition function or the 
pressure P in powers of the activities Za whose first term (3PV — '^Za + --- gives the ideal 
gas law when = 0: it has the form 



where are given (Ruelle 1969, Hansen & McDonald 1976) as sums of contributions from 
connected "Mayer" graphs. In particular ui{£) = 1 so that the m — 1 contribution to the 
product is 



which equals the ideal gas partition function Sidcai,A(/3^^^0) that appeared above in (1.34). 

The higher m terms in (2.42) have the form of additional ideal gas species, e.g., the 
m = 2 term is an ideal gas partition function for a gas whose particles are the possible 
aggregates of two of original particles, with a cluster activity ^z{£i)z{£2)u2{£i, £2) ■ Such a 
cluster describes two particles bound together because the functions Um decay exponentially 
in the separation of any pair of arguments, because the graphs that contribute to them are 
connected. Although these "Mayer" aggregates do capture effects of aggregates of particles 
they are not directly physical because their activities need not be positive. 

This decomposition of the interacting gas into species of ideal gases is only useful if the 
activities of large aggregates are small. In fact we need them to be summable which is the 
same as demanding that the Mayer expansion is convergent. The standard condition (Ruelle 




(2.41) 



Sl,a(/9^/^0) = 

exp ( J] /" n d£,z{E,)e-^''''">''t'^'^''^Urn{£u ■ ■ ■ , <^m) 



\m>l k 



(2.42) 
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1969) for convergence for the activity expansion for a gas with activity z{S) and two-body 
potential v{S,S') is 

sup j \e-^<''^'^ -l\\z{8)e-^""'f'^^^\dS'e^''^ <e-^ 

where B is the stabihty constant, i.e., the best constant such that 

v{Si,Sj)>-BN. 

l<i<j<N 

Suppose for simplicity that there are two species with equal but opposite charges Cq = ±1 
and equal activities 

Using the stability estimate (1-11) we find the Mayer expansion converges if 

|z|/3L2exp(/3/A) < 1 (2.43) 

which we calculated using exp{±VL^x) — 1 ^ ±Vl_x. This condition is sufficient for conver- 
gence of the Mayer expansion, but it is not a very accurate condition when A < /3 < L. In 
section VIII. B we will show that there is a better condition, namely 

zph^ < 1; z(3L'^e^ < 1. (2.44) 

The P/P emphasizes why the estimate is an improvement. For the Yukawa potential and 
hard core, one needs 

ZaX^ < 1 (2.45) 

as well. Now observe that when L — ^Id the second condition in (2.44) drops out, provided 
7 is small, but independent of z, A. The other two conditions are implied by our Debye 
sphere hypotheses (2.11), when combined with the definition of Id- Thus the Debye sphere 
hypotheses say that the gas inside a sphere of radius L = 7/^) is a multi-species ideal gas 
whose activities are small according to 

Y[d8k\z{Sk)\\u^{Sl,...,Sn^)\<0{ — \ \z\\A\ (2.46) 

k \ ^ / 

where z is the largest of the Za- For details on the derivation of this estimate see section 
VIII.B. 

The conditions for convergence of the Mayer expansion when there is an external potential 
are the same except that the exponential of the potential is absorbed into the activity. Thus 
an imaginary potential has no effect in (2.46) because of the absolute value in |2;(£fc)|. 
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E. The Debye-Hiickel limit 

The standard derivation (Mayer 1950) of the thermodynamics, the first part of the Debye- 
Hiickel law (2.3), proceeded by resumming a class of ring diagrams in the Mayer expansion 
or low density expansion. A more systematic basis for this resummation is based on the 
Kac limit, see (Lebowitz & Stell 1968). In this section, following (Kennedy 1982, Kennedy 
1983, Kennedy 1984), we sketch a proof of the Debye-Hiickel law based on the fact that 
the representation (2.29) becomes an explicitly integrable Gaussian integral in the Debye- 
Hiickel hmit (2.14). The Debye-Hiickel mean field emerges as the stationary point for the 
Sine-Gordon action. This line of argument has the advantage of not using any expansion 
which might be divergent. 

Our definition of the Debye-Hiickel limit is low density in the extreme sense that the 
hard core becomes irrelevant and the effects of association into dipoles and other aggregates 
are out of play. The original Debye-Hiickel theory (footnote 6) has terms in Ka that would 
disappear in this limit if a is taken to be the hard core radius. Debye-Hiickel theories (see 
the discussion below (2.14)) have been advocated as a good description for some transitions 
in ionic fluids at higher densities in which phase separation occurs. These are extensions of 
Debye-Hiickel theory in which hard core corrections remain and effective dielectric constants, 
which represent the effect of charges pairing to form dipoles, are used. For consistency con- 
siderations it is desirable to show that these theories are also limiting cases of the Coulomb 
system or obtainable in some systematic approximation. 

Let us work in units where Id = 1- This means that in the right hand side of the 
representation (2.29) we rescale lengths by r = l^r' so that the Debye length becomes 
V D — 1 and r' becomes dimensionless. Then we have exactly the same formulas but the 
parameters are replaced by dimensionless parameters: 

(3' = (3 /Id, X' = X/Id, z'a = zJl, L' = L/Id, A' = Id'A (2.47) 

In terms of the primed parameters the Debye-Hiickel limit becomes 

z'a^oo, /3'^Owith { (2.48) 

( z'J'le-P'/'^' ^ oo 

From now on in this section all formulas involve the primed quantities, so we will save on 
notation by dropping the primes. 
We write (2.29) in the form 

Sa(/3V2^) = J dfXLicP) exp (F(i0 + ^)) (2.49) 

where 

F(i0 + V') = In S^,A,,(.) (i/3 + /3V2^) (2.50) 
We define the action functional S{ip,g) by 
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-Siij, 9) = \j 9V-\gdv + F{g + ^) (2.51) 

The partition function Sa(/3^/^V) does not depend on L in the left hand side of (2.49) so 
we may take the infinite volume and Debye-Hiickel limits at fixed L and then afterwards let 
L — s> 0. The correct order of limits is infinite volume limit followed by Dcbyc-Hiickcl limit, 
but we do it in the opposite order. A serious technical point in the papers we cited is to 
show that these limits can be interchanged. 

In the Debye-Hiickel hmit the action S{'4>, i(f)) becomes a quadratic polynomial in (f). We 
forestall the appearance of some terms linear in (p by using the translation formula (2.26) 
with (j) ^ (f) — ig applied to (2.49) 



Sa(/3V2^) = e^hy-],9dr j ^^^^^^ 



X exp 



j i(t)V-\gdY + F{i(t) + ^ + g) 



(2.52) 



g will be chosen later. By (2.46) the Mayer expansion for F(i0 -\- g -\- is convergent 
uniformly in the Debye-Hiickel limit, provided L <^ 1: indeed the 0{z^^^) term in the 
Mayer expansion is 0{L^)z\K\ uniformly in the Dcbyc-Hiickcl limit. The 2;|A| is a common 
volume factor. Divide it out and fix L <C 1. Then the first term in the Mayer expansion 
dominates the sum of all the others uniformly in the Debye-Hiickel limit. Therefore the 
Debye-Hiickel limit can be taken term by term under the sum in the Mayer expansion and 
under the integral by standard theorems of analysis (dominated convergence). Every 
term becomes a quadratic polynomial in and only the first term in the Mayer expansion 
contributes in the limit L — > 0. To see this in more detail, note that the first term is 



/ drz\^^e-'-^"'"^'^'-'+^^ (2.53) 



By (2.48) terms of order 0{z^a'^ P^^'^) vanish in the Debye-Hiickel limit. The 0{z^^ P'^^'^) term 
vanishes in the hmit by pseudo-neutrality (2.10). Therefore, in the Debye-Hiickel hmit, the 
non- vanishing part of F{i(f) + g + is 

FBK{i(f> + g + ip)^ j^dr + + g + + i-v;^_^(o) j 

+i?2(# + ^ + V') + ^o (2.54) 

where 14o,l(0) arose from the difference between Za and zi^^; see (2.22). R2 is a quadratic 
remainder term 

R2{f) = J y"/(ri)r(ri,r2)/(r2)ciricir2 (2.55) 

Rq is independent of f = i4> + g + ip. Rq and R2 account for all the remaining terms 
in the Mayer expansion. There is no term linear in /: inspection of 0{f) contributions 
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from the Mayer expansion is required to verify that the pseudo-neutraUty condition makes 
them all vanish in the Debye-Hiickcl limit. Furthermore Rq and R2 are independent of z. 13 
and vanishing as L ^ 0. If we substitute Fdh back into (2.52) and regroup into terms 
independent of 0, terms linear in and terms quadratic in 0, the result is 

SA/Sideai,A ~ e-^^^^^'^^ j d/iM e^B^ -^Al-^'-K^.i^W)* g-5'DH(V,9;i0)+fl2W) (2.56) 

Note that the activity in the ideal gas Hidcai,A is z, whereas a naive calculation based on 
(2.24) might draw us into using z^'^\ The ~ means that we keep only the terms that are 
non-vanishing in the Debye-Hiickel limit. We have defined the quadratic functional 

-S^ui^, 9)^\j 9V-\9 dr+j-J{g + i^f dr + R^(^j; + g) + R^ (2.57) 

whose first variation in direction is 

-5'dh(^, g;<P) = J (i>V-\g dr + j-J<P{g + ^) dr + i?^(^ + g; 0) (2.58) 

R2{ip + 9] 4>) is the first variation of i?2(V' + 9) direction 0. It is natural to choose g to 
make S{^^ vanish, but instead we settle for making it almost vanish by choosing g to solve 

-Ag + g + 'ip = Q (2.59) 

because then by (2.21) 

-S'Mip, g; ct>) = ^L" j <pA^g dr + R'^iiP + g; 0) (2.60) 

which means that S'^uiip, 9',4>) disappears at the end when we take L — > 0. However we 
will keep these terms around for a few more lines because we want to derive some formulas 
that will be used in later sections. 
Define 

5dh,a = J dfiL{<P) e-3^/A(<^^-^^,i^{o))*+«2(i^) (2.61) 
which is the normalization for a new Gaussian probability measure 

c^/xdh(0) = 55h,a^/^40) e-8^^A(*'-^-'-(°))'^^+^^(^^) (2.62) 
Then (2.56) is the same as 

SA/Sideai,A ~ e-^°«(^'^) Sdh,a j df,M<f>) 6-^'°^^,.;*^) (2.63) 

As we have explained above, when L — 0, we can drop the 5"dh so that the normalized 
(i/^DH also disappears and we are down to the evaluation of Sdh,a- 
The normalization Sdh,a can be expressed as a determinant using 
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/_ 1 
d/.L(</.)e-^-''A'^'* = det(/ + ^V;o,L) ' (2.64) 

We have dropped the R2 term because (2.46) imphes that 

\R2{i(p)\ < 0{L^) [ (P^dr (2.65) 

J A 

and so it gives no contribution in the hmit L — > 0. To obtain this estimate from (2.46) 
note that R2 involves two variational derivatives with respect to 0. Each derivative brings a 
factor These factors combine with z in (2.46) and become a Debye length Id — 1- We 
want a = 1, but we put it there because taking a = shows that both sides are correctly 
normalized. \4o,l is the operator whose kernel is the potential Voo^ii^i — ^2). The arguments 
ri,r2 are confined to the finite volume A either by boundary conditions on the Coulomb 
potential, if the particles are in a grounded conducting container, or by restriction on the 
range of integration, if the particles are in an insulating container. It can be shown that 
the thermodynamic limit washes out the difference, but the easier case is the first; then 
there are zero boundary conditions on the Laplacian in \4o,l so that the eigenfunctions are 
trigonometric functions when A is a box. Using det(/ + X) = exp[^^ln(l + (^)], where ^ runs 
over the eigenvalues of the operator X, one can show that the leading term in the infinite 
volume limit is 



=^DH,A ~ exp 



-l(27r)-3|A| I dk (lnil + l-V^,,)-l-V^,,^ 



(2.66) 



where the Fourier transform V^^l is given in (1.13). The limit as L — > of this expression 



IS 



exp 



i(27r)-3|A| I dk(ln(l + !)-!) 



= exp[|A|/(127r)] (2.67) 



Returning to (2.63) we have shown that that the leading terms in the Debye-Hiickel limit 
followed by infinite volume limit followed by L — > are 

In (SA/Sideai,A) ^^ + ^J 9i-^)9dr + ^ J (g + i^fdv (2.68) 

where g satisfies the stationarity condition (2.59). When = we recover the thermody- 
namic part of the Debyc-Hiickel law, (2.3). Furthermore the stationarity condition coincides 
with the linearized Debye-Hiickel equation (2.7) in units where = 1 when 

/ = p-y\g + V), V'(r) = eo/3^/V(r) (2.69) 

The strange factors of appear because of the way we normalized the Sine-Gordon trans- 
formation. The Debye-Hiickel result (2.3) on correlations, in units where Id = 1, is obtained 

by taking two functional derivatives of the approximation (2.68) for In (HA/Sjdcai.A) with re- 
spect to —P^^'^ip. This interchange of functional derivatives with Debye-Hiickel limit might 
be hard to justify, but convexity helps in the case of the single charge density expectation 
(Kennedy 1983). 
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There may be inequahties comparing the Debye-Hiickel limit with the exact system. 
Some results follow from Jensen's inequality, but they are limited to the unrealistic case of 
charge-symmetric systems (Kennedy 1982). See also theorem III. 2 in the next section. 

Questions: Are extended Debye-Hiickel theories such as the ones considered by Fisher et 
al. hmiting theories? One should no longer let L — > and should choose g to be the exact 
stationary configuration for S'dh(i--5')- which contains quadratic remainder terms R2 giving 
corrections to the dielectric constant. Also the activity may be renormalized, whereas in 
our case L ^ causes z^^'^ to coincide with z. Is the expansion (Lebowitz & Stell 1968) 
asymptotic in ko? 

F. Charge expulsion 

We continue to use units in which l^, = 1 introduced in section II. E. In our preliminary 
discussion of the Debye-Hiickel theory in section II. A we claimed that the activities should 
be constrained by the pseudo-neutrality condition (2.10). If we start with a system that 
does not satisfy this condition, then even when ■0 = 0, the action S{Q, icp) defined in (2.51) 
is not stationary in at = 0. However, a translation (f) ^ (j) — ig such that S{0, g + i(f)) 
is stationary in (f) at (f) = reveals the physical effect of charge expulsion and an attendant 
renormalization of activities. This is a nonlinear effect, not small enough to fit within 
the approximation of the last section. The physical interpretation is that excess charge is 
expelled to the boundary where it forms a boundary layer of total charge proportional to 
the surface area. This layer modifies the electrostatic potential in the interior in such a way 
that the activities are renormalized to satisfy pseudo-neutrality, which is why the excess 
charge is of order surface area and not volume in this grand canonical ensemble. 

For simplicity we omit the complications that arose in the last section due to the length 
scale A by considering a system with interaction Voo,l with insulating boundary conditions 
and no other forces. In this case we can use the Sine-Gordon representation (2.24) to 
represent the entire interaction. We can also apply the translation formula (2.26) for — > 
(f) — ig. The result is 



exists (Kennedy 1984) a unique g that minimizes this functional. Uniqueness is an easy 
consequence of it being convex. Let t/j — 0. Then in the deep interior of A is close to the 
constant go that minimizes 






Insulating boundary conditions mear 
tends over infinite volume but the S: 
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^^^^,..^MP"'9o) = [ d£ z^^\£)e-^'''^"^^ (2.72) 
Outside A, g solves V^\g — and — > at infinity. Define renormalized activities 2;^^'^^ 

by 

zf'^) = exp(-e«/3^/2go)4^) (2.73) 

and note that 

Sideal,.(^)(/9'/V) = Sideal,.(^.«)(/3'/'[/ ' 9o]) (2.74) 

Since g ^ go inside A the bulk contribution to S 00,1(9) comes from the ideal gas term which 
has the renormalized activities z^^'^\ But since Ag is non-zero near the boundary there 
is a boundary contribution to the pressure in | / gV^\gdr and also in the ideal gas term 
because g ^ go si the boundary. These arise from the surface charge expelled from the 
interior and account for the capacity term in the free energy, discussed in the analysis of 
systems with net charge (Lieb & Lebowitz 1972, page 55), (Graf & Schenker 1995a) 
From (2.70) 

S.(0) = e-^--(^) / d,MeI^^<^^^^ S,,^,,..,.)(/3V^[i0 + ,- go]) 

J =ideaM(i.«)(/3^/^[^-yo]) 

By our choice of translation g we have arranged that the icj) fluctuates around zero so that 
exp {—Soo,l{9)) should be the leading term in the Debye-Hiickel hmit. The next term in 
the approximation would be a Dcbyc-Hiickcl term that comes from integrating over the 
fluctuations using a quadratic approximation as in section II. E. 



G. Symmetries and tunneling 

We continue to work in units in which Id — 1- There is a loss of intuition entailed by the 
Sine-Gordon transformation because the integration extends over imaginary potentials, but 
there are also aspects one can see more easily in the Sine-Gordon language. In particular if 
all charg integral multiples of a fundamental unit of charge, which for simplicity we 

assume to be one, then the action S{il>, g) in (2.51) is invariant under the discrete symmetry 

S{ilj,g)^S{ij,g + i5h), 5h^2T:(5-^'^ (2.75) 

This section is concerned with the breaking of this symmetry, which plays a role in the 
proofs of screening. We will see also that there are tunneling corrections associated with the 
symmetry. 



^^If L = and g is chosen to be constant in the interior of A and harmonic in the exterior, 
then the boundary contribution in (2.71) exactly equals the electrostatic energy stored in a perfect 
conductor with shape A. 
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Let us first note that the Coulomb system in an insulating container in three or more 
dimensions cannot enjoy this symmetry because the Gaussian measure d^L^cj)) satisfies 

j diiLi<i>) 0'(O) = Ko,l(0) < oo (2.76) 

precluding any symmetry of the form (j) ^ + c. 

We have seen in section II. E that the Debye-Hiickel approximation is a quadratic ap- 
proximation around the stationary configuration = 0. The symmetry (2.75) tells us that 
there are infinitely many stationary configurations related by the symmetry to = 0. We 
claim that even though (2.76) implies that the stationary configuration = is favored, in 
the thermodynamic limit there are arbitrarily large regions where (f) is trapped in the other 
stationary configurations. This is analogous to the Ising model at low temperature in which 
there occur arbitrarily large islands of — spins even though the boundary is set at + and 
the majority of spins are +. In the zero temperature limit in the Ising model, the islands of 
— spins disappear, but if one wants to understand low temperature as opposed to zero tem- 
perature, then they are important. In the same way, to obtain screening near, as opposed 
to in the Debye Hiickel limit, the presence of large islands where the = approximation 
breaks down has to be taken into account. 

To take into account the presence of many stationary configurations related by the sym- 
metry we use a functional version of the Villain approximation 

gCos(3;) ^ g-i(a;-27rn)2+l (2.77) 
n 

on the term exp(F) in (2.49). In (8.80) in section VIII. E we explain in terms of this one 
dimensional analogue why taking into account the extra wells in this way improves the 
convergence of perturbation theory. 

exp (F) is invariant under a translation (p ^ (f) + h where h — h{r) is any function of 
r that takes values in the set {n6h} of periods, because 5^^^ ^.(l) clearly has this property. 
Also from (2.54) it has the quadratic approximation 

FMip + ^<^) = ^ + ^(^'^ + + 

+R2{i(j) + ij) + Ro (2.78) 



Thus an analogue of (2.77) is 



i4>) ^FY)Yi{i>-ih+i(j>) (2.79) 



h 

where is an 0(/9^/^) error term representing cubic corrections to each well and put there 
to obtain exact equality, h = h{r) is summed over all functions that vanish outside A, take 
values in the set {nSh} of periods and are piecewise constant on a lattice of unit cubes filling 
A. 

If we think of a typical piecewise constant function h{r) as a height, then the resulting 
landscape splits into plateaus where h{r) is constant but more important are the jumps from 
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one plateau to another. In analogy with the Ising model we call the connected surfaces in A 
where h jumps "contours." Our objective is to argue that near the Debye-Hiickel limit the 
sum over h is dominated by terms in which these contours are small and very dilute. The 
value of /i in a given plateau fixes the well in which is trapped and the contour is where 
(f) tunnels from one well to another. 

We substitute (2.79) into the Sine-Gordon transformation (2.49) and retrace the analysis 
of section II. E leading up to the Debye-Hiickel approximation (2.63) with ip replaced by ip—ih 
and g replaced hy ig 

h 

X J dnY,n{(t>) e-^'°^^'^-'''^'3;i'l>)+R3{f+i9+i4>) (2.80) 

For each h in the sum g{r) can be any function. We choose it to make 5'dh('0 ~ ih,ig) 
stationary so that S"dh vanishes. 

For the moment we drop the term R^, then the (i/^DH integral drops out because it is 
normalized and we are left with an approximation 

5A/5ideai,A ~ Hdh.a ^ e'^'^^^^-^'^'^^) (2.81) 

h 

where 

SDniip - ih, W) = ^J gV'lLgdr +^ Ji^-h- iipfdr 

+R2 + Ro (2.82) 

When ip — g minimizes SDii{—ih,ig) by trying to make {g — hY vanish but cannot 
quite succeed because / gV^\gdv forces g to be smooth whereas h has jumps. Suppose 
h{v) = /ir(i") has only one contour F and h vanishes at the boundary of A. Then F has 
an exterior "sea-level" plateau where h = and an interior plateau where h ^ 0. In the 
interiors of both plateaus g{r) pa h{r). The error in the ~ is 0(exp (— dist(r, F)). 

Any h with h — at the boundary can be decomposed in a unique way 

h^Y.^^ (2.83) 

j 

where each hj has only one contour Fj and hj vanishes at the boundary. Thus hj is defined 
on the whole box A but it assumes just two values: zero outside the contour Tj and /ij(int) 
at all points in the interior of the contour. Across the contour Tj hj jumps by hj{mt) which 
is the same as the jump of h across F^. 

The minimizers are linear functions of h so g — ^ gj where gj is the minimizer for the 
single contour hj. If the contours are well separated then by g{r) — /i(r) -|-0(exp (dist(r, F)). 

'S'dh(V^ - ih, ig) ^ ^ 5'dh(^ - ihj, igj) (2.84) 
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Thus 

SA/Sideai,A ~ Sdh,a J] J] e-'^^^^-""^''^^^ (2.85) 

h j 

The approximation in these two equations was to leave out exponentially small cross terms 
in S'dh involving {gi — hi){gj — hj) and derivatives of Qi times derivatives of Qj. 

The sum over h is equivalent to summing over the number N of contours, the shapes Tj 
of the contours and the jumps hj{mt) at each contour so that 

N ■ ri,...,rjv j ftj(mt) 

The exponential of S'dh('0 — ihj,igj) suppresses the jump hj{mt) at Tj by a factor 

exp[-0(|r,|)/i,.(int)2] (2.87) 

where IF^I is the area of the contour. The ^j(int)^ dependence is because S'dh is quadratic 
in g, h when ip = 0, gj is linear in hj and hj is just hj{mt)x a unit jump configuration. 

In (2.86) the term with no contours, = 0, is our original Dcbyc-Hiickcl approxima- 
tion (2.68). From (2.86,2.87) we see that the contours should behave like an ideal gas with 
an exponentially tiny density near the Debye-Hiickel limit. Nevertheless they represent a 
phenomenon that probably makes the standard perturbation series asymptotic but not con- 
vergent and incapable of proving that there is screening near as opposed to at the limit. 
In contrast (2.86) is the basis for a convergent expansion which does imply screening near 
the Debye-Hiickel limit. The terms we left out in arriving at this representation are expo- 
nentially small interactions between the dilute contours. There are also constraints (hard 
core interactions) in the sum that forbid intersections of contours. Such interactions are 
within the purview of conventional convergent Mayer expansions, which is why it is possible 
to develop a convergent expansion. However this expansion is clumsy in comparison with 
standard perturbation theory. One has to show for example that the error term does not 
destroy the picture we have just developed. One of the principal accomplishments (Glimm 
& Jaffe 1987) of the constructive quantum field theory program was to develop a calculus 
for the analysis of "almost Gaussian" integrals and this is what is used at this point. The 
reader who suffers from residual curiosity can turn to chapter Vlll for a compressed and 
updated account of this calculus and its applications, or to (Brydges & Federbush 1981) for 
an older review. 

Note the role of grounded container boundary conditions is to anchor g and therefore 
h at zero on the boundary. The 1/r potential with no boundary conditions, i.e. particles 
in an insulating container, is more difficult than the case of grounded boundary because 
the anchoring is weaker, since it comes from (2.76). This has been partially investigated 
(Federbush & Kennedy 1985). Their argument is hard because there are long range forces 
at the boundary which make it difficult to control the size of boundary effects. 

In the case of two dimensions with insulating boundary conditions, there is no anchoring 
at the boundary because the Gaussian measure^^ d/iL does not break the symmetry (2.75). 



defined by taking a limit outside the integral as a small mass tends to zero 
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We mentioned in section II. B a result of Prohlich and Spencer (Frohlich & Spencer 1981 &, 
Theorem 4.1) that fractional charges are not screened when insulating boundary conditions 
are imposed. This is because of this symmetry. Intuitively there is no preference for which 
potential well (j){ri) is trapped in, but once it has made up its mind, 0(r2) wants to choose 
the same well, which is a long range correlation. This argument predicts that the correlation 
decays to a constant. On the other hand, if the boundary is grounded, then the well is 
selected essentially everywhere, except for a very dilute gas of contours that enclose regions 
where is in other wells. In this case all local observables are expected to have exponential 
decay, but this is only proved for integer charge observables (Yang 1987). Recall in analogy 
that the Ising model at low temperature with plus boundary conditions has exponential 
decay but this decay is lost if the boundary spins are not tied down because the state is no 
longer pure. It is an open problem to prove that two dimensional systems with insulating 
boundary conditions have screening of integer charge observables. 

Returning to three dimensions when there are irrational charges the symmetry is re- 
placed by a quasi-periodicity, nevertheless there is screening in these systems without exact 
symmetry (Imbric 1983a). The argument requires Pirogov-Sinai theory. 

There is an intriguing difference between fractional charge observables and integer 
charge observables. Suppose the gas consists of unit positive and negative charges but 
there is an additional half-integer charge fixed at the origin r = 0, represented by a fac- 
tor exp (— /3^/^i0(O)/2). The 0(0) fluctuates around g{0) ^ h{0), so there is a factor 
exp (— /3^/^i/i(0)/2) which is ±1 in a way that depends on all contours surrounding the 
origin, which is a more non-local (topological) effect than when an integer charge is placed 
at the origin. To put it another way: a phase factor exp (— /3^/^i/ij(0)/2) must be included 
in the sum over hj{mt) in (2.86) whenever Tj encloses the origin. This does not destroy 
screening in three dimensions, but it might be worth further study. 

We will show, in appendix VIII. E, that the tunneling corrections are of size 
exp {—0{j3 /Id)) so they are not of any consequence in the Debye-Hiickel asymptotic regime. 
However, they should be important if l3/ljj is not particularly small because then the period 
is small and the potential wells that trap (j) are close together and shallow. There could be 
a phase transition in which the imaginary potential icj) ceases to be well localized in wells. If 
as a consequence there are more fluctuations in the typical 0, then the gas will be organized 
into more tightly bound neutral aggregates because when the position r of a cluster or parti- 
cle with charge Cq, is integrated, the rapidly varying phase factor exp(— ■ieQ/9^/^0(r)) will give 
cancellations. Note that the Kosterlitz-Thouless phase transition fits this description. 

Tunneling represents a phenomenon which arises because the system consists of discrete 
charges as opposed to infinitely divisible charge distributions, which would give a single well 
as in the Debye-Hiickel approximation. 



We thank James Glimm for conversations on this point. 
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III. DIPOLES 



The Coulomb plasma is the main focus of this review. However in chapters IV and V 
we will find that a system of point quantum charges has a close analogy to a certain system 
of dipoles and higher order multipoles. The resulting multipole forces are at the origin of a 
breakdown of exponential screening in quantum mechanics. It is therefore of interest to give 
here a short review concerning the lack of screening in classical dipole systems and estimates 
on the dielectric constant. 

Concerning the analogy an important caveat should be kept in mind. In a classical system 
exponential screening is restored as soon as free charges are added to classical dipoles. At 
any non-zero temperature, the quantum gas always has a proportion of free charges due 
to ionization processes, but these free quantum charges are not able to restore exponential 
screening. The analogies and differences between the behavior of quantum point charges 
and classical dipoles are made more precise in section IV.F 



A. The Dipole ensemble 

A dipole is specified by the coordinate S = (r, d) that unites the position r and the 
dipole moment d. We suppose that in the absence of interaction a single dipole experiences 
no preferred directions so that 

Jd£... = J dr JdQ...with J dQ ^ 1, (3.1) 

where dQ is surface measure on the sphere of unit vectors d and d = cid with d a fixed 
dipole moment. 

The interaction between a pair of dipoles is 

V{S,, 82) = (di • Vi) (d2 • V2) V{v^ - v^) (3.2) 

where Vj is the gradient operator that acts on rj. V{yi — Y2) should be the Coulomb potential 
between two unit charges at ri,r2, but the singularity of the Coulomb potential at short 
distance will lead to instability for the dipole system. We choose a length scale L and a 
form factor F and set 

\/(k) = (47r)|k|-2F2(LV) (3.3) 

The form factor satisfies J F(r)dr = 1 and means that point charges are replaced by "charge 
clouds." This should be regarded as an effective interaction which arises by integrating out 
forces on length scales less than L in a more realistic description of the short range physics. 
An effective action of this form is not appropriate unless the gas is dilute so that interactions 
such as hard cores are not playing a strong role. See for example the use of a Mayer expansion 
as in section II. D which led to such an effective interaction. 
The partition function is 



J2'^JdS^...Jd£^e-^''^'--'-\ (3.4) 



N 
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where the interaction includes self-energies, 



1 ^ 

U{£r,...,£N)^-^V{Si,£j). (3.5) 



If instead the interaction is ^ f{£j) where / is an external field, then we obtain an ideal 
gas partition function 



"ideal.A 



W = Em /<Ki...y'<<£7,expUf;/(f,)j 
= exp(./..e-.) 



As in section II. C we can reconstruct the interaction between the dipoles by integrating 
over an imaginary external field, f{£) = ■ Vi0(r), with respect to the measure d/iL 

discussed in (2.19): 

Sa = y" dflM Hideal,A(/3'/M • Vi0). (3.6) 

Formally this is J T>(f)exp (— 5'(i0)) where we define the action S by 

-S{9) gV-'g + zjd£ exp(-/3^/M • V^) 

^^'^^ + z j d£ cosh(/3^/M ■ V^) (3.7) 

where the inverse operator has a kernel whose Fourier transform is y~^(k). 

We shall consider the observable that measures the density of dipoles with coordinate £ 
in a configuration uj — (fi, . . . , £n)-i namely 



N 



p{£,uj)^Y.^{£,£^), (3.8) 



Distribution functions such as the two point function {c{£a)c(£b)) are defined in analogy 
to section LB. Distributions at non-coincident points of c{£,uj), where £ — (r, d), become 
expectations of 2;exp(— i/?^/^d • V0(r)) in the Sine-Gordon language . 
Another useful quantity is 

(e-Z^^^^Iiff) ^ s^i J Si,,,,^(i/3VM . V0(r))e-^^^^/^/ (3.9) 

which measures (/?x)energy of a charge distribution / in the sea of dipoles. By unraveling 
the Sine-Gordon transformation one finds that the charge distribution / has self-energy 
determined by the potential V^(ri — r2) corresponding to (3.3) and it interacts with system 
dipoles according to potential energy (d2 • V2) V{ri — r2). 
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We are using the expected value ( ) to indicate either the finite volume or infinite vol- 
ume limit (s). This infinite volume limit is known to exist for several choices of boundary 
conditions (Frohlich & Park 1980, Frohlich & Spencer 19816, Funaki & Spohn 1997). 

We shall also consider expectations of (j). The two-point function measures the effective 
potential between two infinitesimal charges: as in (1.31), but in the Sine-Gordon language, 
/3x the effective potential between two infinitesimal test charges (r^, ea)a=i,2 is given by 

/ -iei/3V2^(ri) -je2/3i/2,/,(r2)\ 
p-/3eie2V<=ff(ri-r2) ^ 6 )_ 

(e-»ei^'''^</'(ri))(e-»e2/9i/2^(r2)) ^ ' 

The numerator removes the self-energies of the test charges. Since ei, 62 — >^ we have 

yefr(r^_r2) = (0(rO0(r2)) (3.11) 
Furthermore the dielectric constant e is obtained from the large r asymptotics, 

V^^ir) ~ ^ (3.12) 
£ r 



B. A no-screening theorem for dipoles 

The action S{g) in (3.7) is invariant under g ^ g + c, where c is any constant. This is a 
continuous symmetry and the action is local. Locality and continuous symmetry generally 
lead to power law correlations because fiuctuations, — > where /(r) is slowly varying 

as a function of r, are very weakly suppressed . 

No-screening theorems for dipoles have been obtained (Park 1979). Subsequently 
(Frohlich & Spencer 19816) gave the proof we present. This no-screening theorem is a 
little remarkable because it is valid for all activities and all temperatures, but its validity at 
high density is dependent on the way the short distance singularity in the dipole force is sta- 
bilized. If there are hard cores, then crystalline phases are likely. This is proven (Frohlich & 
Spencer 19816) for dipoles on a lattice. For models where the Coulomb potential is smoothed 
out such as our Voo,\ potential Frohlich and Spencer proved further results to the eff'ect that 
the two point function has no oscillations, suggesting that this system is always in a fiuid 
phase. 

Theorem III.l For all activities z > and all inverse temperatures (3 > there is no 
screening in the sense that the two point function has a singularity at \s. — in its Fourier 
transform and 

V^'^ik) > V(k), (3.13) 

where 



^^The Coulomb system possessed a similar symmetry, but with the essential difference that c 
could only have discrete values. 
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V(k)= {v-\k) + ^k^) (3.14) 



is the Fourier transform of the dipole potential with a dielectric correction. 

V(k) is the Fourier transform of the kernel V(r — r') which is the inverse of the Hessian of 

5* at = 0. A version of this theorem holds in any dimension and also for a large class 
of positive-definite interactions in place of dipole-dipole. The factor 3 is the dimension of 
space. We will remark in the proof the main property of the interaction being used. The 
reader will see that the argument is a very general mean field theory bound. 

Prom this theorem and discussion of the dielectric constant in section LB we immediately 
obtain a mean field theory bound for the dielectric constant, assuming there is a dielectric 
constant: 

£ < 1 + ^47r (3.15) 

There is a large literature on the molecular theory of the dielectric constant, with atoms 
and molecules modeled by preformed classical permanent dipoles, using various short range 
regularizations. See, for example, (Stell, Patey & H0ye 1981) and references therein. In 
particular, the Clausius-Mosotti law {e — l){e + 2)^^ = 47rp/3(i^/9 was derived in (H0ye & 
Stell 1974, H0ye & Stell 1976) by taking the Kac limit, in which V{y) is replaced by ^i^V{^i'y) 
and 7 — > oo. However, this law is contradicted by the bound (3.15) if one just sets z = p 
at low density. The reason is that the result of the Kac limit depends on the treatment of 
dipoles in the vicinity of the test charge and the scaling of test charges. The literature would 
benefit from clarification of this point. See the further comment at the end of this section. 

Theorem IH.l is a consequence of 

Theorem III. 2 (The mean field upper hound) 



where 



(e-/^<A/) <exp(^inf-5(^,/)j (3.16) 

-S{g, f) = \j gV-'g -Jgf + zJdS (cosh(/3VM . v^) - l) (3.17) 

Proof. We make a complex translation (f) ^ (j) — ig. This is done using the transla- 
tion formula (2.26), but we arrive at the same place if we do it in the intuitive formula 
J V(f) exp{—S{i(j)) — J i(f)f), which is the numerator of {e^^'^'^^). On expanding out the 
exponent we obtain an upper bound on |exp(. ..)| by dropping all the imaginary terms. 



Thus 



-Re5(# + (7)-Re J{t<p + g)f 
= Re^ I + g)V-\i(j) + g) - j gf + J cosh(/3^/M • V(i(/) + g)) 

^~IJ W^+l j gV-'g- jgf + ^j rf^:cos(/3VM-V0)cosh(/3VM.V^) 

<-S{t<j))-S{g,f) (3.18) 
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so that in the resulting bound on (exp(— J i(f)f)) the J V(f)cxp (—S{i(f))) cancels in numerator 
and denominator and we obtain the mean field upper bound, since g is arbitrary. QED 

Proof of Theorem III.l. We rewrite the mean field upper bound as 

1- (e-'^'^'f'f) >l-exp{-S{ag,af)) (3.19) 
which holds for any function g. Let a — > to obtain 

la J 4>m >-\j gV-^g + j gf-\zjd£ ■ V^)' (3.20) 

This is valid provided g has enough decay so that the terms of higher order than are 
convergent integrals. This and positive-definiteness are the only properties of the dipole 

interaction used in the proof of (3.13). 
Now choose g = Vf. The result is 

(( / 0/)') > / fVf. (3.21) 

This is the same as 



fV'^f > I fVf (3.22) 



so that 



V^^{k)\f{k)\'dk> J V(k)|/(k)|^dk (3.23) 

which, being true for a large class of /, is the result we claimed. QED 

A Kac limit can be defined by replacing V{r) by 'y~^V{'yr) and letting 7 — >^ 0, which is 
the same as taking L — > 00. By scaling one can show that it is equivalent to replacing 5'(i0) 
by 7~^5'(i0) in the Sinc-Gordon transformation (3.6,3.7). This selects the stationary point 
of the action so that bounds such as Theorem III.l are saturated in this limit, provided the 
external charge density / is also scaled suitably. 



C. The scaling limit of the lattice dipole gas 

The results of the last section show absence of screening but one would like to know that 

the system is in a phase characterized by correlations that are asymptotic to canonical power 
laws prefaced by a dielectric constant or tensor. In this section we will describe preliminary 
attempts to delineate when the system is actually in a dielectric phase in this more detailed 
sense. 

Unfortunately the most appealing result is established only for lattice systems, but the 
removal of the restriction to lattice models is a feasible mathematical problem. 

In a lattice model each site r in a simple cubic lattice in a z/ > 1 dimensional box A can 
be occupied by zero, one or more dipoles, each of which can only point in a lattice direction. 
The Sine-Gordon transformation of such systems leads to a partition function of the form 
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V(j)e-^^^^ (3.24) 



where Vcf) — nre(2)dnA ^^(^) ^ finite dimensional integration and the action 5" has the 
general form 

5(0)= W{<p{v)-4>{v')). (3.25) 

<rr'> 

where < rr' >=< r'r > is a nearest neighbor pair of sites. For example, a dipole gas with 
no interactions other than the dipole interaction with lattice regularization is given by 

W{ct>{v) - ct>{v')) = i-(0(r) - cf>{v')f - cos {d(3''\ct>{v') - 0(r))) . (3.26) 

The first term is the finite difference version of (V0)^ and the second term corresponds to 
the cosh(i/?^/^d ■ V0) we have seen before. If there is a hard core interaction preventing 
more than one dipole per lattice site then {z/u) cos is replaced by ln(l + {z/u) cos). Some 
of the results we are about to state were derived assuming periodic boundary conditions on 
the lattice. 

This partition function is also called the anharmonic bedspring, in which case the vari- 
ables 4> are considered to be displacements of the nodes of a bedspring (in a distinguished 
direction) . 

General results (Naddaf & Spencer 1997) have been obtained under the assumption that 
S is convex and is a local (or nearly local) function of the gradient of 0, namely assume 
there is a positive number S > independent of such that the matrix 

of second derivatives obeys 

5] C(r)5"(0; r, rOC(r') > S (CW - C{r')f (3.28) 

<rr'> <rr'> 

for all C- The finite difference Laplacian A associated to the lattice is a matrix defined by 
E (C(r) - C(rO)' = E C(r)(-A)(r,r')C(rO. (3.29) 

<rr'> <rr'> 

Thus this convexity assumption is, by definition, the matrix inequality 

Convexity : ,5"(0) > 5(-A) (3.30) 

This type of assumption has a natural continuum analogue: a lower bound on the sec- 
ond variation of yS" by 5 J (V0)^ or perhaps a quadratic form in —A -|- L^A^, but the next 
assumption brings in the lattice in a strong way. 

The assumption is that there is a constant C such that 
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J2 C(r)5"(0; r, r')C(rO < C CW^- (3.31) 



<rr'> <rr'> 



This is not as natural as the convexity assumption and is false for a continuum model whose 
action contains S J (V0)^ or higher order derivatives. 

Now we come to the scaling limit which concerns expectations of functions of fields 
averaged over a length scale i that is taken to infinity. For any smooth compactly supported 
function / with ^ /(r) = define 

</.(/,)= ^0(r)£-/^+V(^)- (3.32) 

r 

The X] /(i") — means that (j){fe) is really an integral of V0 against a function of compact 
support. 

The theorem will say that when £ is large 0(/^) becomes Gaussian. Gaussian random 
variables are characterized by their covariance, which, in the theorem, will be the continuum 
Green's function C(r,r') solving 

- Yl eijdidjC{r,r')^5{r-r'), (3.33) 

i,j=l,...,d 

where di = d/dri and £ij will be the dielectric tensor. 

For the theorem we need to define the Gaussian measure 

diic{(p) = ^V(pexp ^-i j Y^eijdi(p{r)dj(p{r)dr^ (3.34) 

Such a Gaussian measure exists (essentially the massless free field ) and is characterized by 
Wick's theorem which says 

Theorem III. 3 Define the functional Laplacian 

Let P be a polynomial in fields 4>{ri), . . . , 0(r„), then 

J diic{4>) P{4>) = exp Qa^ P(0)U=o (3.36) 

where exp (|A) is a power series I + |A + • • • which terminates after a finite number of 
terms when applied to a polynomial P. 

An immediate consequence is J (iyUc(0)0(r)(/)(r') = C(r, r') which is why C is called a co- 
variance. The theorem (Naddaf & Spencer 1997) is 
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Theorem III. 4 Under the convexity assumption (3.30) and the uniform upper hound (3.31) 
the continuum scaling limit of the anho,rm,onic bedspring is a continuum massless free field. 
In other words there exists a constant positive-definite dielectric tensor Sij such that, for any 

(11 '^(/^))^ / dl^cWW) (3.37) 

i i 



Thus, in dimension z/ = 3, the lattice dipole gas (3.26) is in a dipolar phase characterized 
by a dielectric constant at least when ^z^d^ii: < 1 because this is the parameter range for 
which (3.26) is convex. 

In mathematical terms this is a new type of central limit theorem. Central limit theorems 
are usually a statement about sums of independent random variables. In the statistical 
mechanics of lattice spin systems in the high temperature phase the spins are somewhat 
independent and one expects a central limit theorem that the sum of the spins in a large 
block of side i normalized by is almost Gaussian and independent of similar sums over 
disjoint blocks. A beautiful theorem of this type has been proved (Newman 1980, Newman 
& Wright 1981, Newman & Wright 1982). The most gentle way to encode high temperature 
( approximate independence) is to require only a finite susceptibility which remarkably is the 
assumption made by Newman, together with a type of ferromagnetism. In the case of the 
anharmonic bedspring the susceptibility is not finite; correlations decay with non-integrable 
power laws and the variables (f) and their gradients are very far from independent. Now we 
no longer have a standard central limit theorem; indeed the right normalization for a block 
is now but nevertheless these block sums become Gaussian. 

From the point of view of the renormalization group the central limit theorem is the case 
of a high temperature fixed point obtained by scaling keeping the mass fixed so that the 
renormalization group dynamics drives the system to a white noise Gaussian fixed point. 
Theorem III.4 instead describes a massless Gaussian fixed point. Notice that the convexity 
hypothesis permits interactions that are not small. 

The proof sharpens the estimates of (Brascamp & Lieb 1975, Brascamp & Lieb 1976), 
in particular the Brascamp-Lieb bounds: 

Theorem III. 5 Let F — F{(t)) he any continuously differentiable function of fields. Under 
the convexity assumption (3.30), 

<^VW^<E(^S'-(,;r,rO^) (3.38) 

r,r 

and 

(e^(^)) < exp |^iX;/(r)(-5A)-^(r,r')/(r')j (3.39) 
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It is a corollary of theorem III. 4 that all reasonable functions of fields converge to 

massless free field expectations. Unfortunately the standard observables are not included in 
this class, being local functions of the unaveraged field 4>{r). For example the dipole density 
is 

p{S; 0) = ^e-^'^'^'^^-^'^W (3.40) 
Theorem III. 4 is more or less equivalent to 

i) as a ^ oo (3-41) 

but to get at the density observables we need information on asymptotics when some but 
not all of the points ri, . . . ,r„ are driven apart by a scaling. Theorem III. 5 provides an 
upper bound but lower bounds are not yet proven. There is no physical reason to doubt 
that they hold, but we are interested in a proof because if the Naddaf-Spencer argument 
could be extended to prove that more general functional of (f) have power law decay then a 
complete proof that there is no exponential screening in quantum Coulomb systems is within 
reach using the strategy in (Brydges & Keller 1994a). An alternative route was started by 
(Brydges & Keller 19946) but it is presently unreasonably complicated. 

One further remark about dipole systems is that the pressure and dipole distributions 
are analytic in the activity in a small neighborhood of zero (Gawedzki & Kupiainen 1983, 
Brydges & Yau 1990). Thus the Mayer expansion is convergent at small activities. Both 
proofs are complicated and despite many efforts, no one has yet obtained the result by direct 
attack on the size of the Mayer coefficients. 



III-9 



IV. SEMI-CLASSICAL COULOMB GAS 
A. The Feynman-Kac representation 

In this chapter, we give an introdTiction to some of the main effects produced by quantum 
mechanics in the Coulomb gas. These effects are most easily seen by perturbing around the 
classical gas, i.e. neglecting the quantum statistics and keeping the first relevant contribu- 
tions to an expansion in powers of the Planck constant h (the Wigner-Kirkwood expansion) 
We have no control on the possible convergence or asymptotic character of these ^-series, 
but the examination of the lowest order terms gives immediately a prediction to the main 
issue: Debye screening cannot survive if the quantum mechanical nature of the charges is 
taken into account. This serves as a gentle training for the reading of chapter V where 
the fully quantum mechanical gas is considered, since most of the arguments and of the 
mathematical structure will be similar. 

The formalism that is best adapted for both a semi-classical and a low density analysis 
of the quantum Coulomb gas is the Feynman-Kac functional integral representation. This 
representation is of course not new and goes back to (Ginibre 1965, Ginibre 1971) with his 
study of the low activity expansion of quantum fiuids with short range forces. But, for quan- 
tum charges, the formalism shows its full capability and becomes particularly operational 
to deal with the long range of the Coulomb potential. 

In section IV. C we discuss the simpler situation of only two quantum charges immersed 
in a classical plasma. The model can be studied without recourse to h expansions and 
illustrates clearly how quantum fiuctuations will irremediably destroy Debye screening: one 
finds an decay of the correlation between the two quantum charges. In section IV. D we 
comment on the relation between this type of correlation between individual charges and 
the conventional van der Waals potential. Section IV.E comes back to a system of infinitely 
many quantum charges and summarizes the results obtained by the semi-classical analysis. 
In the last section, we combine the Sine-Gordon and the Feynman-Kac functional integrals: 
this gives an alternative view on the lack of exponential screening. Again a simplified model 
that captures the essence of quantum fluctuations illustrates how the occurrence of dipole 
(and higher multipolc) fields will spoil the proof of Dcbyc screening. In contrast, mixtures of 
classical particles that carry both charges and dipoles still show strong screening properties. 

We start here and in section IV. B by recalling famihar ideas about the representation 
of the quantum mechanical statistical operator by the Feynman-Kac functional integral. 
For the sake of simplicity, we consider first a single particle of mass m in three dimensions 
moving in an external potential V^(r). According to the original path integration (Feynman 
& Hibbs 1965), the configurational matrix elements of the statistical operator associated 
with the one-particle Hamiltonian 

H = A + V, A = Laplacian in (4.1) 

2m 



^^Note that to have stability of the classical reference system, the semi-classical analysis requires a 
regularization of the Coulomb potential at the origin, except for the one component plasma which 
remains well behaved in the classical limit. 
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read 



(4.2) 



In (4.2) S'(r(-)) is the classical action corresponding to the potential —V (the Euclidean 
action) 



I3h 



dt\- 



dr{t) 



dt 



+ Virit)) 



(4.3) 



associated with the path r{t) starting from ri at "time" t — and ending in r2 at "time" 
(3h. The summation in (4.2) runs over all such paths. 

It is very useful to parameterize the path r{t) by dimensionless variables, making the 
change 



s — 



I3h ' 



< s < 1 



where 



ri,2(s) = (1 - s)ri + sr2 + \i{s) , «(0) = «(1) = 



x = nJl 

V m 



(4.4) 



(4.5) 



is the de Broghe thermal wave length. In (4.4) ^(s) represents now a closed path, called the 
Brownian bridge, starting and returning to the origin within the "time" s = 1. In terms of 
these dimensionless variables, the matrix element (4.2) can be written in the precise form 
(Simon 1979, Roepstorff 1994) 



(ri|exp(-/?(-^A + V 



r2 



27rA 



3/2 

2 , exp 



|ri - r2| 
2A2 



(4.6) 



In (4.6) -D(^) is a Gaussian measure (the Brownian bridge measure) with weight formally 

ds 



represented (up to normalization) by exp | ~^ Jq 
defined by its covariance 



It has zero mean and is entirely 



J D{$,)^^{si)^uis2) = (5/,,^(min(si, S2) - S1S2) 



(4.7) 



where ^^{s) are the Cartesian coordinates of ^(s). 

The representation (4.6) has several advantages: physical parameters ri, r2, /3, h, X 
appear explicitly in the formula. Moreover the diagonal matrix elements 
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have the form of an integral over a certain ^-dependent Boltzmann factor exp(— /3F(r, ^)) 
with l/(r,0 = /;dsy(r + A^(s)), as in classical statistical mechanics. In this view, one 
can think of a quantum point particle as a classical extended object, the closed filament 
^(s) located at r. This closed filament, which plays formally the same role as a classical 
internal degree of freedom, represents the intrinsic quantum fluctuation with size controlled 
by the de Broghe wave length A. The statistical weight of a filament having extension R 
behaves as exp(— ii!^/A^); thus, in the semi-classical regime (A small), only short filaments 
contribute to the integral (4.8). Obviously one recovers the genuine classical Boltzmann 
factor exp{—f3V{r)) in (4.8) as A — > 0. For these reasons, the Feynman-Kac representation 
is particularly suitable to study the semi-classical aspects of equilibrium quantities. 



B. The gas of charged filaments 

The extension of the representation (4.8) to the many particle system is straightforward. 
In the sequel, V{ri — r2) represents a two-body regularized Coulomb potential, for instance 
of the form (1.8) . To each particle labeled by i, i = 1, . . . , TV, we associate an independent 
Brownian bridge ^j(s), distributed with the Gaussian measure -D(^j), and the diagonal 
elements of the N- particles statistical operator read 

N N 

({r.} I exp(-/?i7^,^) I {r.}) = JI (2^A|.)^/^ J ^l^^^^) 



exp I -- 2^ ea^ea^ dsV{ri + XaM^) - rj - Xa,ij{s)) ) (4.9) 

The Dirichlet boundary conditions on i?7v,A are implemented by the constraint that all paths 
in the integrand of (4.9) have to remain inside A for all s, < s < 1. We do not write 
this constraint explicitly since it will be trivially removed in the infinite volume limit. The 
grand-canonical Maxwell-Boltzmann partition function is defined by keeping in (1.40) only 
the term corresponding to the trivial permutation. Denoting it still by 5a in this chapter, 
it is given by the sum 

^ jj((2a,, + l)exp(;i//„))^^" 



{N^} 

X 



'a- 



[ dri... drN{{ri} \ exp {-(3Ha,n) \ {r^}) (4.10) 

J A 

The factor 2a a -|- 1 is the spin degeneracy, since the Hamiltonian (1.36) is independent of 
the spin. 
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A look at the Boltzmann-like factor in (4.9) suggests naturally the introduction of an 
enlarged classical phase space whose elements T are 



^=(«, r, (4-11) 

where a is the species and ^ is the filament shape of a particle at position r. We shall call 
a point T in this phase space a "filament". Since the Hamiltonian is invariant under the 
permutations of the particles of the same species, we can use again the summation identity 
(1.14) to write the partition function (4.10) in the form 

= E / n ^•^fe^(-^fc) exp (-/3C/(^i, . . . , J'n)) (4.12) 
Ar=o ■ fe=i 

where the phase space integration means 

s 



f dT---^Y. f D(^)--- (4.13) 

a=l 



In (4.12) we have set 



Z{T) = (2(7a + l)Za, Zo, = ,^ (4.14) 

(27rA^) 



N 

U{Ti, ...,Tn)^Y. ^oci^c.,V{Ti,Tj) (4.15) 
\=i<j 

with 

V{J^i,T2)= [ dsV{r^ + X^,^,{s)-r2-Xa,Us)) (4.16) 

It is clear that the system defined by the relations (4. 12)- (4. 16) has the structure of a classical 
statistical mechanical system of objects having internal degrees of freedom a, ^ (the species 
and the shape of a filament), with activity (4.14) and two-body interaction (4.16). We call it 
"the system of filaments" and can apply to it all the standard methods of classical statistical 
mechanics. 

In particular, one can define correlations between filaments in the usual way (see (1.19) 
and (1.20)). Introducing the "5-function" which identifies the points J-'i and ^2 in phase 
space, i.e. 

J dT2?>{T^,T2)F{T2) = F(J^i), K^x^T^) = ^a„o.Mx.i2)K^i - r2) (4.17) 

the density of filaments and the two filament distribution are 

p{T) = {p{T)), p{Ta.T,) = (p(^a)p(^.)) (4.18) 
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where 



p(^) = ^5(^„^) (4.19) 

i 

and (■ ■ ■) denotes the grand-canonical average in the system of filaments (the contribution 
of coincident filaments, included in (4.18), will be specified later if needed). However, one 
should remember that the shapes of the filaments are not directly observable quantities: the 
physical particle distributions are obtained by integrating them out in (4.18) and (4.19). 
Thus the density of particles of type a is 

p(a,r) = j D{i)p{J^) (4.20) 

the two particle distribution is 

p(a„, r„, a,, v,) = j D{^jD{^,)p{J^a: ^fe) (4.21) 
and in the same way for the truncated distributions, 

PT{aa,ra,ab,n)= J D ) D pt{T a, T^) (4.22) 

At this point, one has to make a fundamental observation. Although the system of 
filaments can be treated with the rules of classical statistical mechanics, it differs from a 
genuine assembly of classical charged random wires (distributed according to the Gaussian 
measure -D(^) ) on an important point: the interaction (4-16) inherited from the quantum 
mechanical nature of the charges is not the standard electrostatic potential between two closed 
wires of shapes and £2- The latter potential, denoted ^1(^1,^2), is 

Ki(^i, ^2) = / dsi f ds2V{ri + KMsi) - r2 - KMS2)) (4.23) 
Jo Jo 

since every element of charge ^i{si)dsi carried by the first filament has to interact with every 
element $,2{s2)ds2 of the other by the Coulomb law, and the corresponding total electrostatic 
energy is 

k 

t/ei(J^i, . . . , J^fe)) = Yl ea,e,. KK^i, JF,) (4.24) 
i<j=i 

Although the filaments are quantum mechanical in origin, we will use this subscript ci to 
denote the filament system with interaction Ud- Replacing the interaction (4.16) by (4.23) 
would indeed lead to a genuine system of charged filaments obeying exponential Debye 
screening in the appropriate low density or high temperature regime, according to the the- 
orems quoted in chapter II and the discussion in section IV. F. Hence different behaviors 
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between quantTim and classical charges in thermal equilibrium must be traced back to the 
difference between the interactions (4.16) and (4.23) 

W{T^,T2) = e«,e«,(y(J^i,.F2) - V,^{Tr,T2)) 



\ dsi 
Jo Jo 



ds2{6{si - S2) - 1) 



xy(ri + Xa.^Si) - r2 - Xa,^2is2)) (4.25) 

The large distance behavior of W{J^i, JF2) is obtained from the multipolar expansion of the 
potential V written in the form 

V{r, + A„,^i(si) - r2 - \aA2is2)) = V{r^ - r^) + M,{s,)V{r, - r^) 

+M2{s2)V{ri - ra) + Mi(si)M2(s2)V^(ri - ra) (4.26) 

where Mi{s) are the multipolar differential operators 



M,{s)^Y.^^^^M±lllt^ ^ = l,2 (4.27) 



fe=i 



In the decomposition (4.26), the three first terms correspond to the charge-charge and charge- 
multipole interactions, whereas the last term represents the rest of the multipole-multipole 
interactions. Since 



/ dsi{5{si - S2) - 1) = / ds2{5{si - S2) - 1) = 
Jo Jo 



(4.28) 



we see that only this latter term contributes to W 



P^(^i,^2) = e„,e„, / ds^ f ds2{5{si-S2)-l)M,{si)M2{s2)V{r,-r2) (4.29) 
Jo Jo 

The dominant large distance behavior of W 

W{Ti,T2) ~ eaie«2 / dsi / ds2{5{si - S2) - 1) 
Jo Jo 

X(AaiCl(Sl) ■ VrJ(A„2^2(^2) ' V r^)V - Y2) , | Ti - r2 ^ OO 

(4.30) 

is typically a dipole-dipole potential, decaying as | ri — Y2 |~^, due to the dipole elements 
eaiAai^i(si)(isi and eo2A„2^i(s2)c?S2 associated with two filaments of arbitrary shapes. It will 
be shown (chapter V) that the charge-charge and charge-multipole interactions have a fast 
screening in the quantum gas, but not the multipole-multipole interaction. This long range 
term will precisely be the cause for the breakdown of exponential screening in the quantum 
system for all values of the density and of the temperature. For the sake of pedagogy, we 
illustrate the above assertions in a simplified model involving only two quantum mechanical 
charges. This model presents already all the basic mechanisms which will be at work in the 
general case. 
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C. Quantum fluctuations destroy exponential screening 



We consider the simplified situation where all particles are classical but two. For this it 
suffices to single out two filaments Ta — (eta, r^, Th — (ctb, r;,, ^f,) corresponding to two 
particles with non vanishing de Broglie lengths Aa, and to set the de Broglie lengths of 
all the other particles equal to zero. Then the phase space variable of a classical particle 
reduces to (a, r) and its interaction with the quantum particle Ta to 

V{T,, r) = dsV{Y, + A„4(s) - r) = V,,{T,, r) (4.31) 

JO 

which is identical to the classical potential (4.23) with one A set equal to 0, whereas the 
two quantum particles still interact with the potential (4.16). Here the quantum mechanical 
nature of the two specified particles will be treated non-perturbatively. The model is semi- 
classical only in the sense that their surrounding medium is now a classical plasma at density 
p and inverse temperature 

As in (1.17) we denote by a; = (q;i, ri, . . . , «„, r„) a configuration of the classical particles 
and by c(r,a;) = YTj=\ ^afik^ — ^j) the corresponding microscopic charge density. For this 
system, the total interaction (4.15) reads 

U [Ta, Tb-, Lo) = e^^e^^V{Ta, Ti) + e^^ ^ dvV{Ta, r)c(r, a;) 

+ y" dxV[Tb, r)c(r, uj) + U^{u) (4.32) 

where V^iuS) is the Coulomb energy of the classical particles. 

The distribution of two filaments Ta and Th in thermal equilibrium with the classical 
Coulomb gas enclosed in a finite volume A is defined by 

p{Ta,Tb) = i / da;exp(-/3[/(J^„,jr,,a;)) (4.33) 

In (4.33), /^c?a; means (canonical or grand canonical) summation on the configurations uj 
of the classical particles in A and So = (ic<j exp(— /5f/o(^)) is the corresponding partition 
function (in this definition, there is no contribution of coincident filaments). For a single 
filament we define in the same way^^ 

p{T) = ^ [ duex.p{-pU{T,uj)) (4.34) 

with 

U{T, uj) = eo, I dr V{T, r)c(r, a;) + Uo{uj) (4.35) 



/3 ^ Inp(^) and — /3 ^ In p{Ta, Tt,) are the excess free energies when one or two filaments are 
embedded in the classical plasma. 
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Let us also introduce the auxiliary genuine classical distributions of charged filaments 
Pc\{^a,^b), Pc\{^) defined in the same way as (4.34, 4.35), but with V and U replaced 
by Vc\ and U^h Comparing now the quantum and this classical system, we observe that in 
view of (4.31) 

U{T, u) = C/ei(^, UJ) , p{T) = pei(^) (4.36) 
but taking into account the definitions (4.24) and (4.25) 

U {Ta, ^6, a;) = W{Ta, T^,) + U^x{Ta, T^, a;) 

P(^a, ^b) = exp {-mkJ'a. r^)) Pcl(^a, ^b) (4.37) 

We consider now these distributions in the thermodynamic limit, keeping the same no- 
tations. Because of translational invariance, p{T) = Pa{$) is independent of r and p{J-'a,^b) 
depends on |ra — r^l. In the Debye regime, we know that Pc\{^a,^b) clusters exponentially 
fast as |ra — r^l — >• oo (with ^^j, fixed) 

+ O ( exp ( -C inf ^ Jra + KaLisi) - n - \ 

\ \ 0<S1,S2<1 

(4.38) 

Thus the asymptotic behavior of p{J^a:^b) in (4.37) is governed by that of W{Ta-,Ti,), i.e. 
according to (4.30), p{Ta- J^b) clusters as Ir^ — rft|"^. The two-particle distribution is obtained 
by averaging p{Ta, ^b) on the shapes of the filaments. Using (4.36, 4.38) and expanding the 
factor exp {—(iW{TaiJ'h)) yields 

p(r„,r,) = y D{i:)D{i^)^y.v{-m{J'a.Tb))p.ATa.Tb) 

^PaPb-a j D{OD{i,)W{Ta,r,)pa{Up,{i,) 

+ j D{^,)D{QW\j^a, J'b)Pa{ia)Pb{ib) + • • • + O (exp(-C|r„ - v,\)) (4.39) 

where Pa = p{aa,ra). The remainder in (4.39) still gives an exponentially decreasing contri- 
bution since the probability of long filaments is Gaussian small. Let us show that the term 
hnear in W in (4.39) does not contribute to the asymptotic behavior of p(ra,rf,). Indeed 
one notes that the measure -D(^) is invariant under reflections and rotations of the filament 
in three-dimensional space (see (4.7)). The same is true for p(^) in an homogeneous phase 
of the classical plasma. Hence introducing the multipolar expansion of W with the help of 
(4.27), (4.29) one sees the occurrence of the ^-averages 

J D{c)D{^,)piapi^b) 

[XaU^l) ■ '^rf[XbUs2) ■ VrJ V(r„ - T,) k, I > 1 (4.40) 
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In view of these symmetry considerations, the terms with k or / odd vanish, and the terms 

k+l 

with both k and / even are necessarily proportional to (V^) ^ V{r),r = — r^, which is 
a rapidly decreasing function since = 0,r 7^ 0. Since all moments of the Gaussian 

measure -D(^) are finite, one concludes that the W^-linear term in (4.39) decays faster than 
any power of |ra — r^l. The leading behavior is therefore determined by the quadratic term. 
According to (4.39) and (4.30) one finds for the normahzed truncated distribution 



PaPb PaPb 
2paPb 



D{L)mbW{J^a,J'b)Pa{L)Pb{^b)+0{W') 



where B{P,p) > can be expressed, from (4.30), as the Brownian integral of a positive 
functional, depending on the temperature and density of the surrounding plasma. The 
correction | — rf, in (4.41) arises from the square of the dipole-quadrupole interaction. 
Note that the result (4.41) holds in any screening phase of the classical plasma and is non- 
perturbative in h. 

In the classical limit or in the high temperature limit, one finds that the dominant 
contribution to B{P,p) is the same (Alastuey & Martin 1989) 

B{P,p)c^h'^^, h^O or P^O (4.42) 
240 maTrib 

The role of the quantum fluctuations is clearly demonstrated: the dipole-like interaction 
W, which originates from these intrinsic quantum fluctuations, can deflnitely not be shielded 
by the classical charges, even if the latter are in the Debye phase. Nevertheless a partial 
screening remains: the monopole is still screened and this, together with rotation invariance, 
reduces the bare Coulomb decay to r~^. It will be shown in section IV. E as well in 
chapters V and VI that these facts remain true in the full quantum gas. 



D. Origin of the van der Waals forces 

One can interpret the result (4.41) by saying that there exists an effective potential 
$i(r) between two individual quantum charges defined by ^^^'""'^'"''^ — e"^*^*^!''""'"''!^ — 1. It is 
attractive at large distance and behaves as 

/?$i(r) ~ r ^ 00 (4.43) 

One may wonder if this potential between individual quantum charges is related to the van 
der Waals potential, also decaying as r~^. The van der Waals forces are usually computed for 
preformed atoms and molecules in empty space, and are attributed to the dipole fiuctuations 
in these bound entities. 
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To allow for the possibility of quantum mechanical binding in simple terms, we can 
generalize the model of section IV. C by specifying four particles, say two electron-proton 
(e-p) pairs, in equilibrium with a classical plasma (Martin 1996). The information about 
the correlations between these four particles is contained in the fully truncated four point 
correlation function pri^a, i"6, Ri, 1^2) where Tq, designate the coordinates of the electrons 
and Ri, R2 those of the protons. This correlation is defined by the natural generalization of 
the expressions (4.33)-(4.37) to four particles immersed in the plasma. Wc arc interested in 
the amount of correlation that exists between the two (e-p) pairs as their centers of mass arc 
taken far apart, but ignoring the relative positions of the particles. Assuming for simplicity 
that the two protons are infinitely heavy (i.e. classical), the center of mass of the pairs 
coincide with Ri,R2 and the desired normalized correlation is 

^ ' 2/cir„pr(ra,Ri)/cirbpr(rfe,R2)' I 1 21 y ) 

The factor 1/2 takes into account that the electrons, treated here as identical particles, can 
be found in the neighborhood of either one of the protons (the contribution of coincident 
particles is omitted in (4.44)^^). We define the effective (temperature and density dependent) 
potential ^2{R) between the two pairs by 

F(i?) = e-^*^^'^) - 1 ~ -/5$2(i?), R^oo (4.45) 

In contrast to the standard calculation of van der Waals forces, $2 incorporates now the 
effects of the thermal fluctuations of the electrons and of the screening provided by the 
medium. 

As in (4.39), the asymptotic form of $2(-R) is determined by the quadratic term in the 
quantum potential W , namely 



-/3$2(i?)^ j D{i,)D{i,)W\Ta.J'b) 



where now -D(^a), -D(^j) are dressed Brownian measures involving the single (e-p) pair 
correlation. One finds 

/3$2(i?) ^ it: ^00, C(/3,p)>0 (4.46) 

There are two situations of interest, the high temperature limit (/? ^ 0, p fixed) and the 
atomic limit. The atomic limit is obtained by simultaneously lowering the temperature {j3 — > 
00) to favor binding and reducing the density at an exponential rate (p ~ e~^l^ ^0, (5 > 0) 
to get independent atoms. In this limit (e-p) pairs form noninteracting hydrogen atoms in 
their ground states (the atomic limit will be discussed in chapter VII for the full quantum 
gas). 



^^As before, exchange effects due to the Fermi statistics of the electrons are neglected. Exchange 
correlations have a short range. 
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In the high temperature hmit, one recovers the result (4.42), i.e. 




(4.47) 



240m2' 



As /3 — > 0, the electrons tend to be fully ionized and the correlation becomes that found 
between individual charges. 

In the atomic hmit, for 5 > not too large, one has 



where Cyan der Waais < is the usual coefficient found in text books from a second order 
perturbation calculation of the residual Coulomb interaction between two hydrogen atoms 
in their ground states. 

A similar analysis has now been carried out in the full quantum mechanical electron- 
proton gas with the diagrammatic methods of chapter V. One finds indeed that the coeffi- 
cient of the tail of the proton- proton correlation approaches the value Cyan der Waais in the 
scaling limit where hydrogen atoms form (Alastuey, Cornu & Martin 1998). In the medium, 
density and temperature effects modify the amplitude of the usual van der Waals potential, 
but not its range. In particular, the free charges that are always present in the medium 
only partially screen the dipole interaction between atoms, contrary to a naive mean field 
description that would lead to an exponential screening of these interactions. The partial 
screening of protons is due in part to collective effects, and also to individual electrons that 
will bind with the proton. As the temperature goes to zero the former effects vanish, and 
the screening is eventually entirely due to the formation of a neutral bound state. Hence 
the atomic limit coincides with a zero density calculation with two preformed neutral atoms 
in empty space. In fact, the coefficient C(/3, p) interpolates continuously between the high 
temperature and atomic regimes. There is no qualitative difference between the mechanism 
producing the van der Waals forces between bound entities and that causing the tail 
(4.41) between individual charges. Both have their common origin in the same basic quan- 
tum fluctuations put into evidence in the preceding section. One sees that the difference is 
only quantitative. In the atomic regime, the coefficient C(/3, p) has a non-zero value close 
to iC'vander Waaisl whcrcas it bccomcs vauishiugly small at high temperature. 

E. Semi-classical analysis of Coulombic correlations 

We come now to the full quantum mechanical gas, but treated semi-classically as a whole. 
We use again the classical gas as a reference system by sphtting the interaction between the 
filaments (4.15) into 



lim /3~^C(/?, p) = -a 



van der Waals 



(4.48) 



U (J^i, ...,J^n) = Uo{ai, ri, . . . , a 



I"n) + Um{^ 1, ■ ■ ■ ,^n) 



(4.49) 



where 



n 




(4.50) 



i<j 
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is the Coulomb energy of n classical point particles and Um is the residual multipolar interac- 
tion. It is defined as the sum of two-body interactions Vm{J^i, ^2) = ^2) — ^(^i — 
obtained from (4.16) by removing the charge-charge potential: 

Kn(-^i,-^2)= / rfs(Mi(s)+M2(s))V(ri-r2) 

+ / dsMi(s)M2(s)l^(ri-r2) (4.51) 

n 

C/„(^i, . . . , ^„) = ec^fi^ymi^h ^3) (4-52) 

Since the Planck constant occurs only through the de Broghe lengths Aa, it is clear that 
Vm and Um are power series in h starting with a linear term (see (4.5) and (4.27)). Thus 
considering ^ as a small parameter, we are left with the standard problem of perturbing 
the classical equilibrium state (in the framework of the system of filaments) by the small 
modification (4.52) of the potential. This will generate formally power series in h for the 
distributions of filaments, and in turn, after averaging on their shapes, for the truncated 
distributions of the quantum charges 

PT(ai,ri, . . . ,Q;„,r„) = ^ti" pf\ai,Yi, . . . ,Q;„,r„) (4.53) 

fe>0 

We shall not give here the rather complex algorithm for calculating the coefficients pj, of 
this expansion, but only describe some qualitative aspects and state results. 

From now on we assume that the classical reference system is in the Debye phase: its 
correlations Pc1,t(q;i, ri, . . . ,a„,r„) cluster exponentially fast and obey the multipolar sum 
rule (2.16). By expanding around the classical gas, each involves combinations of convo- 
lutions of classical correlations with the multipole potentials coming from (4.51), so that the 
prp are expressed in terms of the supposedly known correlations pci,T of the reference system. 
The main point is that, although the classical correlations have an exponential fall-off, the 
law of formation of the pip ^ is such that some of these convolutions decay only algebraically 
(there are two exceptions, namely if the potential is convoluted with the total charge density, 
then multipolar sum rules insure a fast decay, or if derivatives of the Coulomb potential are 
contracted to form a Laplacian V^V^(r)). 

The results are as follows (Alastuey & Martin 1988, Alastuey & Martin 1989). Obviously 

PT\ai, ri, . . . , an, r„) = Pc1,t(q;i, ri, . . . , r„) (4.54) 

pip\ai,ri, . . . ,an,rn) = 0, k odd (4.55) 

This is because, according to (4.27), odd powers of h in the expansion (4.53) correspond to 
odd moments of filaments. Such odd moments have zero Gaussian average. 

The coefficients (Q;i,ri, . . . ,an,rn) are decreasing faster than any inverse power for 
large spatial separations of the arguments (for the two point function see (Jancovici 1976)). 
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At the order fi^ the truncated two particle correlation behaves as 

(p(ai,ri)p(a2,r2))? = "^i^ "2, 

2 f(3^^^ 



5 (Q / / - ^2 + r - rOGlr, aOGlr', ^2) (4.56) 

In (4.56) /(r) is the square of the dipole-dipole potential 

3 

/(r) = J2 id,d^yi^)f (4-57) 

and 

^(r, a,) = 5^ ^ (p(a, 0)p(a,, r))ei,r, i = 1, 2 (4.58) 



where (• • •)ci means equilibrium average with respect to the classical reference system. Since 
/(r) ~ |r| — > 00, the leading contribution to (4.56) is 

(p(ai,rOp(«2,r2))?) ~ J^^^e^ (/ dvG{v,a,)^ {^j dvG{v,a,)^ (4.59) 

If one examines the correlation of a particle with the total charge density c(r), one 
must replace p{a,r) in one factor G by the charge c(r). Then the coefficient of |ri — r2|~^ 
vanishes in (4.59) because of the charge sum rule (2.16), and one has to carry the asymptotic 
expansion of /(r) further on. For the particle-charge correlation, one finds 

{p(ai, ri)c{a2, r2))T^ = ^ e^^P? ("i, ri, a2, ra) 



02 



1 

Iri — r2|® 16 



I dvG(v,a,)^ Q I dr\r\'K{r)^ 



(4.60) 



and for the charge-charge correlation 

(c(ai,ri)c(a2,r2))T^ = ^ eaiea2PT\ai,ri,a2,r2) 



ai,a2 



1 7/?Vl f dr\r\^K{r)] (4.61) 



with 

^(r)=E— (4.62) 
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In fact, as discussed in sections VI. A. 2 and VI. A. 3, the formulae (4.59), (4.60) and (4.61) 
become asymptotically exact for the quantum Coulomb gas in the low density limit. In this 
limit the functions G{r) and K{r) are given by their Debye-Hiickel approximation^^. 

The formulae simplify in the case of the jellium when one has only one species of charges. 
Then there is no distinction between particle and charge correlations up to a factor e 
(eG(r) = K{r) — ^(c(0)c(r))cLT), so (4.61) applies in all cases. Furthermore, in view 
of the second moment condition (1.33), the correlations at order in the jellium behaves 
as 



For the jellium, (Gombert & Leger 1994) have given the complete formula for (ri, r2) at 
the order including the short range contributions. 

Correlations involving more than two particles also have slow decays at order fi^. For 
instance p^''(«i, ri, 0:2, r2, as, rs) decays also as |r|~^ when the particle 1 is sent to infinity, 
the other being kept fixed. If the three particles are taken simultaneously far apart (say, 

— i"2| = r, |ri — rsl = ar, |r2 — rsl = br, r ^ 00), this function decays as r~^. Moreover, 
the four point function p^'*(q;i, ri, a2, r2, as, rs, a4, r4) decays only as when two pairs of 
particles are sent far away (say, |ri — r2|, jrs — fixed, r = |ri — rsj — > 00). This slow 
decay is due to the breaking of rotational invariance around the points ri and r2 when the 
directions ri — r2 and rs — are fixed. It is precisely because of rotation invariance that 
the pure dipolar decay is reduced to in the two point function, as in the model of 
section IV. C. 

We have reported on the semi-classical tails (4.59)-(4.61) in some details with the pur- 
pose of illustrating the following points. There is definitely no exponential screening in 
the quantum gas, even in the high temperature phase and what ever the density may be. 
Moreover, in contrast with a fiuid with short range forces, various types of decay laws have 
to be expected depending on what type of observables are concerned, and what cluster of 
particles are taken far apart. The occurrence of these different types of decay are due to 
the underlying multipolar forces displayed in the filament formalism in conjunction with the 
effect of screening sum rules (2.16). 

We finally discuss the possible observable implications of these findings. Since the calcu- 
lations pertain to the semi-classical regime, we only consider real systems under conditions 
that the quantum effects are small. We determine the crossover distance ro above which 
quantum algebraic tails dominate the classical behavior by equating (4.63) to the usual 
Debye law 



^''Note that the exchange contribTitions vanish exponentially fast with h, so leading contributions 
in h can be obtained from Maxweh-Boltzmann statistics. 




(4.63) 




(4.64) 
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The classical exponential is normalized in order to obey the charge sum rule. For instance, 
for sodium chloride at room temperature, one finds Tq ~ 60 l^- For a white dwarf under 
prevailing conditions, one has Tq ~ 40 Id- In both cases, the quantum effects on the 
clustering turn out to be very small and the use of the Debye potential is still legitimate 
from a quantitative point of view. However if we apply crudely (4.64) (disregarding its semi- 
classical nature) to the electrons in a metal (replacing Id by the Thomas- Fermi screening 
length Xtf), one finds tq of the order of Xtf- This indicates that the use of the exponential 
potential might be less reliable in this case. 

F. Breakdown of exponential screening in the Quantum Sine-Gordon representation 

In this section we show how the Sine-Gordon transformation of section II. C extends 
to the quantum mechanical system (Frohlich & Park 1978). We saw in (2.24) that the 
classical partition function is a Gaussian superposition of ideal gas partition functions in 
external fields i(j){r). The same holds for quantum mechanics except that one must use time 
dependent fields i0(r,s). 

After the Sine-Gordon transformation is discussed we use it to illustrate the natural 
appearance of dipole forces in the quantum system and make some comparisons with classical 
systems of particles that simultaneously have charge and dipole moments. 

We consider a quantum Coulomb system as in (1.36), but replace the Coulomb interaction 
(1.1) by 

1 ^ 

UocL = - J^e„.e«^.Kx.,L(ri - r^) (4.65) 

where Voo,l is the potential defined in (1.9) whose singularity at the origin is smoothed out on 
length scale L. Note also that self-energies are included in t/oo,L- This choice of interaction 
simplifies the Sine-Gordon transformation. It should be possible to reduce the 1/r to this 
case by a Mayer expansion as was done in the classical case in section II. D because the 
Mayer expansion is known to converge (Brydges & Federbush 1976, Brydges & Federbush 
1977). 

From (1.40) with external field cj) set to zero and using the Trotter product formula 
Ea = TrP exp {-(3 {Ha,n - M • N)) 

= TrP hm ff (e-^lHl.~li-m^-^u^,L\ (4.56) 
1 

where the H^ j^ is the free Hamiltonian and /it • N = 'Ylial^oiNa- We insert 

Each factor e~n'^°°'^ requires its own auxiliary field 0i(r), i = 1, . . . ,n, so that the right 
hand side of (4.66) contains 
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/n 
dl^M ■ ■ ■ d^iL{<Pn) TrPj] (e-^^A.-e-('^/")'''^/'^'-^W^^W) (4.68) 
1 

By scaling the fields (f)j n^^^'^(l)j and uniting them into one time dependent field 0(r, s) = 
0j(r) when s G [ '■-'~^-' , we obtain 

5a = rf[°'VL(</')5ideai,A(/?'/'i</') (4-69) 
where J d)-^'^^ijLL{4>) is the Gaussian defined by 

j 0(r, s)0(r', s') = 5(s - s')Ko,L(r - r'); < s < 1 (4.70) 

and 

Sideai,A(0) = IVP exp ^ ,v " ■ N + ^ dr c(r)0(r, s)^ rfs^ (4.71) 

is the ideal gas partition function for a timc-dcpcndent external field and the exponential 
is time-ordered. In the case of Boltzmann statistics, which means omitting the sum over 
permutations in the projection P, the trace factors over particles so that 



5ideal,Bolt,A(/30) " XI II 7V~! 

{Nc} a 

xexp(^-pj^ (^Hl^^ - ^^■'N + J dr c(r)0(r, s)^ ds^ 

= exp Tre-^-''o ds[h^i')-n.]^ ^4 72) 

where ha{s) = ^-{—A) + ea4>{r,s) is the one-particle Hamiltonian for species a in the 
time-dependent external field and the exponential is still time-ordered. Without making the 
simplification of Boltzmann statistics one has, instead. 
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Sideai,A(/30) = exp J] J] n^T^^-PlSMh^is)-,.] (4 73) 

V a g=l ^ / 

where 0(r, s) is extended periodically in s. This formula will be proved in section V.A.4. 
Notice that this quantum ideal gas is invariant under 

(/.(r,r)^(/.(r,r) + /(r) (4.74) 

where / is any function independent of r such that j^drfij) = 0. The Gaussian rft^'^'yU is 
also formally invariant under this transformation. This "Goldstone mode" is a signal that 
there are long range correlations in the field, whereas in the classical system correlations 
are exponentially decaying. 
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By the Feynman-Kac formula, 



Sideai,Boit,A(/?0) = exp (^j dT z{T)i 

4>{T)= [ <P{r + Xai{s),s)ds (4.75) 
Jo 

where J-', were defined in (4.11, 4.13). There is a similar formula for the ideal gas (4.73) 
without the Boltzmann statistics simplification and it is given in section V.A.4. It could be 
used in the following considerations, but would not change the discussion in any important 
way. Furthermore we set the spins = 0. 

When Aq = the phase space integration dT reduces to J dr and (4.75) reduces to 
the classical ideal gas. This suggests (Brydgcs & Keller 1994a) the following model. We 
replace the integration J d!F over all Wiener paths by integration concentrated on just one 
kind of path which oscillates about the initial point by a distance 0(A,J (the size of the wave 
packet) in a random direction: let da{e) be a spherically symmetric integration on vectors 
e. Then J dT is replaced by 

^ j dr ! da{e) (4.76) 

ce 

and each (r, e) labels a path: r + Xai{s) with ^(s) = esin(27rs). We chose sin(27rs) because 
it is orthogonal to 1. 

The consequence of this and another minor approximation (Brydges & Keller 1994a) is 
that the dependence of (4.72) on 0(r, s) is only through two modes 



' (j){r,s)ds; 02(r) = / 0(r, s) sin(27rs) (is 
Jo 



Since 0i, 02 are Gaussian and J d^^'^^HL(pi{r)(j)j{r') = Voo,l(i" — the integral J d^'^'^^HL 

can be replaced by J diiL{(f)i) J diiL{(f)2) and the partition function becomes 

5a = y" dlJiL{(t>2) J diiL{(pi) 

xexp| J]e'^'^" Jdr J da{e) e'^'^"'^'+'W2)'^"^"'^-^A (4.77) 



Notice that if da{e) is set to 5{e) we revert to the classical Coulomb gas. If 0i is set to 
zero then by reversing the Sine-Gordon transformation we obtain the partition function of 
a classical dipole gas with dipole moments e distributed according to da. Furthermore part 
of the Goldstone mode (4.73) survives as 

02(r) ^ 02(r) + const (4.78) 

If the fields (p = (pi = (j)2 are set equal, this becomes exactly the Sine-Gordon transform of a 
classical gas of particles that are simultaneously carrying charge and a dipole moment. But 
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in this case the symmetry (4-'78) no longer holds and the methods of chapter II prove that 
correlations decay exponentially, near the Debye-Hiickel limit. 

Indeed consider the Dcbyc-Hiickel approximation wherein the exponent in (4.77) is re- 
placed by the quadratic approximation 

Q^-IJ^'^"'" I J (/3^/'0 + (/3/2)^/^A«e-V0)' 

a 

= J + l(^c.V<l>A (4.79) 

The Gaussian diiLexp{—Q) has a covariance whose Fourier transform is 

47r (k^ + L^k^ + + ak^) (4.80) 

where — An fJe^i^" and a — An /^A^e^'*" . There is no pole for any real k which 
indicates exponential decay of all correlations. It has been wrongly said (Martin 1988, 
section II G) that classical ions with structure might have weaker screening properties than 
pure charges. The present discussion shows that such systems have Debye screening near 
the Debye-Hiickel limit. 

The model defined by (4.76) and (4.77) exemplifies again that the breakdown of expo- 
nential screening is really due to quantum mechanics, manifested here by the existence of 
the two fields 0i and 02 (and more generally by the time-dependent field 0(r, s) in the Sine- 
Gordon representation). Reconsider indeed that (f)i ^ 02 ai'c independent Gaussian fields. 
The quadratic approximation is misleading here because the cross-terms coupling (pi, 02 are 
lost in this approximation but they are present at cubic and higher order in (4.77). In 
(Brydges & Keller 1994a) it is argued that these higher order couplings allow the long range 
correlations of 02 to creep into the correlations of both 0i and 02. This gives an alternative 
but more qualitative derivation of the no-screening results. The massive zero frequency field 
01 mediates (parts of) the electrostatic potential defined in (4.23) whereas the non-zero 
frequency field 02 mediates (part of) the non-electrostatic part W defined in (4.25) that 
cannot be screened. 

All these discussions show that screening is destroyed because lack of commutativity 
of momentum and position bring time dependent phenomena in through the back door 
(e.g., the imaginary time ordering mentioned below (4.71)). It is in fact easy to prove that 
imaginary time dependent observables are not screened (Brydges & Seller 1986). 

We should also mention that another class of semiclassical models has been introduced, 
to deal with polarizability effects in dielectrics (H0ye & Stell 1981, Brevik & Il0ye 1988). 
In these models, internal atomic degrees of freedom are treated quantum mechanically (e.g. 
harmonically bound pairs of charges) while the center of mass of atoms behave in a classical 
manner with pair-wise dipolar forces. These models capture physically interesting effects but 
do not incorporate all the effects stemming from the Coulomb potential (binding, ionization, 
collective screening) in an equally consistent and fundamental footing. Hopefully the loop 
formalism developed in chapter V does provide a way to treat all these effects in a fully 
coherent way. 
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V. THE GAS OF CHARGED LOOPS 



A. The statistical mechanics of loops 

1. The magic formula 

We are ready to develop the general Feynman-Kac representation of the Coulomb gas. 
The first step (section V.A) is to include the Fermi statistics of the electrons and the Bose 
or Fermi statistics of the nuclei. This can be done with some combinatorial effort by a 
rearrangement of the sums in the partition function (1.40), leading to the "magic formula" 
(5.16), that has a marvelous effect: as in the preceding section, the gas appears as a classical- 
like assembly of extended objects, the charged loops, and the diagrammatic techniques of 
classical statistical mechanics will apply. The magic formula has been known at least implic- 
itly for a long time, starting from (Ginibre 1965), see also (Ginibre 1971). A combinatorial 
proof of it, following the lines of (Ginibre 1965) can be found in (Cornu 1996a). It has 
been also retrieved recently in (H0ye & Stell 1994). In subsection V.A. 4 we give the main 
lines of a derivation that does not make use of combinatorial arguments. From now on, 
V{ri — r2) = is the exact Coulomb potential on all scales: a short range regulariza- 

tion is no more needed since the latter is henceforth provided by quantum mechanics that 
will be now fully taken into account. 

The rest of the chapter is devoted to the development of a Mayer diagrammatic for this 
system, following mainly (Cornu 1996a). In section V.B the monopole and multipole inter- 
actions are singled out in close analogy to what was done in section IV. B, in order to isolate 
the part W{Ci,C2) (5.36) that will ultimately be responsible for the non-exponential de- 
cay. Then, following the classical treatment of Coulomb divergences, one sums the Coulomb 
chains, providing thus the screening length k,~^ (5.56) (a comparison of this screening length 
with that of the standard RPA theory is postponed to section VI. A. 4). Finally in section 
V.B. 3 the class of resummed diagrams is reorganized in prototype graphs involving four 
types of bonds. Prototype graphs are finite, decay at least as at large separation and 
will serve as the starting point of low density calculations. The content of this chapter as 
well as the next one is mathematically formal since nothing is known about the asymptotic 
nature of the infinite series of prototype diagrams. 

We now come back to the partition function (1.40). We first represent the matrix element 
in (1.40) by a Feynman-Kac integral. Since i^A^A? is independent of the spins, the spin scalar 
product YliLi < '^ap{i) I'^^i ^ factorizes out in this matrix element, and according to (4.6) its 
positional part is represented by the functional integral 

^ / 1 \ 3/2 /I 12 



/ 1 \ ' / 

< {rp{i)}\ exp (-/3i/A,iv) \{ri} >= JJ ( 1 exp i 



2Ai 



^ / n ^(^i) Yl ^"'^"^ ^ dsV{ri^p^i)(s) - rj-,pO)(s)) j (5.1) 

We have kept the notation (4.4) for a path with extremities rj^p(j)(s = 0) = and rj_p(j)(s = 
1) = 

V-1 



ri,p{i){s) = (1 - s)ri + sTp^i) + Xai^iis) 



(5.2) 



The same remark made after (4.9) concerning the implementation of the boundary conditions 
remains vahd: all paths are constrained to stay inside A. 

The main observation to be made at this point is that by a rearrangement of the sums, 
the partition function (1.40) can be formally written in the same classical form as (4.12), 
provided that one introduces a suitably enlarged phase space. This rearrangement exploits 
the fact that any permutation Pa can uniquely be decomposed into a product of cycles. A 
cycle of length q, 1 < q < N^, is a, subset of the particle indices that are permuted 
among themselves in a cyclic way under pa, for instance 

(l,2,...,g-l,g)^(2,3,...,g,l) (5.3) 

The objects constituting the enlarged phase space are precisely associated with such cycles 
in the following way. The permutations p in (5.1) involve now open paths (except for the 
trivial one). However the set of open paths belonging to a cycle of q elements occurring in 
p can be joined together to form a closed loop C (sometimes also called a polymer) with q 
particle coordinates. The loop C corresponding to the cycle (5.3) is made of the collection 
of paths {rfc,ik+i(s)}fc=i,...,q, with g + 1 identified to 1 (see Figure 1). 
Thus a loop is specified by the variables 

{a, q, {rk}k=i,...,q, {^k}k=i,...,q) (5.4) 

It is convenient to unite the q filaments constituting the loop into a single path parameterized 
by a time parameter s running from to g setting 

R(s) = rfc,jk+i(s), k^[s] + l, 0<s<q (5.5) 

where [s] = integer part of s and s — s(mod 1) = s — [s]. One can then locate the loop at 
R by selecting the position of one of the particles, say ri — R, and write 

R(s) = R + A„X(s), X(0)=X(g) = (5.6) 

Then the particles are located at the points 

rfc = R + A„Xfe, Xfc = X(A;-l) (5.7) 

In (5.6), X(s) is the shape of the loop £ at R, as ^(s) was describing the shape of a filament 
at r: it is again a dimensionless Brownian bridge starting and returning to the origin within 
the time q. Its normalized Gaussian measure -D(X) is the composition of that of the q open 
filaments that constitute the loop (see (4.6)) 

D(X) = i2.qf^ n "^^"^^f^yr^'^'^ ^X, . . . dX,Di^,) . . . DiQ (5.8) 

where one sets Xi = X^+i = 0. One can calculate its covariance from the definitions 
(5.5)-(5.6) 
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J D{-X)X^{si)X,{s2) = S^^q 



mm 



£l £2 \ _ Si£2 

q' q J q q^ 



(5.9) 



Thus one can think of a loop C either as a set of particles coordinates (as in (5.4)) or again 
as a single Brownian path, setting 

£=(«, q, R, X(s), 0< s <g) (5.10) 

It is now possible, generalizing also (4. 14)- (4. 16), to define activities and interactions of 
loops. The interaction between two loops £i, i = 1,2, is the sum of those of the particles 
constituting the two loops 



V{Ci,C2) = J2J2 dsV {ri^kiM+i{s) - ^2,k2M+i{s)) 

fel = lfe2 = l"^° 

rgi r<i2 

= / dSi / ciS2(^(Si-52)V^(Rl + A«,Xi(si)-R2-A„,X2(s2)) (5.11) 

Jo Jo 
and for n loops 

n 

U{£,,...,£n) = Yl ea,eay{C^,CJ) (5.12) 

l=i<j 

The activity z{C) of a loop incorporates the effects of quantum statistics, the spin degeneracy 
as well as the internal interaction U{C) of the particles in the same loop 

ziC) = (2.. + 1)^^ ^^i^ exp(-/3C/(£)) (5.13) 



with 



^^lo "^^^ "^^^ ~ '^'''^''''^^ ~ h)V{Xa{X{si) - X(s2))) 



(5.14) 



Clearly, all the above quantities are invariant under cyclic permutations of the particles in 
the same loop. Notice also that particles in the same loop have the same charge; thus their 
mutual Coulomb interactions are positive and this implies the bound 

I -mi < ^^" + ^ exp(/3/x^g) 

^' - (27rgA2)3/2 ^^"^^^ 

Considering now the phase space of loops together with the above definitions, it can be 
shown (see section V.A.4) that the partition function can be written in the same form as 
(4.12) {the magic formula) 
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^^ = E-| / lld£iz{£i)exp{-pU{£,,...,£r^)) (5.16) 



n=0 " i=l 

where the phase space integration means here 

S oo 



A 



^^••• = EE idR I D{^)... (5.17) 

a=l 9=1 



We call the system defined by the relations (5.9)-(5.17) the "system of loops"^^. It reduces 
to the system of filaments (i.e. charges with Maxwell-Boltzmann statistics) if all loops with 
more than one particle are omitted. This view of the quantum Coulomb system as a gas of 
loops is well-suited in the low density regime, when loops remain well localized objects. 



2. Loop and particle correlations 

It is clear that the statistical mechanical system of loops defined by the above relations 
has a classical structure. Therefore, introducing the "5-function" which identifies two loops, 

5(£i, A) = 5„i,a.<^,i,5.<^(Rl - R2)5(Xi, X2) (5.18) 

one can define the density of loops p{C) and the two loop distribution function p(£i, C2) by 
the formulae (4.18)-(4.19) with C replacing J^. 

In order to obtain the particle distributions from those of the loops, some care has to be 
exercised because loops are constituted of several particles. For a configuration of particles 
distributed on a set of loops Ci, we can write the particle density and distributions as a 
summation on the loop index i 

P(">r) = ^(5«^,„^(5(r(fc,i) -r) (5.19) 

i k=l 

where, according to (5.7), r^^ j) = Rj + Aa.Xj(/c — 1), k = 1, . . . ,qi, are the positions of the 
particles belonging to the loop i. For 7^ one writes also 

p{aa, ra)p{ab, r^) = 5ai,a.5a„a, E ^(^(^4) " ra)S{r(^k',j) - n) 

ijtj k=l k'=l 

Qi 

+Saa,a, Yl Yl ^^^(^'^) ~ ^a)Sir^k,^i^ - Yh) (5.20) 

i k^k' 



^^Oiir definitions of the weight of a loop (5.8), the activity (5.13) and the self-energy (5.14) differ 
slightly from those found in (Cornu 1996a). It is natural to define the loop measure (5.8) again as 
that of the Brownian bridge process in time [0, q] . 
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In (5.19) and (5.20), k {k') runs on particles belonging to the same loop. The first term in 
(5.20) refers to particles in different loops and the second term to particles within the same 
loop. Then a comparison between the density of loops p{C) = J^i^i^i^ ^) ^^'^ particle 
density (5.19) shows that 

p(a, r) = Y, dR D{X) J] S{r, - r)p{£) (5.21) 

q=l k=l 

Taking the average of the above relation on the ensemble of loops gives the number density 

oo „ 

p{a, r) = ^ W D{-K)p{a, q, r, X) (5.22) 

g=l 

where p(a, g, r, X) = p{C) with C = {a, q, r, X). We have used the fact that p{C) is invariant 
under cyclic permutations of the q particles in the loop. The two particle distributions are 
obtained by averaging the relation (5.20) where we distinguish again between pairs belonging 
to different loops or exchanged in the same loop (r^ 7^ r^) 

p{aa, r„, ab, r^) = p^"^^) («„, r„, a^, r,,) + Sa^^a.P^^''^ (««, fa, fb) (5.23) 

with 

p^''^''\aa,ra,ah,rh) = 

^^QaQb / D{Xa) / D{Xa)p{aa,qa,ra,Xa;ab,qb,rb,Xb) (5.24) 

and 

00 „ „ q 

p^'''\aa, r„, r;.) = ^ / / D(K) ^ S{n - ra)S{n' - r,)p(a;„, q, R, X) (5.25) 
9=2 k^k' 

For a translation invariant state {p{a, g,R, X) independent of R), one has, using (5.7) 

00 „ q 

p^'^\aa, r„, n)^Y. ^(^) E ^(^«„(Xfe - X,0 + r;. - r„)p(a„, q, X) 

q=2 kj^k' 
00 „ q 

q=2 k=2 

The second line follows again from the invariance of p{aa, q, X) under cyclic permutations 
of the particles in the loop. 

One can also write the corresponding similar representations for the truncated distribu- 
tions. 
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3. Loop interactions 

The same fundamental observation made for the system of filaments remains valid here: 
the interaction (5.11) is not the standard electrostatic Coulomb energy of two charged loops, 
which would be 

Ki(A,>C2)=/ dsi rfs2V^(Ri + A„,Xi(si)-R2-A«,X2(s2)) (5.27) 
Jo Jo 

and the discussion following (4.23) can be reproduced word by word. In particular, the 
Coulomb potential has the multipolar expansion 

y(Ri + A„,Xi(si) - R2 - A,,X2(S2)) = V{R, - R2)+ 

Mi(si)y(Ri - R2) + M2{s2)V(Ri - R2) + Mi(si)M2(s2)y(Ri - R2) (5.28) 

with now 

M,(.) = i^^^i^^^i^^^^, i = l,2 (5.29) 

k=l 

and hence the loop interaction can be decomposed into its charge-charge, charge-multipole 
and multipole-multipole components 

V{Ci, C2) = V^{Cu C2) + V,m{Cu C2) + Vmc{Cu C2) + Vmm{Cu C2) (5.30) 

Kc(>Ci, £2) = gig2l^(Ri - R2) (5.31) 

rQ.2 

Vcm{CuC2) = qi / dS2M2(s2)V^(Rl ' R2) (5.32) 

K^c(^l,^2) = 92 r ds,]Vh{s,)V{R, R2) (5.33) 
Jo 

rn r<i2 

Vmm{C^.C2) = / dSi / dS25{s^ - S2)M^{s^)M2{s2)V {H^ - R2) 

Jo Jo 

(5.34) 

The relevant quantity which will again embody the long range effects is the difference 
between the quantum and classical interaction 

W^(£i, £2) = e„,e«,(V(£i, £2) - KK^i, £2)) (5.35) 
It depends only on the multipole-multipole part of the interaction as in (4.29) 

Vr(£i,£2) = eajCaa / dsi I dS2{S{Si - S2) - 1) 
Jo Jo 

xMi(si)M2(s2)V^(Ri - R2) (5.36) 
and has the asymptotic dipolar character (4.30) as |Ri — R2I ^00 

dsi I ds2 



JO 



x{S{h - §2) - l)(Ac,,Xi(si) • VrJ(A«,X2(s2) • Vn,)V{Ki - R2) (5.37) 
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4. Derivation of the magic formula 

We derive first the formula (5.16) in the special case where there is an external field 
(/)(r) and no two-body interaction. In this case the magic formula is essentially the same as 
(4.73). 

It is enough to consider just one species with charge Cq, and spin Sa because the partition 
function factors. In terms of annihilation and creation operators ao-(r), o^l^) many-body 
Hamiltonian for the external field case is 

al{r)haa^{r) dr; K = ^(-^) + (5-38) 

We begin with the averaged particle number 

W = ^lnSA = i-J] I drIY(a:(r)a.(r)e-'^(^«-'^^)) 

N = J2 f aM^air) dr; 5a = Tre-^^^""''^) (5.39) 

The canonical commutation relations for the annihilation and creation operators d^{v), a^(r) 
imply that 

a.W.- = e-/e-(r.r')a4rO*' (5.40) 

where exp(— /3/iQ,)(r, r') is the kernel of the operator exp(— /?/?,„). Inserting this relation in 
(5.39), using the cyclicity of the trace and the canonical commutation relations gives 

lnSA = ^?7«-^Tre-*-^] (5.41) 



dm 

where Tr is now the one-particle trace and it includes the sum over spin. In the case of zero 
external field this is the standard formula for the density of an ideal Bose or Fermi gas by 
summing the geometric series. 

Integrating (5.41) from — oo to yields 

Y^l—Tre-'^Pi^--^^^ (5.42) 

If we replace q in the exponent by f^--- ds and time-order the exponential then this formula 
and our derivation remain valid for a periodic time-dependent external field 0(r, s), and we 
obtain (4.73). 

The "magic formula" in the case of external field, no two-body interaction and one 
particle species follows by substituting in the Feynman-Kac representation 
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oo q_i „ 



<p{C) = / (j)(R{s),s)ds (5.43) 





where we use the notation described in subsection V.A.I. 

The "magic formula" (5.16) in the case of two-body interactions follows by integrating 

over imaginary timc-dcpcndcnt fields i0(r, s) with a Gaussian measure. According to our 
discussion of the quantum Sine-Gordon transform in section IV. F this reconstructs the two- 
body interaction with self-energies. The formulas (5.13,5.14) follow by absorbing the self- 
energies into the activities. 



B. The loop Mayer expansion 

1. The loop Ursell function 

Since the system of loops has a classical structure, it can be treated by the well known 
methods of classical statistical mechanics. In this section we apply the standard techniques 
of the virial expansion to the two-point Ursell function h{Ca, Cb) together with the partial 
resummations needed to deal with the long range of the Coulomb potential. 

The Ursell function is defined in terms of the two- loops truncated distribution pr^^a-i ^h) 

by 

PT{Ca, A) - p{Ca)p{C,)h{Ca, C,) (5.44) 

Although loops are extended objects, we still call them "points" and represent them graphi- 
cally by points in the diagrammatic language. The non-integrated loop variables Ca-, Cb are 
called root points. We recall that the Ursell function has the following simple loop-density 
expansion (Hansen & McDonald 1976) 

1 r ^ 

In (5.45), the sum runs over all unlabeled topologically different connected diagrams F with 
two root points and N internal points = 0,1,2,..., without articulations (an 

articulation point is such that, when removed, the diagram splits into two pieces, at least 
one of which is disconnected from any root point). To each bond (>Cj, Cj) in the diagram is 
associated the Mayer factor (called f-bond) 

/(£i, Cj) = exp{-(3ijV{Ci, Cj)) - 1, /3ij = /SeaiCaj (5.46) 

and [Yl f]r is the product of all f-bonds in the F-diagram; the symmetry factor Sr is the 
number of permutations of internal points £„ that leave this product invariant. 

It is clear from (5.30) and (5.31) that the loop potential V{Ci,Cj) is Coulombic at large 
distances, so that the bond (5.46) is not integrable. In the classical case it is known that 
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by resumming appropriate classes of graphs (the convolution chains leading to the Debye- 
Hiickel potential), one can introduce new bonds and diagrams (the prototype diagrams) 
that are free of divergences. One applies here the same procedure to the quantum gas in the 
loop formalism. The main guidance comes from the fact that the pure monopole part (5.31) 
would lead to Debye screening. Therefore the subsequent operations aim at disentangling 
the monopole from the multipolar effects, and then extracting as much classical screening as 
possible in order to reduce the range of the interaction. When all these steps are completed, 
no bond will decay slower than |r|~^, thus converting the gas of Coulombic loops into a 
certain gas of dipoles (together with higher order multipoles). 



2. Summing the Coulomb chains 

The basic f-bond (5.46) is spht into the sum of five contributions according to 

f /t ~I~ fee H~ fem ~l~ fmc ~l~ fmm i.^'^'^^ 

with the following analytical and graphical definitions (Figure 2) (abbreviating 

fT{iJ)^f{iJ) + Pi,jV{iJ) (5.48) 
fcc{i,j)^-Pi,jVcc{i,j) (5.49) 

fcm{ij)^-Pi,jVcm{ij) (5.50) 
fmcihj) = -Pi,jVmc{i,j) (5.51) 
fmm{ij)^-Pi,jVmm{ij) (5.52) 

Then, the Ursell function is given by the same formula as (5.45), summing now diagrams 
r with the same topological structure as the F-diagrams, but where the bonds can be any 
one of the quantities (5.48)-(5.52). Points without arrows in a F-diagram and to which are 
attached only two dashed lines are called Coulomb points (Figure 3). A Coulomb point 
involves only interactions with the total charge qCa of the loop (considered as a point like 
object) without multipolar effects. 

A Coulomb chain is a convolution chain where all intermediate points are Coulomb 
points. There are four types of Coulomb chains associated with the four bonds (5.49)-(5.52) 
(Coulomb chains without intermediate points are defined as the bonds themselves). As in 
the classical theory, we sum all Coulomb chains of a given type, and denote these sums 

-^CC) Pcm) ^rncj Pmrn 

(Figures 4 and 5). One has 

■^cci^aj ^b) — fcci^aj ^b) 

oo „ N 

+ ^ Yl ^^iP^^i)fce{^a, J^l) fecial, A) " " " fcci^^N, A) 

Ar=i i=i 

(5.53) 

Since by (5.31) fcc{j0.i,Cj) = —(3ijqiqjV{Iii — R^) does not depend on the configurations 
Xj, Xj of the loops, the summation on the internal variables can immediately be performed. 
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giving a factor / 9) for each intermediate point in (5.53). Because 

of translation invariance the loop-density p{C) = p{a, q, X) does not depend on the location 
R of the loop, therefore the multiple convolution (5.53) in the R variables sum up as in the 
classical chains to 

Fcc{a, b) = -Pa,bqaqb^Ra - Rb) (5.54) 
where $(R) is a Debye-like potential 

$(R) = (555) 
|R| 

with inverse screening length k given by 

- 47r/3 J] J] y" D{^)q'elp{a, q, X) (5.56) 



a q 



From the definitions (5.32)-(5.33) and (5.50)-(5.51), the sum of Coulomb chains correspond- 
ing to the charge-multipole bonds are obtained from (5.54) by application of the multipolar 
operator (5.29) on $(Ri - R2) (Figure 4) 

F^{1, 2) = -/?i,2gi / ds2M2{s2MRi - R2) (5.57) 
Jo 

F^,(l, 2) = -/5i,2g2 / ds,M^{s^)^R, - R2) (5.58) 

Obviously the resummed bonds (5.57) and (5.58) have a short range. Let us have a closer 
look at the multipole-multipole resummed bond (Figure 5) 

-^mm(l) 2) /mm(l) 2) + fmc ^ /'cm(l) 2) -|- fmc ^ fee ^ /'cm(l) 2) -|- • • • 

= -A,2V;„„^(l,2) (5.59) 

-/?1,2 / dSiM^isi) ciS2M2(s2)($(Rl-R2)-V^(Rl-R2)) 

^0 Jo 

= -/?V^(l,2)-A,2 / ds,M,{si) d52M2(s2)$(Rl-R2) 
7o Jo 

(5.60) 

The first term in (5.59) is the multipole-multipole quantum potential that has only 
"equal time" contributions, whereas the rest of the chain can be obtained by an application 
of multipolar operators at both ends (omitting here the term with no intermediate Coulomb 
points). The expression (5.60) results from (5.36), showing the asymptotic dipolar character 
of Fmm (see (5.37)). Notice that the sum of the resummed bonds equals 

F(l, 2) = FUh 2) + FUh 2) + F^c(l, 2) + F^^{1, 2) 

= -/3iy(l,2)-/?i,2$ci(l,2) (5.61) 

where $ci(l; 2) is a short-range interaction defined as in (5.27) 
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<^d{Cl,C2)= dsi ciS2$(Rl + Aa,Xi(si)-R2-A«2X2(s2)) (5.62) 

Jo Jo 

representing the screened classical interaction between loops. 

The screening length in (5.56) can be written in a more familiar form by noting the 
identity coming from (5.26) 

/OO p 

drp^^^\a,r,0) = J] / D(X)g(g - l)p(a, g, X) (5.63) 
q=2 

hence, using also (5.22) 

= 47T(3 el (pa + J drp^^''\a, r, 0)^ (5.64) 

An alternative expression for follows from the quantum charge sum rule (2.16) 

J drcT{r, 0) = ^ ea^Ca^ J dv (pr(Q;i, r, a2, 0) + Pai^oioa'^lr)) = (5.65) 



aia2 



which leads to the following relation, once expressed in terms of the loop distributions by 
(5.23) 

XI (p^ + J drp^^''\a,r,0) \ = - ^ e^^ea^ J drp^^^''\ai,r,a2,0) (5.66) 



011,012 



and thus 



-471(3 X eajCaj J dYp^^^''\ai, r, ^2, 0) (5.67) 

011,012 



In the integrals in (5.66) and (5.67), the contribution of coincident points is not included. 
On physically reasonable grounds, is indeed a positive quantity. One sees from (5.64) 
that reduces to the classical Debye value when exchange effects can be neglected, and 
it will coincide with the RPA value at high density (section VI.A.4), both quantities being 
positive. Moreover, the screening cloud / dvp^^^^\Qi, r, 0:2, 0) of charges of different 

species around the charge should have the sign opposite to Cq^, implying that the right 
hand side of (5.67) is positive. 

3. The prototype diagrams 

The F-diagrams can be divided into classes such that all diagrams in a class lead to the 
same so-called prototype Il-diagram when all the Coulomb points are integrated out. Points 
in a F-diagram that survive as points in a Il-diagram are either convolutions which do not 
involve Coulomb points, or points linked to three or more other points of F. There will be 
four different bonds (called F-bonds) in prototype diagrams: 
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(i) the screened charge-charge bond (5.54) and the two charge-multipole bonds (5.57)-(5.58). 

Since Conlomb points cannot occur in prototype graphs, the use of these bonds is 
subjected to the excluded convolution rule: in a Il-diagrani, convolutions Fcc^Fpc, i^mc* 
Fee, Fee Fem, -^mc * -^cm are forbidden. 

(ii) a bound Fi that incorporates all the bonds that are not taken into account in (i). 
Let us now construct the bond Fi. We distinguish two cases: 

(a) fx and F^^ occur as individual bonds in H-diagrams since they are never attached to 

Coulomb points in F-diagrams, 

(b) single bonds in a Il-diagrani attached to points that had more than two links in a 

F-diagram, as the result of the suppression of Coulomb points. 

The latter situation is obtained when two points in a F-diagram were linked 

(b.l) by any one of the f-bonds together with Coulomb chains in parallel (Figure 6) 

(b.2) by several Coulomb chains (Figure 7). 

The Coulomb chains occurring in (b.l) and (b.2) must have intermediate points since 
there is at most one f-bond between two points in a F-diagram. The sum of such Coulomb 
chains, Fc/i, is obtained by removing from F (5.61) the single bonds (5.49)-(5.52). Thus, with 
(5.30) 

Fehit,j) = F{t,j) + f3,,,V{t,j) = -(3W{t,j)+i3,,j{V{t,j) - $el(^,J)) (5-68) 
The contributions described in (a) and (b) sum up to 

Fi = fT + F„,rr, + f{exp{Feh) " 1) + exp(Fe,,) -l-Feh (5.69) 

The last two parts of (5.69) correspond to the contributions (b.l) and (b.2) (for the com- 
binatorial aspects, see (Cornu 1996a). Combining the definitions (5.46), (5.48), (5.61) and 
(5.68) in (5.69) leads to the final result (Figure 8) 

Fi{iJ) = eM-PW{iJ) - PijMhJ)) - 1 - Feeihj) - Fem{iJ) - F^eih j) (5.70) 
Eventually the loop-Ursell function is given by a sum of prototype graphs 11 

1 r ^ 

h{c,,c,) = J2^ U'^CnPiCnnU^w (5.71) 

n^^-J n=l 

with F-bonds Fee,F^e,Fem and F;. The Il-diagrams have the same structure as the F- 
diagrams, with the additional excluded convolution rule formulated in (i). 

All F-bonds are rapidly decreasing except Fi. Obviously the asymptotic behavior of Fi 
is dominated by that of W{Ci, C2) 

Fi{C^,C2) ~ exp(-/3VF(£i,£2)) - 1 ~ -/3W^(>Ci, A), |Ri - R2I ^ 00 (5.72) 
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which is of dipolc type as shown by (5.37). It is interesting to note that each Fi bond has 
the same form (up to short range contributions) as that found in the model of section IV. C 
when only two quantum mechanical charges are present. The long range part W depends 
solely on the intrinsic quantum fluctuations of the two particles (or the two loops) , and not 
on the temperature and density of the state. 

There arc of course several ways of reorganizing the diagrams with bare Coulomb bonds 
(5.46) into prototype graphs by partial resummations. Other resummations may be more 
adequate depending on the density and temperature regime and the physical quantity to be 
computed. The present reorganization appears to be optimal in the following sense. The 
remaining long range part has been isolated in the single bond F;. The decay cannot be 
improved any more by chain resummations. By (5.72) and (5.37), the longest range part 
of such chains would arise from convolutions of W^-bonds, but such convolutions keep their 
dipolar character (Cornu 19966); the same is true for classical dipoles, see for instance (H0ye 
& Stell 1974). Moreover, our prototype diagrams are the natural generalization of those used 
in the classical case. Indeed, if we set the de Broglie lengths equal to zero, the multipolar 
operator (5.29) vanishes as well as W , Fcm and Fmc- The effects of the quantum statistics 
disappear, and both $ (5.55) and <l>ci (5.62) reduce to the usual Debye potential Hence 
one finds that the only surviving bonds are 

limFee(l,2) = -/3i,2^,,(l,2) 
n—^0 

and from (5.70) 

lim Fi = e-^i-^*°(i'2) - 1 + /?i,2$i?(l, 2) 

They are precisely the two bonds occurring in the Meeron theory of the classical plasma 
(Meeron 1958, Meeron 1961). 



4. Integrability of prototype diagrams 

Let us present some crude arguments for the integrability of H-diagrams; they can be 
substantiated by the more thorough analysis of (Cornu 1996a), but a proof is still missing. 
Each internal point involves a density factor p{C) = p{a,q,X.) which depends only on 
the internal variables of the loop and not on its location in space because of translation 
invariance. Integration on an internal point jO. will be of the form 

/oo » 
dRY,Yl J ^(X)p(a,g,X)G(a,g,R,X;...) (5.73) 

where G{C ; . . .) denotes the /^-dependence of the bonds attached to this point. It can 
be argued that the loop-density p{jC) is bounded in the shape variables X and inherits 
a prefactor expdS/iaQ) from the activity of the loop (see (5.13), (5.15) and (5.74) below). 
Thus for a sufficiently negative chemical potential this factor should insure convergence for 
the q summation on the size of the loops. If these conditions are met, the integrability 
of Il-diagrams, when the distance between the loops becomes large, is determined by the 
the decay of the F-bonds. Because of the absence of articulation points, any cluster of 
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points is linked by at least two bonds to the remaining points of the diagram. Hence the 
product [H-^lnj once integrated on the internal variables, decays at least as |R|^^ in the 
loop distances, showing the integrability for the whole diagram. One should add that the 
splitting (5.47) of the original f-bond may introduce spurious non-integrable singularities at 
the origin in some diagrams. These singularities will cancel out when suitably collecting 
together such dangerous diagrams (the f-bond (5.46) gives a finite contribution at the origin 
because of the smoothening provided by the functional integration on the filaments). 

To summarize the situation, once the topological reduction to prototype diagrams has 
been performed, the two-particle correlation is reduced, in the sense of the formal diagram- 
matic cluster expansion (5.71), to a system that resembles a classical gas of multipoles 
described in the loop formalism by the four types of bonds Fee, Fmc, Fcm and F;. 

5. The density expansion 

An important remark has to be made when one comes down to computations, in par- 
ticular when one seeks to obtain explicit low density expansions of physical quantities. The 
diagrammatic representation (5.71) of the loop Ursell function is a straightforward expansion 
in the density of loops, but not in the particle densities pa themselves, for two reasons: 

• the resummed F-bonds depend on the particle densities through the density dependent 
screened potential $ (5.55), 

• the loop densities entering as weights of the internal points in H-diagrams are not 
simply proportional to the particle densities. 

The first point is not new and occurs of course already in the diagrammatic treatment of the 
classical gas. The second point is more delicate: the loop and particle densities are related 
by the formula (5.22); hence the particle-particle distributions (5.23) calculated from (5.71) 
have a very implicit dependence on the p^- In principle this dependence can be extracted 
as follows. According to the standard rules of the Mayer diagrammatics, the loop density is 
represented by the activity expansion 

N 

p{C) ^z{C)Y,— W dL^z{L^) [H !\g (5.74) 

G ^ n=l 

The G-diagrams are defined as the F-diagrams (with one root point) except that articulation 
points are allowed and weights are given by activities. The integrals diverge in (5.74) because 
of the long range of the /-bonds, but after a suitable topological reduction (resumming 
Coulomb chains), the representation of p{^C) becomes a series of prototype graphs, each of 
them being finite. The procedure is similar to that given above for the Ursell function and 
will not be described here in detail, see (Cornu 1996a). Once this is done, formula (5.22) 
yields in principle the particle densities as functions of the activities 

Pa = Pa(^i,---,^5), q; = 1,...,5 (5.75) 

Because of the overall neutrality Y^ci=\ ^aPa = the system (5.75) consists only of 5 — 1 
independent relations, and so does not determine the uniquely. A unique determination 
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of the Za is usually made by imposing the convenient additional relation (already discussed 
in section I.D) 

s 

^ e^z^ = (5.76) 

Thus by inverting the system (5.75)-(5.76) one can express first pi^C) and then the particle- 
particle correlations in terms of the p^- Explicit calculations in this formalism can of course 
only be performed in the low density regime. Applications to the determination of the 
asymptotic behavior of the correlations and to the virial equation of state will be presented 
in the next chapter. 
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VI. CORRELATIONS AND THE EQUATION OF STATE 



A. Asymptotic behavior of particle and charge correlations 

The first application of the general formalism concerns the determination of the large 
distance behavior of the correlations without recourse to h- expansions, including the quan- 
tum statistics and using the exact Coulomb potential on all scales. The calculation remains 
perturbative in the sense that we examine the large distance behavior of each H-diagram in 
the series (5.71) representing the Ursell function, but, as before, without control on the sum 
of these series. 

The semi-classical analysis of section IV. C has revealed that the two particle correlation, 
at order K^, decay as at large distance. We show that the same result holds for each of 
the H-diagrams contributing to the particle- particle correlation ^^'(aa, r^, a^, r^), i.e. 

Aab{l3, {pa}) 
pTKpCa, Ta, ftfe, r^j ~ — — 

\ra-nv 

where the coefficient Aa&(/9, {Pa}) is the formal sum of the asymptotic contributions of the 
Il-diagrams as |ra — r^l — cxd at densities pa and inverse temperature (3. It is then possible 
to calculate the lowest order contributions in the densities pa to Aab{l3-i {Pa}) leading to the 
exphcit formula (6.22) below. We will describe the general strategy and illustrate some key 
points by explicit calculations, details can be found in (Cornu 1996&). A short review can 
be found in (Alastuey & Cornu 1997&). 



1. Large distance behavior of the loop Ursell function 

The first task is to study the behavior of the Ursell function (5.71) in the space of loops 
as the loops are separated. This large distance behavior will be determined by that of the 
four bonds Fcc,Fcm,F^c and Fi, once integration on the internal loop variables has been 
performed in prototype graphs. Among these bonds, only Fi has the long range multipolar 
character (5.72) that may induce a slow decay of the graph as the two loops are separated. 
Therefore one has to keep track of the effects of F^-bonds (and more specifically of the basic 
quantum interaction W (5.35)) in the convolutions occurring in prototype diagrams. For 
this it is appropriate to recall some facts on the asymptotic form of convolutions in three 
dimensional space. 

The forthcoming analysis depends on the assumption (made throughout the whole section 
VIA) that the leading asymptotic term has a monotone decay, whose behavior is governed 
by that of the small wave numbers k ~ of the Fourier transform. Possible decay with 
oscillation, due to singularities at non-zero values of k, is not considered. Such oscillating 
decays are not expected to occur at positive temperature and low density. 

Consider a (non- oscillating) function g{r) which decays as r = |r| — >• oo as (up to a 
multiplicative constant) 

/ N (Inr)™ , , 

g{r) ~ , n,m integers, n > 3, m > (6.1) 
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This type of decay is reflected by the existence of a non-analytic part QsingO^) in its Fourier 
transform ^(k) as k ^ 0. One can split the low k-expansion of ^(k) into ^(k) = c/regO^) + 
^sms(k) where gregO^) = 00 + J^l^i a,^,k^ + Y.l,u=i^i^'^^fi^'^ + • • • has a Taylor expansion 
around k = and ^sms(k) behaves (up to multiplicative constants) as (Gel'fand & Shilov 
1964) 

~ / |k|"-3(ln|k|)™+i, nodd . _ 

The singular part of Fourier transforms of derivatives 9^-^ ■ • ■d^^g{Y) is 

^/ii ■ ■ ■ k^-Qsingi}^)- Since the Fourier transform of the convolution {gi -k g2){r) is the product 
^i(k)^2(k), one deduces from the above correspondence rules that for ni,n2 > 3 

(«/i*c/2)(r) ~ ^i(r) J dr'c/2(r') + c/2(r) J dr'gi{r'), r ^ oo (6.3) 

In particular, if 712 > rii > 3, the convolution decays as the slowest of the two functions 

(£'i*52)(r) ~ t/i(r) / cir'^2(r'), r ^ 00 (6.4) 



At the border ni = n2 = 3 of the integrability domain one has 

(^i^^2)(r)^^^^^:^ (6.5) 

One sees from (6.5) that convolutions of functions can generate in principle logarithmic 
factors in the power law decays. To simplify the discussion we abbreviate (6.1) by g{r) ^ r~", 
omitting from the notation the possible logarithmic dependence. Since Fi decays as r^^, it 
is clear from the above convolution rules that no H-diagram (or sub diagram) decays slower 
than r"^. The slowest possible decay comes then from diagrams where the two root points 
can be disconnected by cutting a single F^-bond: otherwise chains containing a F^-bond will 
be multiplied by another chain (or subdiagram) decaying at least as ~ so the decay of 
the whole diagram will be at least as ~ r~^. To put this idea at work, one splits the F;-bond 
into two pieces 

Fi = -PW + Fi, Fi = Fi + pW (6.6) 

and enlarge accordingly the set of Il-diagrams to Il-diagrams equipped with the new bonds 
-pW and Fi. One has from (5.72) 

F^l, 2) ~ eM-PW{l, 2)) - 1 + PW{1, 2) ~ ^{PW{1, 2))^ 
and {W{l,2)f ^ |Ri-R2|^oo (6.7) 

|±ti — ±121 

thus Fi is integrable and W is the slowest bond in the new set of It-diagrams. Introduce 
now the subclass of He-diagrams that remain connected when one removes a VF-bond, and 
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denote H{Ca, Cb) the value of the sum of these diagrams. With this definition, the loop 
truncated distribution function (5.44) can be expressed as an exact Dyson equation in terms 
of convolution chains involving H and W 

PT{Ca, A) = p{Ca)p{C,)H{Ca, C,)-pj dCrdC'^G{Ca, Cr)W{C,, C[)G{jC[, A) (6.8) 
+(3^ J dCidC[ j dC2dC'^G{Ca, C^)W{Cu JO,[)G{C[, A)^^!^, 4)^(4, A) + • • • (6.9) 

In (6.8)-(6.9) we have set 

G{C^, A) = p{CMJ^2)H{C,, A) + S{C,, C2)p{C,) (6.10) 

where the term S{jCi, jC2)piJ^i) takes into account the case when VF-bonds are in direct 
convolution with themselves. By construction, no diagram in H{Ca,^b) decays slower than 
~ |Ra — Rf,|^^. Indeed if a H^^-diagram contains no W^-bonds, the slowest possible decay is 
determined by that of Fi according to (6.7). If it contains a H^-bond, this W bond (or the 
chain where it belongs to) is necessarily multiplied by another chain that decays not slower 
than r~^, hence the overall decay of the diagram is at least r~^. In this way we have 
isolated the origin of decays that are slower than ^ in the chains of VT-bonds (6.8)-(6.9). 
Since (6.8) contains in particular the term p{Ca.)p{j^b)^{j^a, J^b) the loop correlation decays 
indeed as ~ r~^. 



2. Large distance behavior of the two-particle correlation 

In the loop formalism, the two-particle distribution consists of the two terms displayed 
in (5.23)- (5.25). The dominant part at large distances will come from p^""^^) involving parti- 
cles belonging to different loops. The contribution p^'^^^ of particles exchanged in the same 
loop will be discussed at the end of this subsection. To obtain p^"^^^^ we have to integrate 
Pri^^a-i ^h) = p{^a)p{^b)h{jCa, jCf,) ovBT the internal variables of the root loops £„, Cb accord- 
ing to (5.24): 

p("-)(a„,r„,a6,r6) = ^ j D{X^)D{X^)p{C,)p{Cb)h{Ca, C^) (6.11) 

The point to be made below is that, once these integrations have been performed, the asymp- 
totic contribution of the chains of H^-bonds (6.8)-(6.9) is strictly faster than pa |Ra — Rfe|~^. 
Hence the decay of is entirely determined by that of the ftc-diagrams constituting the 
function H, i.e. according to the analysis of the preceding subsection. 

Consider for simplicity the diagram consisting of a single VF-bond in (6.8): its contribu- 
tion to p("^^)(ra, r^) is 

oo oc p « 

-f^Y^Yl ^"^f / ^(^«) / D(^b)p{aa, Qa, ^a)p{ab, Qb, Xb)W {C^, A) (6.12) 

qa=l 9b=l 

The same mechanism that was put forward in the model of section IV. C to show that 
the term linear in W was in fact short-ranged works also in the present case. Since in a 
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homogeneous phase of the Coulomb gas the loop-density p{a, q, X) as well as the measure 
D(X) are invariant under global rotations of the shape X of the loop, the expressions 



/ 



D{X)pia,q,X) ^^^'^^^^^\ k = l,2,... (6.13) 



occurring in (6.12) according to the multipole expansion (5.36) of W have to vanish for k 
odd and are necessarily proportional to (V^)'^/^ for k even. Then the contribution (6.13) to 

(6.12) is strictly local since V|j^ (^]^) = 0, R 7^ 0: thus (6.12) decays faster than any inverse 

power. More generally, it can be shown that as a consequence of the rotational symmetry 
and of the harmonicity of the Coulomb potential, all chains (6.9) involving more than one 
VT-bond decay also strictly faster than ^ Ir^ — rtl"^. 

As an example, we give the argument for the term (6.9) involving a two VT-bond convo- 
lution and keeping only the slowest part of these Vl^-bonds, i.e. the dipole-dipole interaction 
(5.37). Once (5.37) has been inserted in (6.9) one obtains convolutions relative to the loop 
positional variables R of the form 



with 



{G^,, -k dy^d^^V -k G^^y^ -k dy^d^^V -k G^^){Yia - Rf,) (6.14) 

G,(R) = j D(Ka) J D(X)X,(s)G'(R,X„,X) (6.15) 
G,,{R) = J Z5(X') J L»(X)X;(s')X.(s)G(R,X',X) (6.16) 



where we write 



^(£1, £2) = ^(Ri - R2, Xi, X2) (6.17) 

because of translation invariance. In (6.14)-(6.17) we have omitted the loop indices a, q and 
the times s as well as the corresponding summations. The loop internal variables have been 
integrated and sums on repeated Cartesian indices are understood. To analyze the long 
distance behavior of (6.14) it is convenient to first consider its Fourier transform 

G., (k) '^G,,^. (k) ^(5,, (k) (6. 18) 

Since the measure -D(X) as well as G(R, Xi,X2) are invariant under spatial rotations , 
Gu{^) and Gfj_i,{^) transform as tensors: hence they are necessarily of the form 

G,{k)=k,a{m) 

G^.(k)=M(|k|) + fc^M|k|) (6.19) 

for some functions a, b and c of |k|. We know that G(R, Xi, X2) and hence ^^.(R), G'^j,(R) 
do not decay slower than ^ |R|~^- Taking the correspondence (6.1)-(6.2) into account, this 
implies that the small |k| expansion of Gi/(k) and G^,^(k) must be of the form 
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^.(k) = K(a + O (|k|2(ln |k|)"^i+i) 

G'^.(k) = 5^.(61 + 62|k|2 + 0(|k|=^(ln |k|)-^) 

+ A;^A;,(c + 0(|k|(ln|k|)-3) (6.20) 

for some mi, 7712,^3 > 0. Upon inserting (6.20) into (6.18) one sees that the Coulomb 
singularities |k|~^ are compensated and (6.18) behaves as C|kp + 0(|k|^(ln |k|)"*+^) as |k| — > 
0, hence (6.14) decays strictly faster than |R|~^ (possibly as ~ |R|~^). The same argument 
extends to all higher order multipoles and chains in (6.9) with arbitrary numbers of VF-bonds, 
showing that p*^"'^^^(aa, r^, a^, r^) decays at least as ~ Ir^ — rh\~^. 

We briefly comment about the exchange part p^'^^\aa,ra,rb) of the particle-particle dis- 
tribution. One sees from (5.25) that this exchange part involves the loop density p{jC) with 
the position of two particles fixed at and rf,. The Mayer expansion of the loop density 
(5.74) has an activity factor z{C) that has the bound (5.15). One can show that after the 
needed resummation of Coulomb chains, all prototype graphs representing p{jC) inherit the 
same activity factor z{£.). Furthermore prototype graphs are bounded functions of the loop 
variables ri, . . . , r^. Taking these facts into account in (5.26), the contribution of a prototype 
graph to p'^^^^Ua, r^, r^,) will be majorized by 

00 ^ 9 n , 

J2q / D{X)\z{£)\J2S{Xa.X, + v^-r,)<^^ 

q=2 k=2 



E e^^-- [dX,... dX, n ^'tw^ ""^''^ E ^(A..X. + r„ - r,) 

q=2 -I k=l ^^^^ k=2 



(2Sa„ + 1) 



^ J dkexp(-ik • (r„ - r^)) 



(2-)^^ "^^exp(/3(£|kp-/.J) 



(6.21) 



The integral in (6.21) is nothing else than the off-diagonal part of the one-body density 
matrix of a free Bose gas (the Bose distribution occurs because we have replaced the fermionic 
sign rja^ by 1 in the bound). At sufficiently low density {exp{f3paa < 1)) this function, 
the Fourier transform of a infinitely differentiable function, decays faster than any inverse 
power of |ra — r(,|. Under the assumption that the same property remains true for the sum 
of all prototype graphs, we conclude that the asymptotic behavior of the particle-particle 
distribution is dominated by its non-exchange part previously discussed. 

The diagrams contributing to the coefficient of the tail can also be analyzed in the 
low density regime (Cornu 1997a). At the lowest order in the densities, one finds that the 
tail has the same form (4.59) as in the semi-classical gas 

PT{aa,ra,ab,n) ~ ^ ^,[6 24o' / / drG''"-{r,ab) (6.22) 

where G^^^^{r, a) is defined as in (4.58), but calculated in the Debye-Hiickel approximation 



/ 



drG"^-^- r, a)= pa\ 5— > ^— 6.23) 
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This is to be expected, since at lowest order in the densities, the particles surrounding the 
two fixed charges at Ta and arc fully ionized and will behave as a classical plasma. It is 
interesting to specialize (6.22) to the electron-proton system. In this case one finds from 
(6.22) and (6.23) with a^electron or proton, e^i = —e^r = e, Pei = Ppr = P 

n2 h*l3*e'^ f I 1 \^ 
PT(«a,ra,a6,r6) ~ ZTJ^^- ( — + — J (6-24) 



ka — Tfol^ 240 4 Vmel m 



showing that in this system all tails have exactly the same decay rate at low density. 

The tail of the particle correlation in presence of a uniform magnetic field B has also been 
calculated (Cornu 19976, Cornu 1998a). Because of breaking of the full rotation invariance, 
it has the slower decay ~ '^^■^"^Js where P4(cos6') is a Legendre polynomial and 9 is the 
angle between B and — r;,. 



3. Correlations of the total charge density 

As in the semi-classical case, the total charge density enjoys better screening properties 
than individual particles. This is refiected by the fact that correlations involving the total 
charge density c(r) = ea/9(a, r) such as the charge-particle correlation 

< c{Va)p{ab, Tfe)) >T= ^ ea^pT{0ia, r„, Qife, Tfe) (6.25) 

Ola 

will have a faster decay than the particle-particle correlation. 

To illustrate the origin of this additional screening mechanism, it is convenient to con- 
struct a certain " dressing"' of the root point Ca in the Ursell function. For this we call a 
Coulomb root point in a H-diagram if Ca is linked to the rest of the diagram by exactly one 
Fee bond or one Fcm bond, and split the Ursell function accordingly 

h{Ca, A) = h%Ca^ £,) + /i"^(£„, A) (6.26) 

h^{Ca,Ch) {h'^'^{Ca, Ch)) is the sum of Il-diagrams where Ca is (is not) a Coulomb point. 
This definition implies immediately the following integral equation 

h'' = Fee + i^ccP * + + F^^p * h (6.27) 

and thus 

h = Fec + Fcm + {FceP + 1) ★ h""" + F^rnP * h (6.28) 

where here the symbol -k means both integration on the internal variables and convolution 
with respect to the positions of the loops; p = p{C) is the weight of intermediate points. 
Note that /i"'^ occurs in the second term of (6.27) because convolutions of F,.,. bonds are 
forbidden in Il-diagrams. As we know the two first terms of (6.28) are rapidly decreasing; 
the third term contributes to priCa, Cb) as 

p(Ca) f dC {F^{Ca, C)p{C) + 5{Ca, C)) h^%C, Cb)p(Cb) 
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= J dCS{Ca,C)p{C)h^%C,Ct)p{Ct) (6.29) 

with "dressing" factor of the point defined by 

S{Ca, C) = p{Ca)F,,{Ca, C) + C) (6.30) 

When summed on the internal variables of the loop and the charges, it is immediately 
seen, using (5.54) and (5.56) that (6.30) gives 

E E ^« / D{^a)S{Ca. C) = CaQa Ut) - k'^^^^) , r = R„ - R (6.31) 
with Fourier transform 



k2 \ |kl2 



Hence the convolution (6.29) will contribute to the Fourier transform of the charge-particle 
correlation as 

ll^l^ J Fourier transform of Il-diagrams constituting /i"'^ 1 , , 



+ |k|2 1 integrated on the internal variables of their root points 

Since the singular term of these diagrams is |k|'^(ln IkD™ (corresponding to a r^^ decay 
as seen in section VI. A. 2), the singular term in (6.33) is |k|^ln implying a ~ r^^ decay. 
The same conclusion can be drawn for the contribution of the last term of (6.28) on the 
grounds of rotational invariance and of the short range of Fcm- The final result is that the 
charge-particle has the following large distance behavior in the low density regime (Cornu 
1997a) 

(p(a„,r„)c(r,))T ~ ^4^^ / (^^ j dv\v\'K^ '' {v)^ (6.34) 

where the last factor is the expression (4.62) calculated in the Debye-Hiickel approximation 

i/*|rri.-(r)^4Et^ (6.35) 

a 

For the charge-charge correlation, both root points and can be dressed as before, 
leading to a f=:i r~^° decay with 

^■•.)«W)-^^(5/*W'^''''w)^ 
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The fact that the space integral of (6.31) vanishes (or equivalently (6.32) vanishes at k = 0) 
is the elementary form of the classical charge sum rule in the Debye-Hiickcl theory (see 
(2.9)). In the quantum gas, this vestige of classical screening can still operate to improve 
the decay of the charge-charge correlation from to r~^°, but not beyond. 

Finally the susceptibility (1.45) (i.e. the response to an external classical charge) decays 
as r~^. Its tail at low density is given by 

We emphasize that (6.22), (6.34), (6.36) and (6.37) are exact asymptotic results for the 
multicomponent quantum gas. 



4. A comparison with the standard many-body perturbation theory 

The standard finite temperature many-body perturbation theory is an expansion in 
powers of the coupling constant (the charge e^). Here we adopt the definitions and nor- 
mahzations of (Fetter & Walecka 1971) and (Cornu & Martin 1991)^2. We recall a few 
facts. The terms of the expansion are represented by Feynman graphs in which bonds 
(«!, ri, Si; «2i S2) are joined cither by free propagators ^a^^.ajG'a {^2 — ^i, Si — S2) or by 
interaction lines Uai,a2{^2 — -Si — S2) (Figure 9). At each vertex {a, r, s) are attached two 
propagator lines and one interaction line. 

A basic link between physical quantities and perturbation theory can be established 
through the imaginary time displaced charge correlation (1.43) by 

CT(r, s, 0) = - ^ CaC^Va-zir, s) (6.38) 

ay 

The total polarization PQ,^(r,r) consists of all connected Feynman graphs in which the 
points {a, r, s) and (7, 0, 0) are joined by internal lines. Va-yi^, t) is a time periodic function 
of period /3. Introducing its space and time Fourier transform V{\i,n), the structure factor 
(1.44) and the susceptibility (1.45) are given by {n indexes even Matsubara frequencies 
cun^ n even) 2^ 

00 

S{k)^- J2 5]eae7PaT(k,n) (6.39) 

n=— 00 aj 



^^There is a minus sign missing in front of the right hand side of (Cornu Sz Martin 1991, formula 

In general, n is odd for fcrmionic propagators and even for bosonic propagators, but the charge 
correlation (1.43) involves only quadratic expressions of the Fermi or Bose fields and only even 
frequencies occur in its time Fourier series. 
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X(k) = (3V{k) J2 eae^-Pa-yiK n = 0), V{k) = ^ (6.40) 



ay 



It is convenient to consider the proper polarization P*^(k, n), the set of all polarization 
parts that cannot be separated into two polarization parts by cutting a single interaction 
line. PQ-y(k, n) and V*^(k.,ri) are related by a Dyson equation 

V* (k n) 

An effective frequency dependent interaction can be defined by 

This definition is justified by the fact that its zero frequency component coincides with the 
static effective potential (1.30) (with e{k) given by (1.45)) 

U^"{k, n = 0) = £-^(k)V^(k) (6.43) 

One has also the relation 

V^^ik, n) - K^{k, n)+pYl ^*c.siK n)ese,U'"{k, n)V;^{k, n) (6.44) 

St] 

The well known random phase approximation (RPA) amounts to replacing 'P^^{r,r) by 
the lowest order proper polarization term, i.e. the product of two free propagators 

nPA,a7(r, S) = -VaSa-y{2Sa + l)G° (r, s)GI{-V, s) = 5^,V*nPA,a{r, s) (6.45) 

Hence from (6.41) and (6.42) (see Figure 10) 



and with (6.39) and (6.40), this defines also SRPA{k) and XRPAik)- 

Let us show that the RPA approximation fails to capture the algebraic tails in the 
correlations^^. The free propagator has the form 



The present considerations apply only if the temperature is different from zero. At zero tem- 
perature, the effective potential shows the long range Friedel oscillations cos(2/c j7r)/r^ due to the 
sharpness of the Fermi surface. 
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Glik, s) = exp (^-Ps - ) K(k) - e{s)] (6.48) 

where n^(k) = ^exp (^^^^^^ /S/^aj ± is the free Fermi or Bose distribution and 6{s) 

is the Heaviside function. It is clearly an infinitely differentiable function of k (if /i^ < for 
Bose particles). One can deduce from (6.45) that the same is true for ^^pa,q('')^) one 
finds the small k behaviour 

PS,,,Jk,„ = 0) = -^+O(|kP)^0 

■P'nPA.c(KnT6 0)=O(\kf), k^O (6.49) 

where p° = {2sa + 1) / c^kn^ (k) is the density of the free gas. With (6.46) this implies that 
the 'PRP^,Q(k, n) are infinitely differentiable at k = and by (6.39) and (6.40) SjiPA(k) 
and XRPA(k) inherit the same property. Thus, in space, these functions decay faster than 
any inverse power of the distance. From (6.46) and (6.49) one sees that XRPA{k.) can be 
approximated for small k by the Debye-like form 

1^1 '^RPA 



with 



'=W = 4./JE4j|5 (6.61) 

SO defining the RPA screening length K^p^. It is instructive to compare it with the loop 
screening length (5.64). For this, we introduce the general relation 

d /■ , . 

-KTB — tP« = Pa+ dvpT{a, r, a, 0) 

^Pa + J drp^^''\a, r, 0) + J drpj'^''^ (a, r, a, 0) (6.52) 

in (5.64), so that can also be written as 



4^/5 (E - E / ^rpr^^(«, r, a, 0) j (6.53) 



In the integrals in (6.52) and (6.53), there arc no contributions of coincident points. The loop 
screening length interpolates between the classical Debye length at low density (as seen from 
(5.64) when exchange effects can be neglected) and the RPA length at high density. Indeed, 
in the high density regime, the kinetic energy dominates the Coulomb interaction energy: 
the gas becomes free so the non-exchange part of the correlation can now be neglected in 
(6.53) and tends to the RPA value. 

As far as the effective potential (6.42) is concerned, one finds in the RPA approximation 
that its zero frequency component U'^/^{k,n = 0) is also Debye-like for small k with a 
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screening length slightly different from K^p^ (for explicit formulae see (Cornu & Martin 
1991, Cornu 1996&)). However its non-zero frequency components f/^'^^(k, n 7^ 0) are not 
shielded and keep a pure Coulombic singularity |k|~^, but within the RPA theory these 
singularities do not induce any algebraic tails in the particle and charge correlations. For 
instance, if 7'*^(k, n) is approximated by Vj^pA^^^(k.,ri) in (6.44), the vanishing of the non- 
zero frequency components as |k| kills the |k|~^ singularity coming from t/^'p^(k, n ^ 0). 

One may however wonder how these tails emerge in the framework of the many-body 
perturbation theory. The role played by W^/^ in Feynman graphs is analogous to that of 
the resummed bond Fi (5.70) in loop prototype diagrams. F; has a long range part W 
(see (6.6)), but this long range part generates tails in the correlations only if it occurs at 
least quadratically in prototype graphs (see section VI. A. 2). In the same way, U^/^ must 
enter quadratically in Feynman graphs to generate tails in the correlations. This leads to 
including as a first correction to the RPA approximation the proper polarization insertion 
Vi made of two propagator loops linked by two RPA lines of Figure 11. The effect of this 
correction has been studied in the one component system (Cornu & Martin 1991). The 
Coulombic singularities of U^Jj^{\<., n 0) come now non-trivially into the game and one 
finds that Pi(k, n — 0) and 'P*(k, n 7^ 0) have respectively a |k|^ and a |k|^ singular term in 
their small k expansion. One can then infer through the Dyson equation (6.41) that these 
singular terms induce r~^° and tails in S{r) and x(r) respectively, in accordance with 
the findings of section VI. A. 3. The graphs of Figure 11 have been considered previously 
in (Maags & Ashcroft 1987, Langreth & Vosko 1987) in connection with calculating the 
cohesive energy of metals and inter-atomic potentials at zero temperature. 

It is not possible to establish in general strict correspondence rules between the many- 
body perturbation algorithm and that of the loop formalism. Since the reference system of 
the many-body perturbation scheme is the noninteracting gas, any approximation there gives 
expressions in terms of free quantities. In the loop expansion,weights and bonds in Mayer 
diagrams contain already the coupling constant and the statistics in a non-perturbative 
way, and they involve the exact particle densities. This difference is clearly seen in the two 
screening lengths K^p^ and that depend respectively on the free densities and exact 
densities Pa- One can however establish a precise relation between the loop chain potential 
and the RPA effective potential in the Boltzmann limit when statistics are neglected (see 
(Cornu & Martin 1991, section IV)). 

The two formalisms play complementary roles. The many-body formalism is suitable 
for performing low temperature and ground state calculations. By analytic continuation 
in time it can also be used to compute real time displaced correlations. Information on 
ground state quantities in the loop formalism would necessitate applying the Feynman- 
Kac formula backwards at some point before taking the /3 — >^ 00 limit. However the loop 
formalism provides a better insight into Coulombic correlations and screening. For instance 
the reduction to a gas of multipoles is not visible in the many-body framework where the 
effective potential U^Jj^ keeps an unscreened bare Coulomb part. Low density behaviors are 
much more conveniently obtained from the loop expansion: it would not be straightforward 
to recover the exact results of section VI. A from an analysis of Feynman graphs. The next 
section offers another example of the usefulness of the loop formalism. 
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B. The virial equation of state 

The derivation of the equation of state of a quantum mechanical plasma is an old problem 
of practical and theoretical interest. In particular, a precise knowledge of this equation plays 
an important role in astrophysical conditions, and it is not possible to give credit here to 
the large literature devoted to this question. We quote the thorough study of the quantum 
Coulomb gas undertaken by Ebeling and collaborators by the usual methods of the many- 
body problem. This work is reported in (Kraeft, Kremp, Ebeling & Ropke 1986) and some 
comments will be given at the end of this section. In the sequel, we put forward that the loop 
representation is also a convenient tool for calculating the equation of state at low density 
and fixed temperature. In particular the formalism enables simultaneously keeping track, in 
a coherent and systematic way, of the quantum effects and of the long range of the Coulomb 
potential, a notoriously difficult problem. A detailed presentation of all calculations can be 
found in (Alastuey & Perez 1992, Alastuey, Cornu & Perez 1994, Alastuey, Cornu & Perez 
1995, Alastuey & Perez 1996) and a review on this subject is presented in (Alastuey 1994). 
Here we outline the general computational scheme and illustrate some salient points. 

To evaluate the density expansion of the thermodynamic functions, it is appropriate to 
start from a standard identity relating the free energy f{(3, {pa}) per unit volume to the 
two-particle correlation with a variable coupling parameter 



Here pT,(/(Q;a, r, ccb, 0) is the truncated two-particle correlation (without the contribution of 
coincident points) for a system of charges interacting by the Coulomb potential g/\r\ with 
a dimensionless coupling parameter g, < g < 1. For each g, prg(aai Tj ctfci 0) must be 
evaluated at the same values of the densities and temperature that occur as arguments of 
f{P,{pa})- The particle correlation PT,g(cto, i"a, ctft, r;,) is in turn related to the loop corre- 
lations through (5.23), (5.24) and (5.25) and the loop correlations have the diagrammatic 
expansion discussed in section V.B. Then one has to organize the contributions of this 
expansion in terms that are of increasing order in the particle densities and insert them 
in (6.54). This will yield the desired expansion for the free energy. Finally the pressure 
P{P, {pa}) follows from the usual formula 



1. The Maxwell-Boltzmann pressure 

At this point one can make the following comment about the treatment of exchange 
effects. It is clear from (5.13) that g-particle exchanges are of decreasing order exp{j3pa(l) 
as q grows {pa negative); hence they can also be treated perturbatively at low density. This 
is the view point adopted in (Alastuey et al. 1995) and that is reported here. One considers 




(6.54) 




(6.55) 



The calculation has been carried out up to order p^/^ . 
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first the Maxwcll-Boltzmann (MB) CoTilomb gas, and one includes the exchange effects in a 
second perturbative stage. The MB system is simpler since loops reduce to filaments. Also 
the exchange part (5.25) of the correlation can be ignored in the first stage. Dropping here 
the Fermi statistics should not alarm the reader. All terms of the low density expansion 
remain well defined, since these terms involve only finitely many body Hamiltonians, and 
thus the stability of matter plays no role at any finite order. As already noted, the local 
singularity of the Coulomb potential is smoothed by the usual laws of quantum mechanics 
(the uncertainty principle) as in atomic physics. 

In the works (Alastuey et al. 1994, Alastuey et al. 1995), prototype diagrams are con- 
structed with a slightly different definition of resummed bonds than in section V.B.3. One 
has here, keeping similar notations as in section V.B.2 and V.B.3 

the charge-charge bond Fee, the charge-dipole bonds Fed, Fdc 

the dipole-dipole bond Fdd, the residual bond Fr (6.56) 

The bond Fee is identical to (5.54) (with qa = qb — '^); Fed and Fde are defined as in (5.57) 
and (5.58) but retaining only the dipole part 

Fed{l,2) = -^e,e, f ds,{\,Us,) ■ V^)^^^^^^^^^ (6.57) 

The bond Fdd is equal to the dipole part of the second term in (5.60), i.e 

Fdd(l,2) = -/?eie2 



while Fr incorporates all the other contributions (higher order multipoles and quantum 
effects). Note that in this decomposition we have the two types of long range r"^-bonds Fdd 
and Fr (in contrast to the prototype graphs of subsection V.B.3 where only Fi is dipolar). 

As already explained in subsection V.B.5, there are two sources for the density depen- 
dence: the weight p{a,^), that has itself an expansion in the particle densities Pa, and the 
screening length n occurring in the resummed bonds. One performs first the activity expan- 
sion of Pa with the help of appropriate prototype graphs (here the bonds are as in (6.56) 
but weights are activities), and one eliminates activities in favor of densities with the help 
of the relation Pa — J D{^)p{a, ^). The analysis shows that p{a, ^) is given as double power 
series in p^^^ and In p with first terms 

p(«, ^)^Pa + Yl PcPi fdrf D{^,) (6.59) 
X exp (^-PeaCy dsV{\r + Xy^^{s) - Xa^{s)\)^ 

-J D{^^)exp(^-peaey dsV{\r + Xy^^{s) - Xa^^{s)\) 
+ Oip'/') 
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These series are introduced in the statistical weights of the expansion of the Ursell function. 
FinaUy one obtains also the free energy (6.54) and the pressure (6.55) as double power series 
in p^l'^ and In p. Half powers of p come from the Debye length ~ in the screened 
potential, while In p terms occur because the r~^-bonds and Fj. are at the border line of 
integr ability in three dimensions. 

The density dependence contained in the Debye length can be extracted by a scaling 
analysis. Introducing the dimensionless variable a; = k|Ri — R2I one sees that F^c scales as 
Fed {Fdc) as and F^^ as k'' . More generally a higher order multipole term of the form 

|rLi — rL2| 

scales as «;'=i+'=2+i_ in the quantum effects (associated with the de Broglie length Aq,) 
and classical screening (associated with the Debye length k"^) can be disentangled by a 
decomposition of the form 

Fr = F,(l + kE^ + K^E2 + ■■■) + ^^G2 + K^Ga + • • • (6.60) 
where Fq is the quantum (truncated) bond 



F,(l, 2) = exp (^-^2 ^ dsV{\v + - Aa.^2(«)l) 

- 1 + V ~ ^ + '^^^ ^5[Aa,li(«) • V - K^iM ■ V] (^i 



(6.61) 



and the and Gi are scaled functions involving contributions of higher order multipoles. 
We distinguish three different types of effects in this game: 

• classical screening manifested by the occurrence of the Debye length in the scaled 
bonds 

• quantum diffraction effects due to the coupling of quantum fluctuations ^ to the clas- 
sical Debye potential as in Fee, Fed and F^a 

• quantum bound and scattering states appearing through F^. 

As an illustration the table in Figure 12 gives a (non exhaustive) list of examples of graphs 
contributing to different effects in the MB free energy. The calculations are tractable up 
to order p^/^. There are two simplifications up to this order: internal weights p(tt,^) can 
always be replaced by their lowest order p^, and only low order multipoles contribute in view 
of the above mentioned scaling arguments. The MB pressure is given by the terms (6.66) to 
(6.69) in the final result stated in the next subsection. 
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2. The exchange contributions and the final result 



To include the exchange contributions perturbatively, one has to come back to the ex- 
pression (1.40) of the grand-partition function with statistics and collect together all the 
terms involving the exchange of exactly n particles 



= sr + E-f ^ = -a' 1 + E 1^ (6-62) 




where S^f^ is the Maxwell-Boltzmann partition function and the contribution with Na ex- 
changed particles has necessarily an activity factor exp{f3iJ,aNa)- The expansion (6.62) 
generates in turn an expansion of the grand-canonical pressure 

oo 

P = /3-^ lim — lnSA = P^^ + y"^^"^ (6.63) 

with P^^ the Maxwell-Boltzmann pressure and P^"^ the contribution of n exchanged 
particles. The matrix element in (1.40) corresponding to a two particle exchange, 
say (r2, ri, r3, . . . , rjvl exp(— /Si^A.Ar) ki, r2, r3, . . . , rjv), will involve N ~2 closed filaments 
rj^j(s), i = 3, . . . , N together with two open filaments ri 2(5) and r2,i(s). Keeping these two 
open filaments fixed, but performing the configuration integrals and grand-canonical sums on 
all the remaining closed filaments leads to the consideration of the inhomogeneous Maxwell- 
Boltzmann partition function 2^^^'"'^""^^ It is the grand-partition function of a system of 
(closed) filaments in presence of the potential due to these two open filaments. Similarly, 
higher order exchange terms can be viewed as partition functions of closed filaments in pres- 
ence of an external potential created by the open ones. The density of filaments pC'"/*"™) (JF) 
in the inhomogeneous system can be expanded in terms of the correlations p{J-), Pt{^i, ^2) 
of the homogeneous gas by standard methods. Finally, one can apply to the latter correla- 
tions all the diagrammatic techniques previously described to obtain the density expansion 
of the MB gas. There is however one caveat. The particle densities are given in principle 
by the usual relation 

Pa = ^^({/^7}) (6-64) 

where P is the total grand-canonical pressure (6.63) including the statistics. But 
has been computed directly as function of densities that we call now within the MB 
approximation. Then the can be eliminated in favor of the true densities pa by noting 
that both P and P^^^ are evaluated at the same values of the chemical potentials Pa- Hence 
we have the thermodynamic relation 

= aliB/^H^r}) (6-65) 

with /^''^ the free energy of the MB gas. Using (6.64) together with (6.65) enables us to 
express the Pa and p^™ in terms of the desired particle densities p^. The final result is 
(Alastuey & Perez 1996) 



VI-15 



^ 247r 

a 



(6.66) 



(ln2-l) /I 3, 1, , ^ 
^ ^ + f - - -ln2+ -ln3 ) /3«:e„e^ 



a'yS a'ySf] 



a 



(6.67) 
(6.68) 



V2 

>i'rQ{xa^) + — /3^e^e5ln(fi;A„^) 



+ 



(6.69) 
(6.70) 



with 



+ 0{pHnp) 



07 



m„m^/(m„ + m^), = 15.205 ± .001, C2 = -14.733 ± .001 



In (6.69), Xa-y = —V^Pcae^/ Xa-y and Q(:x„^) is the quantum second order virial coefficient 
first introduced in (Ebeling 1967, Ebehng 1968) 

QM = -tJ-^ lim I / dr [(27rA^^)3/2(r| exp(-/3i/«^)|r) 



1 + 



r 2r^ 



while E{xaa) is the exchange integral 



£;(xaa) =2V7r y (ir(-r| exp(-/3i7c,a)|r) 



(6.71) 



(6.72) 



In (6.71) and (6.72) Ha'y is the relative one-body Coulomb Hamiltonian for a particle of 
mass vria'y submitted to the potential and C — .577216 ... is the Euler-Mascheroni's 
constant. 

The formula (6.66)-(6.70) incorporates and generalizes several earlier results. The per- 
fect gas contribution YliaPo^ (^^^ dominant one) has first been established rigorously in 
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(Lebowitz & Pena 1973), completed by (Hughes 1985) The next term is the lowest order 
correction to the free gas: one recognizes the familiar Debye-Hiickel contribution (2.4). The 
quantum diffraction effects appear in the term (6.68) at the order p^/^. The term (6.69) 
incorporates the total contribution of bound and scattering states of the two-body Coulomb 
Hamiltonian^^. The regularization of the integral in (6.71), which insures that Q{xa'y) is 
finite, is not arbitrary but comes from the structure of the bond Fq (6.61). One has here 
a contribution, the proper second order virial term, and a p^/^ multiplicative correction 
that arises as a many-body mean field effect. Finally, the exchange term (6.70) has a similar 
structure. 

The terms up to order p^/^, including the exact -virial contribution (6.69) are found 
in the work of Ebeling (Ebeling 1967, Ebeling 1968) and related references (Kraeft et al. 
1986). The full formula (6.66-6.70), up to the diffraction term (6.68) appears in (Kraeft et 
al. 1986). These works use the method of the effective potential (or method of Morita). In 
this method, one associates to the quantum Coulomb gas an equivalent classical system of 
point particles interacting with an effective many-body interaction. In the case of the MB 
gas, this effective many-body potential \&iv is defined by interpreting the diagonal element 
({ri}| exp(— /3i7jv)|{rj}) as a classical Boltzmann factor 

N . .3/2 

({r,}|exp(-/3i/^r)|M) = n ( ^) exp(-/?*jv(ri, . . . , r^v)) (6.73) 

1=1 ^ ' 

One can decompose n into a sum of two-body, three-body and higher order interactions, 
the two-body potential being 



*2(ri,r2) = -/?-Mn 



2 

,2 \3/2 



n {^^>^lf (ri,r2| exp(-/3i/2)|ri,r2) 



.1=1 



(6.74) 



Calculations up to order p^^^ have been performed by keeping only the two-body potential 
(6.74), see references in (Kraeft et al. 1986). The long range Coulombic part of (6.74) can 
again be eliminated by summing the convolution chains and one can use the Abe-Meeron 
diagrammatics. This procedure gives the above formula for the pressure except (6.68). 
It turns out that one has to take the three-body interaction into account to recover the 
diffraction term. 

One can specialize the equation of state to the one component plasma and take the clas- 
sical limit^''': these classical terms coincide with those calculated in (Cohen & Murphy 1969). 



Hughes studies the continuity of the free energy at vanishing density in relation to the techniques 
of Fefferman discussed in the next chapter. 
^^The ln(«;Ao!7) t^rm can be recombined with Q{xa'y)- The splitting is introduced to have Q(xa-y) 
defined as in (Kraeft et al. 1986). Note also that our de Broglie length differs from that in (Kraeft 

et al. 1986) by a factor ^2. 
^^This can be achieved by considering an asymmetrical two component system and letting m2 

oo, 62 ^ 0, p2 ^ oo with ^2^2 finite and p2e2 + Pi&i = 0. The one component plasma remains well 
behaved in the classical limit because of the spreading of the positive charge. 
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One retrieves also the fi-corrcction njS'^e'^p'^h'^/G'm to the classical equation of state (Pollock 
& Hansen 1973). This correction has to be of order p^: indeed one checks that the h^p^^"^ 
contributions coming from (6.68) and (6.69) compensate exactly. Finally high temperature 
series are recovered by expanding Q{x) and E[x) in powers of (such expansions are 
derived in (Kraeft et al. 1986), and these series agree with the expressions obtained by the 
methods of (DeWitt 1962, DeWitt 1966, DeWitt, Schlanges, Sakakura k Kraeft 1995). Ex- 
plicit formulae for various thermodynamic functions as well as their specialization to these 
particular cases can be found in (Alastuey & Perez 1996). 

The generalization of the low density equation of state to a uniformly magnetized plasma 
is given in (Cornu 1997&, Cornu 1998&). 
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VII. THE ATOMIC AND MOLECULAR LIMIT 



A. The electron-proton gas in the Saha regime 

The diagrammatic technique developed in chapter V has provided a systematic low den- 
sity expansion of the equation of state at fixed temperature. In particular, the lowest order 
term PP = Yl'^=iPa (6.66) represents a mixture of perfect Maxwell-Boltzmann gases 
constituted by the S species of fully dissociated charges: it is a free plasma state. Possible 
two-body bound states appear only as corrections in the second order virial coefficient. In 
this chapter, we will treat another regime where the basic constituents are now chemically 
bound entities behaving as ideal substances. This regime ( called the Saha regime) is char- 
acterized by a joint limit where both the density and the temperature go to zero in a coupled 
way. The (exponentially fast) rate at which the density is reduced as T — > determines 
a certain energy-entropy balance, selecting in turn the formation of some specific chemical 
species. Indeed low temperature favors binding over ionization, whereas low density, by 
increasing the available phase space per particle, favors dissociation. This limit is called the 
molecular (or atomic) limit. It was first formulated in precise terms by (Fefferman 1985) 
and further studied by (Conlon, Lieb & Yau 1989, Graf & Schenker 1995 &). Equilibrium 
ionization phases in this context are discussed in (Macris & Martin 1990). The Saha regime 
is also described in the usual language of many-body perturbation theory in (Ebeling, Kraeft 
& Kremp 1976, Kraeft et al. 1986). 

In order to understand the issues involved in a simple setting, we consider the electron- 
proton (e-p) system and the possible formation of hydrogen atoms (see also (Fefferman 
1986, Martin 1993) for general discussion and background). Let us adopt for a moment 
the thermodynamic view point that a chemist would take in presence of three preformed 
species, the electrons (e), the protons (p) and the hydrogen atoms (a) in their ground state. 
As a first approximation he takes into account the binding energy of each chemical species, 
but otherwise treats them as perfect gases of point particles. In this approximation, the 
grand-canonical densities are 

(771 ■ \ ^^'^ 
^-^J e^p{-f3{Ej - Hj)), j ^e,p,a (7.1) 

where fij are the respective chemical potentials, rria = rrie + mp. Eg = Ep = and Ea < 
is the ground state energy of the hydrogen atom. The law of chemical equilibrium for the 
dissociation reaction e+p a requires 

IJ-a^ l^e + l^p (7.2) 

and one also must have the neutrality pe — Pp- Taking (7.2) into account and introducing 
the combinations 



^^This density pj = pidcalj SiS well as the densities (7.28) are those of ideal gases to which the 
interacting system will be eventually compared. For brevity, we drop the mention ideal here and 
in (7.28). 
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it is easily seen that the neutrahty condition imposes the choice 



and hence (7.1) becomes 



Pe = Pp= [ o /3fc2 exp(/5/i), p„ = - exp{-l3{Ea - 2/i)) (7.5) 



Finally, since the species are treated as non-interacting, the pressure obeys the law of perfect 
gases 

/3P = Pe + Pp + Pa = P (7-6) 

In (7.6), the pressure is expressed in terms of the total number density p of protons and 
electrons p = Pe + Pp + '^Pa and we have introduced the Saha coefficient for the degree of 
ionization 

7 = — < 7 < 1 (7.7) 

Pe + Pa 

A low density-low temperature regime can clearly be obtained by choosing p negative and 
letting P ^ oo. Let us examine various cases noting from (7.5) that pa/ Pe ~ exp(— /3(£'a— /x). 

(i) Fix p < Ea- The atomic density is exponentially small compared to that of the electrons 

and protons; there is full dissociation (7-^1) and the equation of state reduces to 
pp = p as in the low density expansion (6.66). 

(ii) Fix p > Ea (but p < Ea/2). One obtains the opposite situation where the pure atomic 

phase dominates; now 7^0 and pP = |p as /3 — > 00. 

(iii) At the borderline p — Ea, Pe and pa are of the same order. More precisely, replace p 

by 

p{P)^Ea + Xp-' + o{p-^) (7.8) 

for some A, —00 < A < 00, i.e. approach the point {Ea,0) in the {p,T) plane along a 
direction having a finite slope A. Then one obtains ionization equilibrium phases that 
interpolate between the fully ionized and the atomic ones. Their degree of ionization 
found from (7.5), (7.7) and (7.8) is 



VII-2 



and they obey the equation of state (7.6). Written in terms of the pressure, the Saha 
coefficient reads after a short computation 



rriemp 



3/2 \ ~V2 



^^H^^p,MSK±^, ,,p(-/Ji;.)| (7.10) 



This discussion makes clear that to an increase of /j, corresponds an increase of the densities 
that favors the electron-proton binding when /i crosses the value Ea- How can we justify 
this picture from statistical mechanics? Consider now the exact infinite volume pressure 
P(/3, fie, fip) of the e-p system in the grand-canonical ensemble, and recall from section I.D 
that it depends only on /j, — ^fs±fe 

Pif3,fle,l^p)^PW,l^) (7.11) 

One may expect that the equation of state associated with P{(3, fi) becomes close to (7.6) 
as /3 — > cxD for appropriate values of /x as in cases (i)-(iii). The theorem on the atomic and 
molecular limit will gives a precise formulation of this assertion. 

In the elementary considerations just presented, we have not taken into account, among 
other things, that electrons and protons can form other complexes as, for instance, hydrogen 
molecules. Let us examine this issue from the view point of statistical mechanics. If typical 
configurations of charges are indeed those of a dilute gas of hydrogen atoms (case (ii)), these 
configurations will be mainly formed of e-p pairs, the extension of a pair being of the order 
of the Bohr radius gb, whereas two different pairs are at distance pa {as ^ pa ^^^)- If we 
choose at random in this configuration a region D of linear extension R with 

as < i? < p-^'^ (7.12) 

we will observe that this region is empty most of the time, but if it contains something it is 
exactly one hydrogen atom, except on rare occasions. In the grand-canonical formalism, it 
means that one must be able to find chemical potentials pe, l^p such that for /5 3> 1 and all 
(7V„Ar^)^(0,0),(l,l) 

Poo{(3,pe,Pp,D) :$>pn{f3,Pe,Pp,D) :$>pN,Np{(3,Pe,Pp,D) (7.13) 

where 

PNeNj,W, Pe, fJ'p, -D) = J- exp{p{peNe + PpNp))TYexp{-(3HD,NeNj,) (7.14) 

is the probability to find exactly Ng electrons and Np protons in D. Thus the condition (7.13) 
means indeed that if D is not void, it is more probable to find an e-p pair than anything 
else. The dominant term in Tr exp{— (3 Ho, NeNp) as /3 — > oo will be exp{— (3ED,NeNp) where 
ED,NeNp is the ground state energy of HD,NeNp- Hence to obtain the inequalities (7.13), one 
must be able to find values of pe and pp such that 

Ed,n,Np - (PeNe + PpNp) > Ed,!! - {Pe + Pp) > ^ for all (A^e, Np) ^ (0, 0), (1, 1) (7.15) 
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Since pa tends to zero exponentially fast as /9 ^ oo we can also let the region D grow with 
j3 (but maintaining the inequalities (7.12)) so that £^£),iVeArp will differ from the bottom of 
the spectrum of H^^n^ in infinite space 

En.Np = inf spectrum(i?;VeJVp) (7.16) 

by a vanishingly small error. Moreover, introducing //,, i> as in (7.3) 

HeNe + iipNp = ii(Ne + Np) + u{Ne - Np) (7.17) 

we remark that the part n of the chemical potential controls the total number density and 
u the total charge density. We can take advantage of the fact that the infinite volume limit 
of the pressure will not depend on z/ to make a convenient choice. As it was done in (7.4), 
wc fix u to ensure neutrality in the ideal gas when the lowest energy state corresponds to 
charged complexes. Then u will be of order 0(/3~^) as in the dissociated e-p system and can 
be dropped from (7.17) as P gets large. Therefore, in view of these remarks (that will be 
made mathematically more precise in section VII. C) the relevant inequalities to be satisfied 
to produce the situation (7.13) as /3 — > oo are (with En — Eg) 

En^n, - l^{Ne + Np)> Ea-2ii>0 for all {N^, Np) (0, 0), (1, 1) (7.18) 

for some values of 

Let us examine some implications of (7.18). If (A^e, Np) — (1, 0) (single electron) Eio — 
implies in (7.18) 

Ea<ii<0 (7.19) 

If {N,,Np) = (2,2), (7.18) gives 

PL<^{E22-Ea) (7.20) 

where £'22 is the ground state energy of the hydrogen molecule. One can find fi satisfying 
(7.19) and (7.20) only if £„ < 1(^22 - K) or equivalently if \E22 - 2Ea\ < \Ea\. The latter 
inequality means that the binding energy gained by the formation of an hydrogen molecule 
must be less than the binding energy of the atom itself, a well known fact. It allows the 
entropy to win over molecular binding at low density. However the numerical values of 
(£■22 — Ea)/2 ~ — 9ev and of /i ~ — 8ev, k^T ~ 0.026ev under prevailing Earth conditions 
give an exceedingly small value of the ratio pa/ pu^ ~ 6xp(/9(£22 "-E'a ~2/i)). Only molecular 
hydrogen appears on Earth, while the formation of the atomic gas discussed here can only 
occur in very dilute extraterrestrial medium. These considerations on molecular hydrogen 
do not suffice to guarantee the validity of the full set of inequalities (7.18): this necessitates 
a simultaneous examination of all ground state energies En^^n^, a highly nontrivial problem. 
We shall discuss it further at the end of the next section. A geometric ilhistration of the 
inequalities (7.18) is given in Figure 13. Notice that the binding energies of the ion H~ 
and the ionized molecule H2 are too weak: they do not meet the inequalities analogous to 
(7.18), appropriate to these cases, for any value of /i. 
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Coming back to the e-p pressure, we may conjecture that if fi is increased just beyond 
|(i?22 — Ea), the chances for binding also increase, and the pressure P{P,^) will become 
close to that a free gas of hydrogen molecules as /3 — > oo. More generally, if one fixes ^ < 0, 
one expects that P{P,iJ,) approaches the pressure of free gases of some complexes that are 
selected at this particular value of the chemical potential. 

B. The main theorem 

With these preliminaries in mind, we are ready to formulate the main theorem about the 
Saha regime for electrons and an arbitrary number of nuclei with masses ma, charges Cq,, 
particle numbers N^, and chemical potentials /la, — 1, . . . ,S {a — 1 refers to electrons). 
As in section I.D we adopt a vector notation for these. Set 

m = (mi, . . . ,7715), e= (ei,...,e5) 

5 

N^(N,,...,Ns), N = J2Na, /x = (//i,...,//5) (7.21) 

a 

As we know (see (1.49)), the infinite volume pressure depends only on the component /x 
of the chemical potential orthogonal to the charge vector. As a consequence, the kind of 
complexes that will occur for a certain value of the chemical potentials fi depends only on 
the value of fi. We will use this freedom to make a convenient choice of v along the proof. 

The study of the molecular limit of P(/3, fj,) (i.e. /9 — > 00, /x fixed) will be carried out 
under the basic assumption on fj,: 

(A) there exists k. > {k independent of N) such that 

Hf^- H-N > kN for all N (7.22) 

The complexes (atoms, ions, molecules) that will be formed in this limit are defined as 
follows. Consider the lowest energy state of the A'"-body Hamiltonian i^N 

£;n = inf spectrum(i/N) = inf{(*, H^i^) : * G Hjv, ||*|| = 1} (7-23) 

and set 

Eifi) = inf (E^-fi- N) 

A value for which this infimum is taken, i.e. 

E^.-fi-N.^E{fi) (7.25) 

is called the composition of the complex (j) having ground state energy Ej = Ej^.. Notice 
that (7.22) imposes the bound Nj^a < K,~^E(fjL) for all a, so complexes are made out of a 
finite number of constituents, hence their number g is also finite. The complex (j), made of 
Nj^a, a — 1, ... ,S, elementary constituents, has total particle number, mass and charge 
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(7.24) 



s 

Nj = ^ Nj^^, rrij = ni- Nj, = e • Nj (7.26) 

a 

Strictly speaking, has no ground state because of translation invariance. However for 
N >2 one can decompose i^N — + H^^ into the sum of the kinetic energy of the center 
of mass and the relative Hamiltonian, expressed in the coordinates relative to the center of 
mass. Since the spectrum of starts from zero, one has also £'n = inf spectrum(ifN^). 
Consider now the threshold a^H^^) of the continuous spectrum of H^^. We know from the 
HVZ theorem (Reed & Simon 1979) that 

a,(H'^')^ inf (E^^+E^^), N>2 (7.27) 

where the infimum is taken over all partitions of the particles into two clusters a and b 
having Na and A*";, particles, and Na = 1, . . . , N — 1. From the very definition of E{fi), one 
concludes from (7.27) that a^H^ - M " N) > 2£^(/x), N > 2. This means that the part of 
the spectrum of H^^ — /x -N located in [E{iJ,), 2E{n)) consists of isolated eigenvalues {E{n) 
is strictly positive) . In particular Ej^. has to be an eigenvalue of if^' and the corresponding 
eigenstate is the the wave function of the bound complex (j). When N — 1, complexes 
reduce to the individual electrons or nuclei themselves with pure kinetic energy and we take 
then i^N = 0. 

When there are charged complexes, we add to (A) the assumption 

(B) if a complex has a charge Cj ^ 0, there is another complex having a charge of sign 
opposite to Cj (but not necessarily equal to —Cj) . 

insuring that an overall neutral gas can be formed. 

We describe now the mixture of free gases made of these complexes. To each complex 
(j) one associates the density of a perfect gas of point particles 

^'^^'^^ = (2^) exp[-/3(£;, - /X • N,)] (7.28) 

and define the pressure of an ideal mixture of such gases by 

9 

The corresponding number density of electrons and nuclei is 

5 Q 9 
Pideal(/5,M) = 5^-— Pideal(/5,/i) = iV,p,- (/3, /x) (7.30) 
a=l j=l 

If some complexes are charged, Pidcai(/3,M) for a general may belong to a non-neutral 
ensemble. Then, under assumption (B), the neutrality J2j^jPj = is implemented by 
setting 
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-Pideal,neutral(/3, = inf Pideal(/3, + i^e) 



(7.31) 



Then Pideai,neutrai(/3, /i) = -Pideai,neutrai(/3, A*) Verifies (1.49), and one can formulate the theorem 
on the molecular limit as 

Theorem Let fi be a chemical potential satisfying the condition (A) and determining com- 
plexes (j), j = l,...,g, as described above, fulfilling the condition (B). Then there 
exists £ > such that 

P{(3, fx) = Pideal,neutral(/9, M)(1 + 0(exp(-/?£))) (7.32) 

as (3 ^ oo and /i — >• /z. 

The theorem deserves several comments that will be again illustrated in the e-p system. 
One should first appreciate that all the residual effects due to thermal excitations and 
interactions between complexes when the temperature and the densities are not strictly 
equal to zero are controlled by the 0(exp(— /^e:)) correction. The issue of the theorem is 
precisely about obtaining this control. 

Concerning the assumption (A), it is a difficult task to find what are the complexes 
associated with a chemical potential /j, and to verify (7.22). Conversely, if we ask for the oc- 
currence of a certain chemical species, it is still beyond the present possibilities to rigorously 
determine the range of /j, (if any) allowing for the formation of this molecule. Suppose that 
we ask for the formation of hydrogen atoms in the e-p system. This means that there should 
be an interval G (/Ui,/U2) such that (A) holds and inf(Arg^Arp)^(o,o)(-E'iVeArp — A*(A^e + ^p)) is 
taken at {Ne,Np) = (1, 1) with E{/j) = Ea — 2/j. This is equivalent to the fact that the 
inequalities (7.18) hold when fi G (/xi,/i2) (if (7-18) is true for such a fi, one can pick k, > 
small enough so that the same inequalities are true for n + k; then k, provides the constant 
in the lower bound (7.22)). A httle bit of thinking reveals that the validity of the inequal- 
ities (7.18) for a range of /i {/i slightly above £"„) is equivalent in turn with the following 
statement about the stability of matter: there exists a constant B with < B < \Ea\ such 
that 

Hn^n, > -B(N, + Np-1) for all {N„ N,,) ^ (0, 0), (1, 1) (7.33) 

The point in (7.33) is that the constant B can be chosen strictly less than \Ea,\ for all cases 
except of course for the hydrogen atom itself. Therefore a complete proof of the existence 
of the atomic phase relies on exhibiting a sufficiently small stability constant. Although 
(7.33) must hold on experimental and numerical grounds, a rigorous proof has not yet been 
provided (see (Fefferman 1986) for a more detailed discussion). 

Taking from now on (7.33) as a working hypothesis, we recover easily the previous cases 
(i)-(iii) when applying the theorem to the e-p system. If n is shghtly above Ea, there is a 
single complex, the hydrogen atom (a) with E{fi) = Ea — 2(1 (case (ii)). If = Ea, we have 
three complexes (e,p,a) and E{/j) = \Ea\ (case (iii)). If < Ea, the complexes are (e,p) and 
E{iJ,) = \n\ (case (i)). When there are charged complexes, the minimizer u of PideaiiP, fJ'+i^e) 
has the value (7.4) so that the result of the theorem is precisely the equation of state (7.6) 
up to an exponentially small correction. 
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For arbitrary nuclei, the theorem describes generically similar situations involving more 
complicated atoms, ions and molecules. If /x selects a single complex, we have a free gas of 
that chemical species. If fx selects several complexes, we can have coexisting free gases as 
T ^ 0, generahzing the situation (iii) of the section VILA. As a corollary of the theorem, we 
can specify the relative proportion of these gases by taking jj, a function of the temperature 
that approaches hnearly as T — > (as in (7.8)) 

^iJ, + X(3-' + o((3-') (7.34) 

The densities of the gases will have weights proportional to exp(A • Nj), leading then to 
general Saha formulae. 

Although the determination of the kind of complexes occurring for a given value of fj, 
remains unsolved, it is pleasing to see how the thermodynamics of ideal substances emerges 
in principle from the basic statistical mechanics of electrons and nuclei in the molecular 
limit. 

C. Elements of proof 

The theorem was first proven (Fefferman 1985) for the e-p system under the assumption 
(7.33) (with a weaker control 0{f3~^) of the error term). All the main ideas and subsequent 
developments summarized here have their roots in this work. 

The general strategy of the proof does not rely on low density expansions, but on an 
analysis of typical equilibrium configurations of dilute gases of atoms and molecules. It 
proceeds along the following qualitative ideas. Let fi be fixed and satisfy the condition 
(7.22), and take /? very large. Then, in such configurations, if we look into a region of linear 
extension R in the range (7.12), according to the same discussion as that following (7.12), 
this region is empty with high probability, but if it contains something, this will essentially 
be one of the complexes determined by fi. The idea is now to decompose the total domain 
A = IJ^. Dr into subdomains Dr of diameter R and to let R grow with /? in a suitable way. 
One has to solve the two following problems. 

Problem 1 Determine the size of Dr such that the partition function S/)^(/3, ft) of a single 
subdomain is that of a system having at most complexes of the type (j), j = 1, . . . , g, 
in Dr. 

Problem 2 Show that the residual interaction between the different regions Dr is negligible 
as /? — > cxD so that the total partition function ~ Ylr'^DriP, IJ') is close to 

that of independent subdomains Dr. 

We make the preliminary decision (valid throughout the whole proof) to choose u = i^o(/5) 
as the minimizer of (7.31) insuring the neutrality of the ideal gas. This is equivalent to the 

pseudo-neutrality condition met earlier on several occasions. If some complex is not neutral 
and the condition (B) holds, uoiP) exists, is unique and I'oiP) = as can be checked 

from the explicit formulae (7.29) and (7.30) (see (7.4) in the e-p system). If all complexes 
are neutral one takes z/ = 0. Thus in the sequel we keep in mind that 

/z = /x + 0(/?-^)e (7.35) 
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1. Towards a solution of problem 1 

We envisage now problem 1, following mainly (Conlon et al. 1989). We write the grand- 
canonical partition function for a finite region D (ZM.^ oi volume |D| as a sum on all particle 
numbers 

g 

5^(/3,m) = 1 + $^Sb(/3,M,N,)+ Y1 ^DW,fi,m (7-36) 

i=i Nt^Nj 

where we have singled out the contribution of the g complexes that correspond to the chosen 
value of n and 

Sd(/3, N) = Trn exp {-(3 {Hd,t, - M " N)) 

= Trn exp {-f3K^^^) Tm exp (//-'^ - M ' N)) 

(7.37) 

In (7.37), Tin means that the trace is taken on states having the Dirichlet conditions appro- 
priate to the finite volume Hamiltonians Ho as well as the appropriate particle statistics. 
For R large enough, one expects that the center of mass contribution is close to that of a 
free particle in space, and that the finite volume ground state energy Edj of the complex j 
is close to its value Ej in infinite space. Indeed one can establish 

TV„exp {-PKS:^;) = (^)"' PI (l + O (i)) (7.38) 

and 

Now let 

\D\ = \Dp\ ~ e''^, c> (7.40) 

grow exponentially fast as the temperature tends to zero, and consider the contribution of 
the ground state of -ff^Jf^N. to the partition function (7.37) for N = N^. Taking (7.38), (7.39) 
and (7.40) into account, this contribution is 

= ft(/3,M)|£>;j| (1 + 0(6-!")) (7.41) 

By (7.25) and (7.35), the density of the ideal gas pj{P,fjL) (7.28) decreases as 
^-3/2g-/3£;(/x)_ E{fj,) > 0; thus one can choose c > small enough so that (7.41) is still 
exponentially small for (3 large. 

We like to show that, under the condition (7.22), all the other terms in the sum (7.36) are 
negligible compared to (7.41), provided that the rate of growth of |Z^/3| is suitably chosen. 
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(7.39) 



As \D\ — > oo, part of the spectrum of becomes continuous, corresponding to the 

scattering states of all the sub-complexes belonging to i^^^N- The basic idea is to ensure 
the convergence of the traces by the pure kinetic energy: one will dispense with the full 
Hamiltonian with the help of (7.22) and the ensuing gap property. Then the strict positivity 
of the constant k, in (7.22) and of the gap will be used to control the volume dependence of 
these traces. For this one borrows a small fraction of the kinetic energy writing 

H^ = r]K^+{l-7])K^ + V^, r]>0 (7.42) 

where = — Y2iLi \^ ^^^^ Vn is the total Coulomb interaction. Clearly the Hamilto- 
nians (1 — ifiK-^^f + V-^ and (1 — ri)^^{K-!<i + Vn) = (1 — v)~^H'n are related by the scaling 
transformation Xi — > (1 — r])xi. Hence they are unitarily equivalent. This implies, using 
(7.22) and(7.35) 

(1 - 7])K^ + V^N - M • N > (1 - ^?)"^inf spectrum(ifN - M • N) 

+ (l-r7)-V-N + 0(/5-^)Ar 
> (1 - 7])-\k - r^max \i^i\)N + 0{(3-^)N = c^N 

i 

(7.43) 

Thus 

- M • N > ryi^N + CiAT (7.44) 

where Ci > provided that t] is small enough and f3 is large enough. The same inequality 
holds when Hf^ is replaced by -f^D,N since Dirichlet boundary conditions increase the kinetic 
energy; hence one concludes that 

N . .3/2 

Sd(/3,M,N) < eM-Pc^N)Trne-^^^^ < \D\^l[ (^^^^j e"^'^^^ (7.45) 

If now \D\ — iDpl grows at the rate e'^^ with < c < ci, Eo^iP, A*, N) is bounded by e~'^^^^ 
for some C2 > 0. Thus by choosing Nq large, one can make the exponential decay of the 
whole sum of terms with N > Nq in (7.36) fast enough to have 

J2 S^,(/3,M,N)=p,.(/3,/i)P/j|0(e-^^), £>0 (7.46) 

f^,N>No 

For the remaining terms, i.e. the terms with N 7^ N^-, N < Nq and the contribution 
of the excited states to S£)(/?, /it, N^), one uses essentially the fact that E{fj,) is separated 
from the rest of the spectrum by a gap, according to the discussion following (7.27). On the 
subspace orthogonal to the ground states of H^^ — /j, ■ N • one has, arguing as in (7.43) 

(1 - 77)Xn + % - M • N > (1 - v)~\E(tJ,) + 25) 

+ (l-77)-V-N + 0(/3-')A^ 

> {l-r)y\E(iJ,) + 25-r)msiK\iii\N) + 0((3-'^)N 

i 

(7.47) 
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for some S > 0. Since now N < Nq, one can choose rj sufficiently small and /3 sufficiently 
large to have on this subspace 

Hf^- fi,--N>riK^ + E(n) + 5, 5 > 0, N < No (7.48) 

Hence for N 7^ N^-, N < Nq, one finds as in (7.45) 

N . X 3/2 

Since |-D^|^ grows at most as e'^^^^^ we can again have 

5D,(/3,M,N) = p,-(/3,/x)|L>;5|0(e-^^), N^N-, 7V<7Vo, £>0 (7.50) 

by taking c small. Notice that a volume factor \Df}\ corresponding to the free motion of the 
overall center of mass has been explicitly maintained in the estimates (7.46) and (7.50) (at 
the expense of a possibly smaller constant 5). Finally one deals with the excited states of 
— ju, ■ in the same way. Combining this with (7.41), (7.46) and (7.50) in the grand 
canonical sum (7.36), one obtains 

9 

H^,(/?,/x) = l + J2PjW,E)\Dp\ (l + 0(e-^)) (7.51) 

i=i 

Since the second term in the right hand side is exponentially small as /? ^ 00 (see the 
comment after (7.41)) one has also 

lnS^,(/3,M) = (1 + (e"'^)) (7-52) 

Taking into account the definition (7.29) together with our choice of v, this says that the 
pressure PopiP, A*) in Dp is given by 

PPoMt^) = ^lnSz,,(/3,M) = /3Pideai(/9,/i) (l + 0(e-^^)) (7.53) 

This solves problem 1: in a volume {Dg] ~ e'^^ with c > sufficiently small (depending only 
on fj,), the pressure is exclusively due to the non-interacting complexes (j), j = l,...,g. 
However this does not prove the theorem yet. In (7.53), Dp is constrained by the growth 
condition (7.40) and not an independent thermodynamic variable as it should be in the 
grand canonical ensemble. One needs to solve problem 2. 

2. Towards a solution of problem 2 

Problem 2 is technically considerably more sophisticated, so let us give first some very 
heuristic ideas. Decompose = [J^ D^ into a disjoint union of domains Dr with character- 
istic functions xor, XOri^) = 1 if x is in Dr and XOri^) — otherwise. One wants to show 
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that the interactions between the different domains D^. become neghgible as /3 — > oo. To 
this end, define an uncorrelated Hamiltonian by 

^uncor ^Hk- (7.54) 

where V^^"^ is the potential energy between the different domains 

U'°' - Y.('-^e^y''°'{xi,Xj) (7.55) 



\Xi — X2\ 

U^^ (as well as H^^'^°'^) will be an effective temperature dependent energy via the condition 
~ e^^ determined in problem 1. Note that if — a;2| < R and both Xi and X2 belong to 
the same subdomain, then V^°^{xi,X2) = 0; also V'^°^{xi,X2) is Coulombic at large distances. 
This resembles the properties of a regularized Coulomb interaction. Suppose that we are 
able to define the uncorrelated Hamiltonian in a less crude way than in (7.54) and (7.56) 
with the effect of replacing (7.56) by a smooth positive definite regularized potential, finite 
everywhere for |xi — X2I < it! and asymptotically Coulombic . A possible candidate is 

^ J J \xi - yi - X2 + y2\ 

where (pR^y) = R^^(j){^) with (j){y) a smooth function, {(piy) = 0, \y\ > 1, / dy(j){y) = 1) 
representing a charge density supported in a region of extension R with total charge equal 
to 1. Other candidates closer to the form (7.56) are 

T/:-(xi, X2) = I ^ ^1 ^ (7.58) 

with h{x) a smooth short range function, h{0) — 1, and ^— j^j^ positive definite. Then we 
can use the basic positivity argument already introduced in section LB 



t^reg = I E ea,e„,y-(x„ Xj) - ^^-(O) ^ ^ 

ij i 

1 ^ 



i=l 



because of the positive definiteness of V^^^{xi, X2). Note that by scahng, in both forms (7.57) 
and (7.58), one has Keg'^(O) = O (-^) . Therefore if one can find such a regularized correlation 
energy, one will have (being still very sketchy) that H^'^°^ provides a lower bound to the 
full Hamiltonian 

//a - • N = nr"^' + U^l - /A • N > Hr^' - (1 + 0{R-^))fj, ■ N (7.60) 
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Thus, up to a small correction to the chemical potential, one finds that the partition function 
is dominated by that of the uncorrelated domains 

Ha(/3,m)< n ^Dr{P:Kl + 0(R-'))) (7.61) 

We set now 1-0^1 = \Dr,f3\ ~ e'^^, (3 sufficiently large, and apply to each of the Dr the result 
(7.52) of problem 1. Taking the logarithm of (7.61) and absorbing the error in the chemical 
potential in the 0(e^^^^) correction, one obtains in the thermodynamic limit 



/3P(/3,iu)< hm 

— |A|^oo 



r,f3\ 



/3Pideai(/3,/*)(l + 0(e-^^)) 



= /3Pideai(/3,M)(l + 0(e-^^)) (7.62) 

The technically most elaborate parts of the works (Fefferman 1985, Conlon et al. 1989, 
Graf & Schenkcr 1995 &) consists in constructing an uncorrelated Hamiltonian having the 
requested positivity and smoothness properties, and comparing it with the true Hamiltonian 
H\. In its main lines the construction of (Conlon et al. 1989, Graf & Schenker 19956) goes 
as follows. 

The decomposition = |J^ Dr is reahzed with regions Dr that are translates and dila- 
tions (with the scale parameter R) of a fundamental domain Dq. Let T^^, be the Fock space 
for a Coulomb system localized in Dr, having Hamiltonian i//)^ with Dirichlct boundary 
conditions on dDr- The most obvious way to define states of independent Coulomb systems, 
each of them confined in a region Dr, is to introduce the product space JFJ^'^'^" = H^^nA^S -^^r 
and the Hamiltonian JZornA^^ on it. One wants to compare this obviously strictly un- 
correlated Hamiltonian with the true Ha acting on when the scale parameter i? — > oo. 
For this, one introduces a partition of the unity on of the form XlrXri?!^) — 1 
X G where Xr,R are suitably chosen smooth functions essentially supported in the region 
Dr of linear size R. With the help of this partition, one defines an isometry J : J-'\ — > J^^^^'^'^ 
which maps the states of the system onto the uncorrelated ones. Then one compares the 
image J~^i'^DrnAj^$-^Dr)J of the strictly uncorrelated Hamiltonian with H\ on The 
calculation gives 



(7.63) 



where U^^I^ has the form (7.55) and (7.58), and 



h{x) = J dyx\x + y)x\y) (7.64) 

is a smooth function localized in the neighborhood of Dq: (7.63) is a more precise version of 
(7.54). The O (-^) term is the price of the increase in kinetic energy paid for the additional 
Dirichlet conditions imposed at the boundary of each Dr in the strictly uncorrelated Hamil- 
tonian J2r Hor- This price can be made small as i? — cxd because of the differentiability of 
the localization functions Xr r- 
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The main issue is now the lower bound (7.59), i.e. to have ^^^"'^(xi, X2) positive definite. 
In (Conlon et al. 1989), the regions Dj. are cubes, and this positivity is obtained by replacing 
the Coulomb potential in ^DrnAj^$ -^Dr by a Debye potential ^-j^ with small £ (e ~ e~'^^). 
Then it turns out that 

is similar to Voo,x (1-8). Its Fourier transform can be shown to be positive for e and R in 
appropriate ranges by using properties of h{x) (7.64), and the arguments (7.59)-(7.62) can 
eventually be cast in a rigorous form. However, to complete the proof, one has to show that 
the assumption (A) in the main theorem as well as the results obtained in problem 1 are 
stable when the Coulomb potential is approximated by the Debye potential. In particular, 
if (A) holds with a constant then one has also for the Debye Hamiltonian 

- A* • N > K{e)N with K{e) ^ k, e^O (7.65) 

This theorem is proven in (Conlon et al. 1989) with the same localization techniques 

In (Graf & Schenker 19955), the decomposition = IJ^ -Dr is realized with simplices. 
Then the stability (7.59) follows remarkably from the geometrical properties of a simplex. 
Namely for the spherical average h{\x\) of h{x) (7.64) (here x = Xdq is the characteristic 
function of a simplex), the function has a positive Fourier transform. 

The original proof (Fefferman 1985) uses certain coverings of (called Swiss cheese) 
made of small cubes and large balls of different radii, and a positive definite potential of the 
form (7.57). The delicate point is to estimate the difference V^"°^ = y^™' — V^j.™'' between 
(7.56) and (7.57). V^"°^ still contains a local singularity as well as non-positive definite 
Coulomb energy terms. By averaging on the different radii of the balls, Fefferman shows 
that these dangerous contributions in V^"^°^ are majorized by a Coulomb potential ^^o-^ 
with a small amplitude, and that these dangerous terms become eventually negligible. Then 
the stability (7.59) follows from the positive definiteness of Vj,^™ (7.57). 

In all cases the proofs involve fixing a number of important details not mentioned here. 

Finally, (7.62) must be completed by the converse inequality; the job is fortunately easier 
and relies on an application of the variational principle. We refer to the original papers for 
this final step. 
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VIII. CONVERGENT EXPANSIONS 



In previous chapters we have alluded to the Mayer expansion and results on the conver- 
gence of Mayer expansions. This is an old subject reviewed in numerous places including 
(Brydges 1984) but since that review some beautiful new combinatoric formulas have been 
discovered. The main objective of this summary is to update (Brydges 1984) to bring these 
to the attention of our readers, together with some applications and ramifications. This is 
done in sections VIII.A, VIII.B. 

Related ideas convert the considerations of section II. G into convergent expansions and 
thereby prove there is screening. We outline these arguments in sections VIII. C, VIII. D, 
VIII. E. This updates an older review (Brydges & Federbush 1981). 

We arc accustomed to expansions in theoretical physics where the terms are labeled by 
connected graphs. It is usually the case that the number of connected graphs at order N 
grows more rapidly than A^! so that such expansions appear to be divergent at least until 
cancellations between graphs have been taken into account. The older review (Brydges 1984) 
was centered on rigorous bounds that say that the sum over all connected graphs at order 
N is actually smaller than the sum over just the least connected graphs, i.e. tree graphs. 
The new part of this review is theorem VIII. 2 which makes the relation between all graphs 
and tree graphs much clearer. 



A. Tree graph formulas and the Mayer expansion 

We begin with a little review of the theory of activity expansions (Ruelle 1969, Hansen 
& McDonald 1976, Mayer & Mayer 1977). Suppose particles labeled 1,2,...,N have 
two-body potentials Vy. The total potential of a subset X of these particles is 

Ux=Yl (8-1) 

where ij = ji is an unordered pair of distinct particles and ij e X means i,j e X. Consider 
the formulas 

= E n - 1) 

G ijeG 

E IK"-""-!) (8-2) 

G connected ij&G 

In the first the sum is over all graphs G, connected or disconnected, on vertices in X. It is an 
immediate consequence of expanding the product over ij of exp {—Vij) — 1 + 1, remembering 
that a graph on X is, by definition, an arbitrary subset of {ij : ij G X}. The second formula 
is the usual way to define the connected part exp (— t/x)^, but there is another equivalent 
definition, namely the connected part is the (unique) recursive solution of 

--"'-Em E n^^-"^'! 'fi^i>i 



M>1 Xi,...,Xm j=l 
UXi=X 



c 



(e-^^)^ = 1 if IXI = 1 (8.3) 
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where Xi, . . . , Xm are disjoint. The equivalence of the two definitions is not a hard exercise 
and is the main step in obtaining the Mayer expansion, which is 

Theorem VIII. 1 Let Un — U{i,...,n} where Vij — V{ri,rj) is the pair interaction between 
particles in a grand canonical ensemble S. Then 

1 r ^ 

i"-=E]^/n^^.-^(^.)(^~'''^)c (8-4) 

N>1 ' j=l 



The proof that this formal expansion results from (8.3) is in most textbooks in statistical 
mechanics. The classical proofs that this expansion is convergent for small activity were 
given in (Groeneveld 1962, Penrose 1963, Ruelle 1963). Our line of development begins 
with (Penrose 1967) who showed that when the interactions are repulsive the sum over all 
connected graphs at order is dominated by the sum over tree graphs at the same order. 
By Cayley's theorem the number of tree graphs on vertices is N^'"^ which is comparable 
to A^!. The convergence of the Mayer expansion at small activity will be seen to be an 
immediate consequence. 

Notation: We consider the class T of graphs called forests: A graph G is a forest if it has 
no closed loops or equivalently if each connected component of G is a tree graph. A single 
isolated vertex which is not in any bond is considered to be an empty tree graph. The empty 
graph which has no lines, so that every vertex is isolated, is therefore a forest. A forest G 
has the property that for any pair of vertices ij, either there is a unique path that joins i 
to j consisting of bonds in G or there is no path at all. To every bond ij e {1,2,..., N} 
is associated a non-negative parameter Sij. The vector whose components are all these 
parameters is denoted by s and U is the vector where each Sij — 1. Let G be a forest. For 
each bond ij we set 

. „ ^ _ ( infjsf, : b G path in G joining i and j} , . 

^^^^^''^-\Oifnopath ^^■^> 

Thus if ij belongs to the forest G then — Sy. Also, aij{G,s) depends only on the Sij 
parameters assigned to bonds in G. We let cr{G, s) denote the vector whose components 
are aij{G,s). We set d'^ = YlbeG^/^^b- The symbol cPs means that for each b E G Sf, 
is integrated over the interval [0, t]. We write / o a{G, s) = /(cr(G, s)) so that the notation 
(d'^F) o cr(G, s) will mean: first do the derivatives with respect to with b e G and then 
evaluate at <t(G, s). 

The following theorem is a type of fundamental theorem of calculus. It reached this 
simple form through the series of papers (Glimm, Jaffe & Spencer 1974, Brydges & Federbush 
1978, Battle & Federbush 1984, Brydges & Kennedy 1987, Abdesselam & Rivasseau 1995). 

Theorem VIII. 2 Let F(s) be any continuously dijferentiable function of parameters 
where ij — ji and i 7^ j e {1, 2, . . . , N} and N >2. Let t > 0. Then 

F(tll) = J2f ° ^) (8.6) 
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where 

J d^s{d^F)ocr{G,s) = F{0-Sl) 

when G is the empty forest. 

Note that the theorem is vahd when t = because all terms in the sum over G vanish 
except when G is the empty forest and in this case the right hand side is F(OII), by defini- 
tion. Also, \i N — 2, the theorem reduces to F{t) — F(0) + /* 9i2-F(si2)dsi2, which is the 
fundamental theorem of calculus. A proof of this theorem is given in section VIII. F. 

To study the Gibbs factor cxp {—[5 ^ Vij) by this theorem we introduce parameters Sij 
in such a way that the Gibbs factor is the value of 

F(s) = e-^^-(^) = exp (^(3 j (8.7) 

when s = U. We choose Vij{sij) to vanish at Sij = so that s = corresponds to all inter- 
actions being switched off. We will call such a choice of dependence on sij an interpolation. 

Assume we have sTich an interpolation. Given a forest G on vertices {1, 2, . . . , A^}, we 
can decompose it into connected components. These connected components are subgraphs, 
trees, on subsets of vertices Xi, . . . , Xm- Thus G determines a partition of the set of vertices 
into "clusters" Xi, . . . , Xm- Vertices label particles and from the definition of <t = cr(G, s) 
there are no interactions in UN{cr) between particles in different trees so that (d'^F) o cr 
factors across the partition Xi, . . . , Xm- By comparing the definition of connected part (8.3) 
with the result of the theorem one obtains 

Theorem VIII.3 For N >2 

{e-^u^)^ -Y.f {d^e-''''^) o a{T, s) (8.8) 
where T is summed over all connected trees T on N vertices. 

Tree Graph domination: Suppose that the potential is repulsive, Vij > 0. Then let T^j(sy) = 
SijVij or indeed any other interpolation with V-j{s) > 0. By 

Um{ct{T,s)) > V.,(a,,(T,s)) = J2 y^Ji^^J) (8-9) 
ijeT ijeT 

and theorem VIII.3 

I (e-'^^)c I < E /' n ■(%)e-^^'^-(-^-) - E - ^"''^") (8-10) 

T ijeT T ijeT 

which says that the connected part — the sum over all connected graphs — is dominated 
by the sum just over tree graphs, which is the remarkable result (Penrose 1967). 

However there are also tree graph domination results when the interaction is not purely 
repulsive. <t(G, s) has the following unobvious property to be proved in section VIII. F. 
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Theorem VIII.4 // 



Y,V,^>-BN (8.11) 

then 

^aij{G,s)Vij>-BN (8.12) 

ij 

The theorem remains valid when self-energies 1/2^^^ V^j are added to both left hand sides. 
Tree Graph Domination (8.10) now becomes 

\{e-^^^)J<e^'^J2Uf^\^^^\ (8-13) 
T ijeT 

Tree Graph limit: Consider the limit 

z^oo, /3 ^ with fixed (8.14) 
By theorems VIII. 1 and VIII. 3 the leading terms are 

1 r ^ 

AT ' T k=l ij£T 

where T is summed over all connected trees T on N vertices. There are — 1 lines in 
a connected tree graph on N vertices so every term in the sum is 0{z^ j3^~^) = 0{z). 
The estimates in section VIII. B show that the Mayer expansion converges uniformly [after 
dividing out an overall factor of ^ x volume] as this tree graph limit is taken so the limit 
can be taken under the sum over A^. This series can even be summed exactly. It is identical 
with the series in powers of z for 

^ J h{S)V-^h{S)dS + J z{S) d£ e-^^^^^ (8.16) 

evaluated at the h that makes it stationary. is the operator inverse of of the operator 
whose kernel is V . This is the well known folk theorem that Feynman tree graphs sum to 
the classical action. 

More general interpolations besides SijVij are actually useful because sometimes one has 
interactions which are really the sum of several interactions at different scales and then it is 
useful to use more complicated interpolations that turn them off one after the other. Since 
there will be examples of this in the next section we also record a more general stability 
result 

Theorem VIII. 5 If the interpolation satisfies 

J2v:,{t) > -B{t)N (8.17) 

ij 
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then 



f/„(<T(G,s)) > 



■2E 

ijeT 



B{s)ds 



U8) 



Tree Graph Domination (8.10) now becomes 

I (e-^^^)cl < E fd^' n mM^)\e''^^^"''^'^'' (8.19) 



B. Convergence of Mayer expansions 

We demonstrate some uses of tree graph domination by giving estimates on the radius 
of convergence for the Mayer expansion for a variety of systems including ones discussed 
earher in this review. 

We will use the following two principles: A tree graph on N vertices whose lines represent 
factors /(r — r') contributes exactly 

j f{v)dv^ X Volume (8.20) 

because tree graphs, being connected, have — 1 lines. The integrals over vertices can be 
evaluated in order starting with ones on lines that are outermost branches (Trim the tree!). 
Furthermore the number of tree graphs at order N is iV^~^ (Cayley's theorem). 



1. Repulsive potentials 

As an instructive exercise we recover the result in (Penrose 1967). Consider for simplicity 
one species. By tree graph domination (8.10) the Mayer expansion is bounded by 

N ' '' 

^Ei^E(w /(l-e-'"''")<ir)"'|.||A| 

AT ■ T ^ / 

l-l / (1 - e-^^(^^)rfr ) N|AH^^Q^-N|A| (8.21) 

N ' ^ / N ' 

which is convergent provided 

Q := 1^1 J dr{l- e-^^W) < e"^ (8.22) 
This estimate is also used in the next two examples 
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2. Debye spheres 



In section II. D we encountered conditions (2.44). The system consisted of charged par- 
ticles interacting either by two body potentials 

VL,x{r) = |r|-i (e-l"-!/^ - e-'^-l/^) (8.23) 

or by the Yukawa potential VL{r) and hard cores of radius A. First we consider two species 
interacting by (8.23). Suppose we can find an interpolation as in theorem V111.5, then we 
can repeat the argument given above for repulsive potentials, but now using the tree graph 
domination (8.19) instead and find convergence when 

Q := 2\z\ J dr\pV'{r,s)\ exp ^2/3 B{s')ds'^ < e'^ (8.24) 

The factor 2 is there because there are two species. 

The interpolation is chosen to be V{r,t) = Vi.^(t)(r) where i{t) is a length scale which 
decreases from L when t = to A when t = 1. Then we show below that 

f) f,-\r\/m 

V^'(r,t) = -^^^ (8.25) 
obeys a stability estimate (8.17) with 

Also 

V\r,s)\dr^-g / '-^dr ^ -^^Aniisf (8.27) 



/ 



because |y(r, s)| = V'{r,s). Therefore (8.24) reads 

^16n\z\p deeexp(^(^-^^ <e-^ (8.28) 

When X < P < L we break up the range of the d£ integration into [A, /3] and L] and find 
that this criterion is satisfied if 

\z\(3L'^ < 1, Izlp^ef^/^ < 1 (8.29) 

which is the condition (2.44). The Ruelle estimate (2.43) is not as good because it does 
not distinguish the length scales and essentially associates the bad stability factor exp(/?/A) 
to all scales. In the Debye-Hiickel limit Q = 0{L'^ /l^) is determined by the first criterion 
because \z\(3'^e^/^ — >• 0. By this remark and (8.21) estimate (2.46) follows. Results related 
to (2.44) were given in (Gopfert & Mack 1981, Imbrie 19836, Brydges & Kennedy 1987). 
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To obtain the stability estimate (8.26), note that by the argument given near (1.12) the 
stabihty constant B{i) for V^g is half the self-energy, which can be calculated using the 
Fourier transform to be 




(8.30) 



Now take L in (8.30) equal to £ + di . 

It is not a hard exercise to rework this calculation for two dimensions. One recovers a 
result of (Benfatto 1985) that ii P < 2 [ An using the units in (Benfatto 1985)] that the 
Mayer expansion is convergent for \z\ small uniformly in the short distance cutoff A. In 
other words for (3 smaller than the first threshold the Mayer expansion converges for the 
two dimensional Yukawa gas without any cutoff. This is notable because this interaction is 
singular at r = and is not stable in the Ruelle sense. The reader will find this exercise 
done in (Brydges & Kennedy 1987). 

The Yukawa gas with hard core is convergent under the conditions (Brydges & Kennedy 
1987), see also (Imbrie 19836), 

|^|A'<1; \z\PL^<.l; < 1 (8.31) 

Briefiy, this is derived by combining the hard core example with the Yukawa by letting s 
vary over [0, 2] so that V{r, s) has no interactions at s = and has the complete interaction 
when s = 211. For s G [0, 1] we choose V as in the repulsive example and for s G [1, 2] we 
have V given by (8.25) with i{s) = L at s = 1 and i{s) = at s = 2. The main point 
is that the hard core is already there when the Yukawa is turned on so that 2B{s) for s 
corresponding to l{s) ~ is the stability constant including the hard cores, which prevents 
a divergence. 

Tree graph formulas were applied iteratively by (Gopfert & Mack 1981) to expand succes- 
sive scales in the interaction in a renormalization group analysis. They showed how to adapt 
stability estimates to scales and particular clusters of particles. They were motivated to look 
for better estimates on the convergence of Mayer expansions in order to prove confinement 
at all parameters for three dimensional Euclidean lattice quantum electromagnetism by ex- 
ploiting a connection with screening. The iterated tree graph formulas of Gopfert and Mack 
were then simplified (Brydges & Kennedy 1987) by passing to a limit in which the scales are 
replaced by a continuous parameter (the t in Theorem VIII. 2). This was partly motivated 
by the study of renormalization (Polchinski 1984). (Abdesselam & Rivasseau 1995) then 
found Theorem VIII. 2 as a corollary of the results in (Brydges & Kennedy 1987). 

The ideas in (Gopfert & Mack 1981) were independently developed by Gallavotti and Ni- 
colo and applied by Benfatto to the two dimensional Yukawa gas (Gallavotti 1985, Gallavotti 
& Nicolo 1985a, Benfatto 1985, Gallavotti & Nicolo 1986). Note that the trees in these pa- 
pers express the hierarchical structure of the clusters on different scales and should not be 
confused with our trees which are Feynman or Mayer graphs. 

3. The Mayer expansion for a polymer gas 

This application is technical and should be omitted on first reading. It will be used later 
in this chapter. Consider the grand canonical partition function 
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^(^)=Em e n^w (8-32) 

N Xi,...,XjvCA j 

UXj=A 

where A is a finite set. In applications to lattice models A is all lattice points in a big 
box. A typical element of A is denoted by rj. Swimming around in this big box are the 
"polymers" Xi, . . . ,Xn, which are disjoint subsets of A. For lattice models they are often 
nearest neighbor connected but instead of making this assumption we assume that the 
activity A{X) of polymer X obeys a bound of the form 

|^(X)|<al^lJ] n ^(^'^') (8-33) 

G -qri'GG 

where |X| is the number of elements r] in X, G is summed over all tree graphs on X. For 
example when A is a box in a lattice, we could take w{ri, rj') to vanish whenever rj, rj' are not 
nearest neighbors and to be one otherwise. In this case A{X) — unless X is connected. 
Now we show that the Mayer expansion is convergent when 

Q:=aJ2 (V + ^iv, V')) < (8.34) 

where 5,,,,,/ is the Kronecker delta function on the elements of the set A. Suppose there is a 
notion of distance between elements in A and w{rj,rj') decays exponentially as dist (?],?]') — >• 
oo. If eQ < 1 then correlations obtained by differentiating InZ(A) with respect to external 
fields ^ with A{X) = A(X, -0) decay exponentially. 

For example 7] could be a point in a lattice in which actually represents a unit box 
centered on r] so that X represents a union of these boxes. If A[X, ip) is independent of 
variations of iIj{v) at points r outside X then 



and higher variational derivatives will decay exponentially , uniformly in the size of A as the 
points r, r' are separated. 

We drop the constraint that polymers are disjoint from the sum in (8.32) and impose it 
instead by a hard core Gibbs factor exp(— f/jv). Let 

J(Xi, . . . ,Xjv) = A{X,) ■ --AiXM) (e-^-(^i'-'^-))^ (8.36) 

In the Mayer expansion 

E J{Xu...,X:,) (8.37) 



N>1 Xi,...,XnCA 



we write 



VIII-8 



Then we set M = X^i^j ^^te that M!/(]^nj!) is the number of ways to partition 
{1, . . . , M} into subsets 71, ... , 7jv. Therefore (8.37) is the same as 

M>\ N>1 li,--;lN Vi,-,Vm 

= {rij ■ j e 7^} (8.39) 

with the constraint on the sum over rji. . . . .rjM that rjj with j in the same partition are 
distinct. We substitute in the tree graph domination bound (8.10) 

<Y.\[ |e-^('''^')-l| (8.40) 

T ijeTneXi,r,'eXj 

and (8.33). Note that the hard core interaction V{ri,r]') means that exp {—V{7],7]')) — 1 is 
—Sf^^r)'- The trees with w bonds and the tree with 5^ ,,/ bonds hnk into one connected tree G 
on {1, . . . M} and the respective sums are equivalent to a sum over all possible connected 
trees on {1,...M}. The partition {7i,...,7jv} determines uniquely which bonds in this 
tree are V bonds and which are w bonds. Indeed the sum over partitions, including the 
1/N\ is equivalent to summing over the choice w bond or V bond for each bond in the tree. 
Therefore the sum over partitions is the same as assigning to each bond ij 

K,rij+^iVi,Vj) (8.41) 

Therefore the expansion (8.39) is bounded by 

M 

^Em E En(^--.+^(^^'^.)) (8-42) 

M>1 vi,-,VM G ijeG 

dropping the constraints on the distinctness of 771,..., r^M- As above, (8.20-8.22), this is 
convergent if (8.34) holds. 

The additional claim we made concerning exponential decay of correlations is an imme- 
diate consequence: suppose we perform variational derivatives such as 



on InZ where is an external field that enters through dependence A{X) — A{X, -0). Then 
all terms in the Mayer expansion (8.37) vanish except those where the union of the sets Xj 
contains both points ria,rib, because 6A{X,ip)/ {6tp{rij) = when r] ^ X. When ria,rih are 
far apart the surviving terms are exponentially small since they are connected graphs with 
exponentially decaying propagators. This argument that is being made term by term in the 
Mayer expansion is vahd when all the terms are summed because we estabhshed that the 
expansion is convergent, provided the convergence has some uniformity in ■0. 



VIII-9 



C. Gaussian integrals and truncated expectations 



In this section, following (Brydges 1984, Brydges & Kennedy 1987) we derive tree graph 
formulas for truncated Gaussian expectations. The point of the formulas is that all the loops 
in the associated Feynman diagrams are resummed back into a positive Gaussian measure, 
leaving only tree diagrams, which carry the connectedness information. 

Truncated expectations for arbitrary observables will be defined later in (8.60) but for 
products of polynomials they can equivalently be defined by applying Wick's theorem and 
then discarding all but the connected graphs. 

Suppose ( ) is a Gaussian average over variables 0j, i — 1, . . . N, 

i 

(Mj) = A? = (8.43) 

A is a matrix with positive eigenvalues. Suppose Pj = Pi{(f)i) are polynomials each of which 
depends only on one of the variables, then by Wick's theorem 

(ni^.) = e^S..^'^J]p,|^=o (8.44) 
I I 

where the exponential is defined by power series and 

d d 

' <90i d(j)j 

We apply theorem VIII. 2 with 



^ij = Cij — — (8.45) 



F(s) = exp -'^^^^^ P (8-46) 

with P = Y\iPi and the conventions Sij — Sji and Sa = 1 . The right hand side of theo- 
rem VIII. 2 will then contain 

exp I ^ J] a^j{G, s) A,, J A^P|<^=o (8.47) 

where A'^ = Ylij^c'^ij- '^^^^ ^ Gaussian expectation ( )'^'^ characterized by 

{Mj)''''^(Tij{G,s)Cij (8.48) 

The key point is that it truly is a Gaussian expectation: we shall show below that Cijaij 
has positive eigenvalues, so that it has an inverse which defines the expectation {(piCpj)'^'^ in 
parallel to (8.43). We define 



G,s 



(8.49) 
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and conclude from theorem VIII. 2 that 



i G i 

where G is summed over all forests. 

Notice that the previous equations contained exponentials of A which only make imme- 
diate sense on polynomials, but now that we have got rid of these we can use this formula 
when the Pi are arbitrary smooth functions. 

By the discussion above theorem V111.3 there is a factorization ( = nrcG( )^ where 
T runs over trees in the forest G. Therefore by definition of the truncated expectation 
( )t given in (8.60) we obtain 

N 

i=l G i 

where G runs only over connected trees on A'^ vertices. This argument easily generalizes to 
allow each Pi to depend on more variables 0a;, x & Xi. This formula expresses a truncated 
Gaussian expectation as a sum only over tree graphs. The loops that would appear in a 
standard application of Wick's theorem are resummed into the expectation (— defined in 
terms of the Gaussian expectation with the altered propagator Cijaij. 

To see that the matrix CijCFij has positive eigenvalues: since A^j has positive eigenvalues, 
Cij has positive eigenvalues which is equivalent to 

X]ry,C,,7y,- > A„,in^7;| (8.52) 

i,j 3 

By theorem VIII. 4 with Vij — r]i[Cij — Xm.mSij]r]j, B — and self energies included, 

E^? (8.53) 

ij j 

which implies that the matrix A{cr)~^ = CijCJij has positive eigenvalues. We define the 
expectation ( )'^'^ by replacing the matrix A by A[(t) in (8.43). 

D. Polymer representations and exponential decay 

In this section we describe how to obtain a convergent expansion for the logarithm of 
functional integrals of the form 

{Z{K)e'^^f)= [ dfx{<f))Z{A,<f))e'^'f'f (8.54) 



^^A superscript T refers to a tree graph and a subscript T means the expectation is truncated. 
We apologize for the notational collision. 
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where dfi is a Gaussian measure with a propagator [covariance] C(r, r') that has both infra- 
red and uhraviolet cutoffs so that it is smooth on the diagonal and has exponential decay 
as r — r' — >^ cxD. The A in the functional Z{A, 0) means that variational derivatives at points 
outside A vanish: 

Z^{A, 0; /) ^ -^Z{X, (t> + = (8.55) 

if / = in A and the same is true for all higher derivatives. 
Functionals such as 

Z(A, 0) = exp(z / cos0(r) dr) (8.56) 

factorize 



Z{X, U ^2, 0) = Z(Xi, 0)Z(X2, 0) (8.57) 

but we will work with a less restrictive property, which is that Z{A., 0) has a polymer repre- 
sentation 

N ' Xi,...,XjvCA j 

UXj=A 

where the interiors of the sets Xi, . . . , Xn are disjoint. In other words ^(A, 0) is at least a 
sum over contributions that factorize. This type of representation is not useful if the sets 
Xj are very complex so we insist that all sets are finite unions of unit cubes in M^, centered 
on points with integral coordinates. 

Recall the well known combinatoric miracle that when the partition function has a graph- 
ical expansion, the logarithm has the same expansion except that only the connected graphs 
appear. This comes about as a result of the following easily verified relation (cumulant 
expansion) between expectations and exponentials 

(eO = exp( (8.59) 
\Ar>l ■ / 

where the truncated expectation (Fi • • • F/v)r is recursively defined by solving 

(Fi---F^) = J]J](F^)r if A^>1 

(Fi)t = (Fi) (8.60) 

where tt is summed over all partitions of {1,2,..., N} and F'^ = Y{j^^ Fj. For the case (8.56) 
the expansion for the logarithm of (8.54) based on these relations is convergent at small z, 
essentially because the cosine has the special property that it and all of its derivatives are 
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bounded uniformly in The next steps can be viewed as a reduction of more general 
functionals such as cxp(— A J 0*^) to this bounded case. 

Relations similar to (8.59,8.60) hold for the polymer representation (8.58). This is plau- 
sible once one realizes that the polymer representation is an exponential in a different com- 
mutative product o defined by 

A o B{X) = J2 A{Y)B{X \ Y). (8.61) 

YCX 

As always Y, X arc finite unions unit cubes. Y, X arc permitted to be the empty set. We do 
not distinguish between open and closed cubes, that is we identify two cubes if they have 
the same interior and regard them as disjoint if their interiors are disjoint. Define 

Sxp{K)^I + K + ^KoK + --- (8.62) 

X(X) = unless X is the empty set, in which case it is one. X is the identity for this 
product. We only define Sxp{K) for functionals K that vanish when applied to the empty 
set. It is easy to verify that 

Sxp{K) = K ^ ^ (8.63) 

Thus Z(A, (f)) — £xp{K){A, 0). Furthermore £xp{A+B) — Exp{A)o£xp{B). As an example, 
consider the function □(X) that equals one if X is a cube and vanishes in all other cases. 
Then by (8.63), for any set X, 

£xp{n){X) = 1 (8.64) 
Returning to our functional integral (8.54), let us define 

i^T(X) = J] ^(ir°^(X))r, =OifX = (8.65) 

iV>l 

More explicitly 

Xi,...,XjvCX j 

UXj=X 

where the truncated expectation ( )t is defined as above. Then Kt is the logarithm in the 
sense 

{Sxp{K) {X)) = Exp{Kt) {X) (8.67) 



^''By reversing the Sine-Gordon transformation the expansion becomes the Mayer expansion for 
±1 charges interacting by two body potential C(r, r'). 
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which is the o equivalent of (8.59). However these series are terminating after a finite 
number of terms that depends on the sot X because there are at most a finite number of 
partitions of any given X. As an extreme example, if X = A is a single cube, then there 
are no partitions of X into proper subsets so only the N — 1 term in (8.65) survives and 
accordingly Xt(A) = {K{A)). 

This analysis is unaffected by the external field term exp( J (f)f) in the partition function 
(8.32) because it can be absorbed into the polymer activities K. From the definition of Sxp 

£xp(K(cf>)){A)J-^f ^£xp(K(<l>,f))(A), K{X,<t>J) -.^ K(X,cl>)e^x^f (8.68) 

We want an expansion for the standard logarithm of the partition function. Comparing the 
definition of £xp in (8.63) with the polymer gas partition function (8.32), we see that the 
standard logarithm is given by the Mayer expansion described in section VIII. B. 3, provided 
we can weaken the constraint UXj = X that occurs in £xp{Kt){X) to UXj C X. There is 
a simple way to achieve this: for any function of sets A[X) declare that A is normalized if 
A{X) = 1 whenever X is a unit box, then 

£xp{A){X) = £xp(D)o£xp(A - D){X) 

= E]^ E £xpin)ix\uX,)l[A{X,) 

N ■ Xu-,XnCX j 

-E]^ E n^M' ^---A-n (8.69) 

N ■ Xu-,XnCX j 

The constraint UXj = X has been lifted because by (8.64) £xp{0) = 1 on X \ UX,-. If A is 

A{X) = Kt{X)/ H KriA), (8.70) 

AcX 

where / has its usual division meaning, then it is normalized. Also, by the same principle 
as was used in (8.68) 

£xp{Kt){X) = I H Kt{A)\ £xp{A){X) (8.71) 
\acx ) 

so 

£xp{Kt){X)={^Kt{A)\Y^^^ (8-72) 

\ACX / N ■ Xi,...,XivCX 3 

In section VIII. B. 3 we gave a criterion (8.33,8.34) for how small the activity A{X) must 
be for convergence of the Mayer expansion and exponential decay of correlations. A is 
built out of Kt so we need to see that obeys a bound similar to (8.33,8.34). This 
is where section VIII. C comes into play. By the tree graph formula (8.51) Kj^ is small 
provided variational derivatives with respect to of the initial polymer activity K(X^ 0) in 
the partition function (8.54) are small. We will give a partial explanation of this through 
the following example. 
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Example: suppose 



Z(A) = e-^(^), V{K,(t>) = X f (p" (8.73) 

Because this factors into a product over boxes in A, by the same argument that shows that 
£xp{0) — 1 this equals £xp{K) with 

When A is small variational derivatives of K{X, 0) are small, uniformly in 0, because 
da 

because 0^ controls J (p^g when is large. This is somewhat of an oversimphfication because 
the dependence of the estimate on the test function g is not uniform. Unfortunately a gentle 
reference to cover these specific missing details is not available, but (Lemberger 1995) is 
detailed and close in spirit to the discussion here. More general background references are 
(Glimm & Jaffe 1987, Rivasseau 1991, Brydges 1992). The ability to work with the bounded 
exp(— A J (p*) in place of J 0^ is half the reason why we obtain a convergent expansion 
when A is small. The other half is that all formulas involve at most tree graphs whose 
numbers are compensated by A^! factors, see the discussion at the beginning of section VIII. B. 
Exponential decay of correlations for functional integrals of the form (8.54) follows as in the 
discussion of the polymer gas in section VIII. B. 

It is possible that polymer representations could be useful in numerical renormalization 
group calculations. This has been considered (Mack & Pordt 1985, Mack & Pordt 1989). 
Their program can be viewed as a natural step beyond methods that do work well for 
hierarchical models (Koch & Wittwer 1991, Koch & Wittwer 1994). 



E. The polymer representation 

In this section we indicate how exponential decay of correlations for the Coulomb system 
is deduced from the considerations of section VIII. D. The missing step, outhned in this 
section, is to show that the Sine-Gordon representation is an integral of the form (8.54) 
with Z{A, (f)) = Sxp{K){A, (p) for polymer activity K[X, 0) which becomes very small as the 
Debye-Hiickel limit is approached. 

As introduced in section II. E we use units where Id — i- We also simplify the presentation 
by discussing a Coulomb system with two species with equal but opposite charges and equal 
activities z interacting by a smoothed Coulomb potential of the form Voo,l with L = Ijj. 
Then there is no need for the Mayer expansion of section II. D, which complicates the original 
proofs. 

With these simphfications the left hand side exp (F) in the "Villain" approximation 
(2.79) becomes the ideal gas partition function 

exp Q d£ z^''\S)e-^'^'''"''t'^ (8.74) 
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We make the dependence on A explicit by writing exp {F(A, i(f))). The left hand side of the 
Villain approximation (2.79) factors 

^F{XUY,ict>) ^ ^F{X,i4,)^F{Y,i4,) ^g y^^j 

whenever X and Y are disjoint. The sum over h in the right hand side of the Villain 
approximation also factors, if X, Y are unions of unions of cubes. From this there follows 
the same property 



for R3. 

Prom (2.80) 



SA/(5ideal,A5DH,A) — / C^/^Dh(0) ^6 



^ -SY)ii(-ih,ig) 

h 

y^^-S' Y3Yi{-ih,ig;i4))+Rz{ig+i4>) 



For simplicity we set the external field = 0. In the formula (2.82) for S-D^l—ih, ig), R2, Rq 
vanish for this simplified model. The propagator for c?/iDH(0) decays exponentially which is 
the property we need in order to use the expansion from section VIII. D. 
Now we write the d/i integrand in the form 

£xp{K){A,<l>)^J2jp_ E U^i^^^'^) (8-78) 

N ' Xi,...,XjvCA j 

UXj=A 

The rough idea will be that each Xj is a contour Tj or several conglomerated contours 
together with a collar of width M which is thick enough so that the effect of the contour is 
locked inside Xj. To take care of regions without contours there will also be the dominant 
possibility that X is a single cube. 

The outcome of our manipulations will be a polymer activity K{X, 0) whose p— th vari- 
ational derivative with respect to 4> is of size 

^exp{nmj,-Oia)\X\), a = [zll^/^ \X\ > 1 (8.79) 

is the number of particles in a Debye sphere and is very large near the Debye-Hiickel 
limit. The exp(— 0(a) |X|) is the tunneling factor promised in section II. G. |X| is the 
number of unit boxes in X. \\ \\x is a norm involving |0p and integrals of derivatives 
of 4>. The constant k, is small enough so that exp(fi;||0||^) can be integrated with respect 
to d//DH- Since polymers never overlap this is enough for a proof that the expansion of 
section VIII. D is convergent. The non-overlapping of polymers is the reason why they are 
useful for obtaining convergent expansions. 

To understand why this crucial limitation on k cannot be achieved unless the tunneling 
effects represented by contours are taken into consideration consider the one dimensional 
version of the Villain approximation. 
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^o?cos{4>/a) _ ^Rz{4>) 



If we restrict Rz{(t>) to the period [— tto, tto] and note that it is largest at = ott we find, by 
throwing away all but the n — Q term, an upper bound 

ei?3W<el-0^ ^=i_i^<i (8.81) 

Since R^{(j)) is a periodic function the bound holds for all (p. R3 is analytic with more or less 
the same Gaussian bound in the strip (p + iy with \y\/a <^ 1 so 



e^'^'f' (8.82) 

K < \ means that cxp {R■^) and all derivatives can be integrated against the one dimensional 
analogue of d/iDu, namely the Gaussian exp (—0^/2) 



«2iLe«.»W = |0(';-"):fP^° (8.83) 
dxP W 1 if p = ^ ^ 



2n J dxP 1^1 if p = 



as a — > 00. This fails if R3 is instead defined as the error in the single well approximation 
a2cos((/)/a) fti 1 -(/)72 + i?3. 

The tunneling phenomena contained in the Villain approximation arc invisible in per- 
turbation theory but they are the reason why perturbation theory will not converge. The 
appearance of exp(— a|X|) in (8.79) quantifies the physical idea that the tunneling effects 
are very rare in the Debye-Hiickel limit when a is large. This bound (8.79) together with 
the considerations of section VIII. D implies the exponential decay of correlations. 

The first task on the way to (8.78) is to obtain an exact version of (2.84) so that the sum 
over h will factor as in (2.86). Let M be a (large) integer. Given h we choose the translation 
g to minimize SB}i{—ih,ig) subject to the constraint that g{r) = h{r) at every point that is 
distant by more than M — 1 from any contour in h. The price of imposing this constraint 
is that in (8.77) we no longer have S'j:,}i{—ih,ig;i(f>) = as was the case in section II. G. 
Now S'-D}i{—ih, g;i(f)) is an integral of (p and derivatives of (j) times derivatives of g over the 
surface of the region where the constraint holds. These derivatives are of order 0(exp(— M)) 
because the constraint is almost irrelevant since the minimizer without the constraint satisfies 
g{r) — h{r) + 0(exp (— dist(r, F)). We choose M large so that S'^ni—ih, g'-jicf)) is as small 
as the error term Rz{ig + 

With this choice of g wc have an exact version of the additivity property of S-^ni—ih, ig) 
in (2.84): whenever h = hi + h2, 



SDui-iKig) = SBui-ihj^igj) 
7=1,2 



provided every contour in hi is separated by a distance of 2M or more from every contour 
in /i2. See the discussion below (2.85) and note that our constraint forces the cross terms to 
vanish. Under the same conditions we also have 
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S'BK{-ih, g; i(t)) = ^ S'^ai-ikj, gf, i(t)) (8.85) 

3 

Define K{X, i(f)) to vanish if X is the empty set or not a connected union of cubes, 
otherwise 

The second fine apphes when X is a single cube. The constraint X(K) = X is defined 
below. The constraint forces X to hold contours when it is not a single cube and con- 
tours are responsible for the exp(— 0(a)|X|) in (8.79). This definition leads to the desired 
representation 

SA/(Sideal,ASDH,A) = / ci/.DH(0) E M E 11 ^^^^^ (^.87) 

N ' Xi,...,XfjCA j 

UXj=A 

It is convenient to measure distance using 

dist(r, r') = max |rj — (8.88) 

i 

SO that the unit cube centered on the origin is {r : dist(0,r) < 1/2}. X{h) is defined to be 
the open set consisting of all points closer than M to some contour in h. To derive (8.87) 
note that any h occurring in the sum (8.77) determines X{h) which decomposes uniquely 
into disjoint open connected subsets Xj. Each Xj determines a unique hj such that 

X{hj) = Xj] hj = outside A 

hj = h + a constant when restricted to Xj 

h^Y.^j (8.89) 



3 



The first two conditions determine hj and the last is implied by them. Therefore the sum 
over h in (8.77) is equivalent to summing over iV, Xi, . . . , Xjv followed by /ii, . . . , h^- By 
(8.84) and the symmetry 

^R3{Xk,iY,j gj+i<t>) — QR3{^k,igk+i<t>) j'g gg'j 

we can factor the summands. This symmetry is the one wc remarked in (2.75) and it holds 
in this form because for j ^ k Qj = hj is constant on Xj.. Incidentally a fractional charge 
immersed in the system destroys this symmetry. See the remarks on fractional charges in 
section II. G. 



F. Proofs for theorems in chapter Vlll 



Proof of theorems VIII. 4 and VIII. 5: Following (Brydges & Kennedy 1987). Let x(a > r) 
denote a function of r which is 1 when a > r and zero otherwise. Then 
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UMoa{T,s) = J2 drVr^r) 

= / c?r5]K;(^)xK(r,s)>r) (8.91) 

A value of r and a tree T determine a partition of the particles = vertices of T into clusters: 
simply erase all bonds ij e T for which Sij < r. The rest of the tree falls into connected 
subtrees. Each cluster consists of the vertices in each subtree. The bonds ij E T we erased 
are in fact eliminated by the x the above equation. The key point is that if Z A; is another 
bond not in T with k and I is separate clusters that used to be connected by ij then it too 
is eliminated by x because Sij < r implies aki{T,s) < r. Therefore at each value of r we 
have a r dependent partition of the vertices into clusters Ca such that 

E V:,{T)x{'7iAT, s) > r) = 5] 5; yA(r) (8.92) 

ij a ijeCa 

Thus a stability estimate J2 ^iji'^) — —B{t)N is inherited by Un o <t(T, s), because we can 
apply it to each cluster in the last equation and then integrate over r. Theorem VIII.4 
follows easily. 

Theorem VIII. 5: Following (Gopfert & Mack 1981) there is a systematic way to improve 
stability estimates by using numerically calculated constants for small clusters. The simplest 
instance is to note that for clusters of single particles we have no interaction to bound so 
that 

E V{,{T)x{a,,{T, s) > r) > -B{t)N, (8.93) 

ij 

where Ni is the number of particles not in single clusters, 

Ni — \{i : i e some Cj with \Cj\ > 1}\ 

< |{i : Sij > T some ij in T}\ 

< Xisij >r)^2j2 Xisij > r) (8.94) 



Therefore 



Um o o-(r, s) > - / dr 2B{t) V x{sij > r) 
•^0 ij 



(8.95) 



which is (8.18). QED 

Proof of Theorem VIII. 2: We know the case N — 2. We assume N > 2 and make the 
inductive hypothesis that the theorem holds for all cases with less than N vertices. It is 
sufficient to prove that both sides have the same t derivative and agree at t — 0. The t 
derivative of the right hand side is 
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(8.96) 



G&T b€G 

We interchange the sums and obtain where 

Xb^ Yl f d''\'sd''Fo(T{G,s)s,=t (8.97) 

G&T,G3b 

Now we fix h. To use the inductive hypothesis, we reduce by one by rewriting Xf, as a 
sum over all forests G on N — 1 vertices. The N — 1 vertices are obtained from the original 
vertices {1, 2, . . . , A^} by identifying the vertices at either end of b. The possible bonds c are 
viewed as sets of identified bonds on {1,2,..., N}. For example, ii b — 12, then li and 2i 
become the same bond which we would denote by li or 2i. 

s denotes a vector of parameters Sg assigned to bonds c. Given such parameters, we 
determine parameters s by assigning Sc to all bonds c e c and t to bond b. Thus define 
s = j (s) where the c component is given by 

To each forest G 3 b on N vertices we associate a unique forest G/6 on — 1 vertices by 
identifying the vertices in the pair b. Thus 

G/b^{c:ce G} (8.99) 
For any function / of parameters s, 

(^''^V)oj(s) = a'5(/oj)(s) (8.100) 

G:G/b=G 

To understand this, consider 

/ U=S \ / v=s 

Also, for any function / of parameters s and G 3 b, 

j'd''\''sfoa{G,s)s,=t = j' <f'hfo-^oa{G/b,s) (8.101) 

The key point here is that (Tc(G, ^)sb=t depends only on c. 
We return to (8.97). By summing over G with G/b held fixed 



.102) 



G G:G/b=G 
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By (8.101) this equals 

J2 J2 f d^sd^Foioa{G,s) (8.103) 

G G:G/b=G 

which by (8.100), with / = d^F is 

J2 I d^s{d^dbFoi)o&{G,s) (8.104) 
G 

which by the inductive hypothesis is 

{dbF)o-^{tl)=dbF{tV) (8.105) 

Therefore the t derivative of the right hand side of (8.6) is dbF(t1) which is the t 

derivative of the left hand side of (8.6). QED 
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FIGURES 




FIG. 1. A q particle loop 
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FIG. 2. Graphical definitions for bonds 



FIG. 3. A Coulomb point 
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FIG. 4. Resummed Coulomb chains 
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= •-> 

FIG. 5. Multipole-multipole resummed bond 
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FIG. 6. f-bonds with Coulomb chains in parallel 
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FIG. 7. Coulomb chains in parallel 
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FIG. 8. Graphical symbol for right hand side of (5.70) bond 




FIG. 9. Vertex with two propagator lines and one interaction line 
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FIG. 10. The random phase approximation 
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FIG. 11. A first correction to RPA 
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FIG. 12. Examples of contributions to Maxwell-Boltzmann free energy 
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FIG. 13. The inequalities (7.18) 
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